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Executive  Summary 


This  report  documents  the  latest  version  of  the  Stochastic  Lanchester  Air-to-Air 
Campaign  Model  (SLAACM),  developed  by  LMI  for  the  Tactical  Air  Division,  a 
component  of  the  Office  of  the  Secretary  of  Defense,  Program  Analysis  and 
Evaluation  (OSD  PA&E).  During  the  past  year,  optimized  offensive  air  cam¬ 
paigns,  including  suppression  of  enemy  air  defense  and  the  impact  of  electronic 
warfare  and  low  observable  technology,  have  been  added  to  the  basic  air  defense 
SLAACM  to  produce  a  new  “attack”  version  that  includes  both  defense  and  of¬ 
fense  scenarios.  The  “classic”  air  defense  SLAACM  has  been  retained  as  a  sepa¬ 
rate  version. 

SLAACM  is  a  fast  and  flexible  model  designed  for  analysts  who  need  to  consider 
many  combat  scenario  options  quickly  and  who  wish  to  have  indications  of  the 
uncertainties  in  military  outcomes.  SLAACM  models  both  the  defensive  and  of¬ 
fensive  counter-air  battle,  including  order  of  battle  optimization  by  both  attackers 
and  defenders.  The  current  version  can  address  issues  of  fighter  combat  effective¬ 
ness,  defense  battle  management  capabilities,  time-phasing  of  offensive  and  de¬ 
fensive  forces,  use  of  indigenous  defensive  forces,  and  basing  options  for 
deployed  forces.  Unique,  efficient  mathematical  algorithms  and  fast  integer  pro¬ 
gramming  tools  allow  SLAACM  to  generate  optimized  results,  including  statisti¬ 
cal  uncertainties  for  numerous  types  and  large  numbers  of  combatants  for 
campaigns  of  many  days,  with  sufficient  speed  to  allow  statistical  variation  of 
campaign  scenarios  and  deployment  options.  Included  tables  and  charts  display 
the  day-by-day  development  of  the  campaign.  In  addition,  the  Microsoft  Excel, 
Visual  Basic  implementation  of  the  model  allows  wide  flexibility  for  storing  and 
displaying  results  in  user-preferred  formats. 

This  report  includes  descriptions  of  the  fundamental  mathematics  and  analytical 
structure  of  SLAACM  and  users  guides  both  the  classic  and  attack  versions. 

The  authors  wish  to  thank  Mr.  Bradley  Berkson,  the  director  of  OSD  PA&E,  for 
his  support  of  the  development  of  SLAACM. 
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Chapter  1 

Introduction 


The  Stochastic  Lanchester  Air-to-Air  Campaign  Model  (SLAACM)  was  devel¬ 
oped  by  LMI  for  the  DoD  Office  of  the  Director,  Program  Analysis  and  Evalua¬ 
tion,  Tactical  Air  Division  (TACAIR).  This  parsimonious  probabilistic  air 
engagement  and  optimized  campaign  model  is  designed  to  be  a  flexible  tool  for 
analysts  who  need  to  consider  many  cases  quickly  in  the  PC  environment  and 
who  wish  to  have  indications  of  the  uncertainties  in  military  outcomes.  SLAACM 
combines  stochastic  Lanchester  engagement  modeling  with  integer  programming 
campaign  optimization. 

SLAACM  is  an  analytic,  probabilistic  model,  as  compared  to  a  Monte  Carlo 
simulation,  and  provides  complete  and  repeatable  traceability  between  input  and 
output  variables.  The  stochastic  Lanchester  engagement  method  used  in 
SLAACM  was  developed  to  address  shortcomings  in  deterministic  Lanchester 
analysis  when  modeling  campaigns  involving  multiple  few-on-few  engagements 
of  heterogeneous  forces,  i.e.,  involving  many  flights  of  many  types  of  aircraft. 
Stochastic  Lanchester  analysis  can  directly  model  large  engagements  consisting 
of  multiple  few-on-few  engagements  and  provides  full  statistics  of  the  outcomes. 
LMI  has  developed  a  unique,  high-speed  algorithm  for  engagement  calculations 
that  is  key  to  practical  stochastic  Lanchester  campaign  analysis. 

SLAACM  is  implemented  in  Microsoft  Excel  and  Visual  Basic  for  Applications 
(VBA)  and  typically  runs  a  fully  optimized  10-day  campaign  on  a  standard  PC  in 
less  than  2  minutes.  SLAACM ’s  speed  and  flexibility  are  due  to  use  of  the  fast 
engagement  analysis  algorithm,  the  capability  of  modern  PCs  and  the  use  of  either 
a  fast  heuristic,  or  a  fast  commercial  IP  solver  for  optimization. 

SLAACM  outputs  in  benchmarking  tests  compare  closely  with  those  of  some 
large-scale  simulation  models,  but  SLAACM  is  intended  to  complement,  rather 
than  replace,  such  models.  The  basic  SLAACM  is  unclassified.  Analysis  of  classi¬ 
fied  scenarios  is  performed  by  loading  the  SLAACM  on  classified  computers  and 
running  it  with  classified  weapon  system  and  scenario  input  data.  The  defensive 
air  campaign  capabilities  of  SLAACM  been  developed  in  a  boot-strap  fashion 
over  the  past  8  years,  adding  features  in  response  to  analysis  requirements.  The 
“classic”  defensive  air  campaign  version  of  SLAACM  has  been  used  to  model 
and  analyze  several  large  defensive  campaign  scenarios  of  interest  to  DoD. 

During  the  past  year,  we  have  expanded  SLAACM  to  model  offensive  air  cam¬ 
paigns.  We  have  expanded  the  offense  to  include  typical  U.S.  attack  packages, 
including  suppression  of  enemy  air  defense  (SEAD)  escort  aircraft.  We  have 
added  surface-to-air  missiles  (SAMs)  to  the  defense.  The  model  expansion  has 
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been  accomplished  through  a  major  transfonnation  of  the  VBA  code  structure 
from  procedural  to  object-oriented,  allowing  straightforward  inclusion  of  new  of¬ 
fense  and  defense  packages.  The  new  structure  supports  the  ultimate  goal  of  mod¬ 
eling  an  integrated,  optimized,  campaign  using  both  offensive  and  defensive  air 
power  in  order  to  inform  acquisition,  deployment,  and  employment  decisions. 

This  report  describes  the  analytical  and  statistical  methods  used  to  generate 
SLAACM  outputs,  including  the  new  offensive  capabilities.  It  also  gives  detailed 
instructions  for  operating  SLAACM,  with  illustrative  examples. 

SLAACM  Overview 


SLAACM  treats  campaigns  between  two  opponents,  called  Attack  and  Defense.1 
The  Attack  side  uses  fighter-escorted  bombers  to  attack  assets  defended  by  the 
Defense  side’s  fighters  and  SAMs.  In  each  day’s  operations,  the  Attack  side  de¬ 
termines  optimal  allocations  of  its  fighters  and  bombers  to  attack  packages, 
maximizing  a  payoff  function  that  considers  the  bombers’  payload  and  weapon- 
delivery  accuracy,  the  value  to  Attack  of  destroying  Defense  fighters,  and  the 
penalty  to  Attack  of  losing  fighters  and  bombers. 

Defense’s  fighters  respond.  Each  Defense  fighter  can  be  “smart”  or  not.  Smart 
fighters  can  determine  the  makeup  of  oncoming  attack  packages  before  engaging 
and  can  coordinate  their  operations  to  make  optimal  defenses,  maximizing  a  pay¬ 
off  function  that  considers  the  value  to  Defense  of  destroying  Attack  fighters  and 
bombers,  and  the  penalty  to  Defense  of  losing  a  fighter. 

Defense  fighters  that  are  not  smart  encounter  attack  packages  randomly.  An  ad¬ 
justable  parameter  varies  the  effectiveness  of  Defense’s  battle  management  for 
not-smart  fighters,  reducing  the  probability  of  engagements  below  that  implied 
simply  by  the  numbers  of  attacking  packages  and  defending  flights. 

SAMs  are  included  in  the  new  Attack  SLAACM  scenarios,  but  not  in  the  Air 
Defense  SLAACM  scenario.  Some  Attack  SLAACM  scenarios  include  SEAD 
aircraft  to  suppress  the  SAMs. 

In  Attack  SLAACM  scenarios,  each  attack  package  engages  one  target,  which  is 
defended  by  one  SAM  site.  If  not  destroyed  by  SEAD  aircraft,  the  SAM  site  is 
capable  of  firing  a  user-input  number  Nm  of  missiles  at  the  bombers.  SAMs  en¬ 
gage  SEAD  aircraft  and  bombers,  but  not  attacker  escort  fighters.  Defending 
fighters  engage  attacker  escorts  and  bombers,  but  not  SEAD  aircraft. 


1  In  previous  SLAACM  documentation,  friendly  forces  are  identified  as  Blue  and  Green,  and 
enemy  forces  are  Red.  With  the  ability  of  friendly  forces  to  defend  or  attack,  we  have  changed  the 
nomenclature  in  this  report  to  “Attack”  and  “Defense.”  The  previous  version  has  been  called  “Red 
SLAACM,”  “Air  Defense  SLAACM,”  and  “Classic  SLAACM,”  and  the  new  SLAACM  version 
incorporating  attack  has  been  called  “Blue  SLAACM”  and  “Attack  SLAACM.”  Despite  the 
names,  the  new  version  includes  both  defensive  and  offensive  scenarios  and  is  a  superset  of  the 
older  version. 
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Introduction 


SLAACM  does  not  use  iterative  simulation.  Rather,  it  calculates  the  outcome 
probabilities  using  analytic  probability  algorithms,  including  the  full  loss  statistics 
of  each  day’s  combat.  The  mean  losses  are  propagated  and  used,  along  with  any 

replenishment  aircraft,  to  detennine  Attack  and  Defense  orders  of  battle  (OOBs) 

2 

for  the  next  day’s  operations. 

SLAACM  tracks  the  destructive  effectiveness  of  the  bombs  delivered,  the  total 
tons  of  bombs  delivered,  as  well  as  the  tons  of  smart  bombs  (guided  munitions 
with  higher  destructive  capability  per  ton).  SLAACM  also  tracks  which  defense 
aircraft  killed  which  attack  aircraft  and  vice  versa,  and  standard  deviations  of 
Defense  and  Attack  losses. 

SLAACM ’s  worksheets  and  analytic  structure  provide  complete  audit  trails  for 
the  input  parameters  used  in  the  calculations.  The  kill-rate  ratios  used  in  the  calcu¬ 
lations  may  be  obtained  in  several  ways.  The  kill-rate  ratios  may  be  derived  from 
probabilities  of  detection  and  kill  (Pd  -  Pk  probabilities)  of  the  combinations  of 
combatants  involved.  They  may  also  be  derived  from  loss  ratios  found  from  simu¬ 
lations,  including  man-in-loop  simulations,  or  from  data  on  actual  air  combat,  or 
from  combinations  of  all  these  sources.  This  flexibility  of  inputs  is  an  important 
strength  of  SLAACM. 

Obtaining  the  kill-rate  ratios  from  simulation  results  gives  SLAACM  results  a 
basis  of  consistency  with  the  simulations.  (Again,  we  point  out  that  SLAACM  is 
intended  to  complement,  not  supplement,  large-scale  simulations.)  Because  many 
SLAACM  exercises  have  been  based  on  loss  ratios,  the  input  loss  ratios  are  dis¬ 
played  in  the  model.  In  classified  versions  of  SLAACM,  the  loss  ratios  are  based 
on  documented  historical  data,  simulations,  and  military  expertise.  The  analytic 
nature  of  the  model  ensures  that  the  effects  of  parameter  changes  are  consistently 
and  individually  reflected  in  the  output  and  are  not  convolved  with  other  parame¬ 
ter  effects  as  happens  in  iterative  simulations.  These  features  are  quite  helpful  to 
understanding  the  implications  of  different  scenario  options  and  weapons  parame¬ 
ters. 

Report  Organization 


This  report  is  organized  as  follows: 

♦  Chapters  2  through  5  describe  how  SLAACM  works.  Specifically,  those 
chapters,  respectively,  describe  how  the  model  is  used  to 

>  calculate  air-to-air,  SEAD,  and  SAM  engagement  results, 

>  calculate  bomber  effectiveness, 


2  The  term  “day”  here  refers  to  a  single  set  of  campaign  sorties.  There  could  be  more  than  one 
such  set  on  a  single  calendar  day,  or  there  could  be  calendar  day  gaps  between  sets. 
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>  specify  the  Attack  and  Defense  payoff  logic  and  calculate  optimal  at¬ 
tacks  and  defenses,  and 

>  calculate  campaign  results. 

♦  Chapters  6  and  7  are  users  guides  for  operating  Classic  SLAACM  and 
Attack  SLAACM,  respectively. 

The  appendixes  contain  supplemental  information  of  interest  for  combat  model¬ 
ing: 

♦  Appendix  A  discusses  the  mathematics  of  the  fast  engagement  analysis 
algorithm. 

♦  Appendix  B  addresses  the  use  of  NASA-developed  Markov  tools  for  aux¬ 
iliary  engagement  analysis  and  scenario  development. 
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Chapter  2 

How  SLAACM  Works: 

Engagement  Models  and  Calculations 


SLAACM  does  not  simulate;  rather,  it  makes  complete  calculations  of  outcome 
probability  distributions  for  specific  probabilistic  engagement  models.  That  ap¬ 
proach  is  significantly  different  than  the  approach  used  in  other  campaign  models, 
so  it  is  appropriate  to  begin  our  discussion  by  explaining  the  SLAACM  approach. 

Basic  M  vs.  N  Probabilistic  Engagement  Model 


The  heart  of  SLAACM’s  calculations  is  a  set  of  probabilistic  models  of  air-to-air 
engagements.  The  simplest  engagement  model  used  in  SLAACM  is  the  M  vs.  N 
probabilistic  engagement  model.  We  describe  that  model  in  considerable  detail,  to 
clarify  some  fundamental  aspects  of  SLAACM’s  operation. 

In  previous  reports,  defenders  were  friendly  and  designated  as  Blue  or  Green  air¬ 
craft,  while  attackers  were  hostile  and  identified  as  Red  aircraft.  With  the 
modeling  of  offensive  operations  by  friendly  forces,  this  nomenclature  no  longer 
applies;  instead,  whenever  we  deal  with  specific  attack/defense  scenarios,  we  will 
refer  to  Attack  and  Defense  aircraft.  However,  for  discussions  of  basic  engage¬ 
ment  modeling,  we  retain  the  Red  and  Blue  nomenclature. 

In  the  basic  M  vs.  N  probabilistic  model,  M  Blue  aircraft  engage  N  Red  aircraft. 
Throughout  the  engagement,  the  time  between  kills  by  each  Blue  aircraft  is  taken 
to  be  identically  and  independently  distributed  exponentially  with  parameter  kb, 
and  the  time  between  kills  by  each  Red  aircraft  is  assumed  to  be  identically  and 
independently  distributed  exponentially  with  parameter  kr.  Specifically, 

Tr  ~  kre_krTr  [Eq.  2-1] 


and 


Tb  ~  kbe-kbTb  [Eq.  2-2] 

where  means  “distributed  as”  and  xr  and  ib  are,  respectively,  the  time  between 
kills  made  by  each  Red  aircraft  and  the  time  between  kills  made  by  each  Blue  air¬ 
craft. 

If  the  time  between  kills  by  an  engaged  aircraft  is  exponentially  distributed  with 
parameter  k,  then  the  probability  that  the  aircraft  makes  a  kill  in  a  time  interval 
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(t,  t  +  5t)  is  k5t  +  0(8t~),  independent  of  total  elapsed  time,  t.  (This  property  is 
sometimes  called  the  “memoryless”  property  of  exponential  distributions.) 

If  two  such  aircraft  are  present,  then  by  virtue  of  the  assumption  that  their  times 
between  kills  are  statistically  independent,  the  probability  that  the  pair  makes  a 
kill  in  the  interval  (t,  t  +  5t)  is  k5t  +  k5t  +  0(5t2)  =  2k8t  +  0(8t2).  If  m  such  aircraft 
are  present,  the  probability  that  the  set  of  them  makes  a  kill  in  the  interval 
(t,  t  +  8t)  is  mk8t  +  0(8t2). 

If,  at  a  given  time  during  the  engagement,  M  Blue  aircraft  and  N  Red  aircraft  are 
present,  we  will  say  that  the  engagement  is  in  state  (m,  n).  Figure  2-1  diagrams 
the  possible  transitions  that  the  engagement  can  make  into  and  out  of  state  (m,  n). 


Figure  2-1.  Engagement  State  Transition  Diagram 


Now  let  P m,n(t)  be  the  probability  that  the  engagement  is  in  state  (m,  n)  at  time  t. 
Let  us  consider  the  probability  Pm,n(t  +  8t).  The  engagement  can  reach  state  (m,  n) 
at  time  (t  +  8t)  in  just  three  ways: 

♦  The  engagement  was  in  state  (m,  n)  at  time  t,  and  no  kills  happened  dur¬ 
ing  the  interval  (t,  t  +  St). 

♦  The  engagement  was  in  state  (m,  n  +  1)  at  time  t,  and  the  set  of  Blue  air¬ 
craft  made  a  kill  during  the  interval  (t,  t  +  8t). 

♦  The  engagement  was  in  state  (m  +  1,  n)  at  time  t,  and  the  set  of  Red  air¬ 
craft  made  a  kill  during  the  interval  (t,  t  +  8t). 


1  0(5t2)  signifies  “higher  order  terms”  of  order  5t2  and  higher,  which  are  very  small  and  can  be 
ignored. 
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By  the  assumptions  that  all  the  kill  events  are  independent,  these  three  events  are 
statistically  independent.  Accordingly,  the  probability  that  one  of  them  occurs  is 
the  sum  of  their  individual  probabilities,  so  that,  neglecting  terms  of  0(5t“), 

pm,„  (t  +  St)  =  Pm  n  (t)  [l  -  mk  b8t  -  nk  r8t]  +  Pm  n+1  (t)  [[mk  b8t  +  Pm+Kn  (t)nk  r5tJJ  .[Eq.  2-3] 

Subtracting  Pm,n(t)  from  each  side  of  Equation  2-3,  dividing  by  8t,  and  taking  the 
limit  as  St  — >  0 ,  gives 

Pm,n(t)  =  _(mkb  +nkr)Pmn  +mkbPm  n+1  +nkrPm+l  n  .  [Eq.  2-4] 


With  the  result  (Equation  2-4),  we  can  write  down  a  general  initial  value  problem 
describing  the  evolution  of  the  statistics  of  an  M  vs.  N  engagement,  in  which  the 
Blue  forces  break  away  when  they  have  fewer  than  bmm  aircraft,  and  the  Red 
forces  break  away  when  they  have  fewer  than  rmjn  aircraft: 


Pm,n  (mkb  +  nkr)Pmn  +  mkbPm  n+i  +  nkrPm+1  n, 
vpmin  <  m  <  M  and  rmin  <  n  <  N; 

P0,n  =  nkrPbmin,n  /  rmin  <  n  <  N 
Pm,0  —  ^kbPm,rmm'  bmin  —  ^  ^ 

Pm,n  (°)  =  0  except  PM/N  (0)  =  1 


[Eq.  2-5] 


(The  probabilities  that  the  system  has  more  than  M  Blue  aircraft,  or  more  than  N 
Red  aircraft,  are  of  course  always  zero.) 

The  initial  value  problem  (Equation  2-5)  presents  a  system  of  linear  ordinary  dif¬ 
ferential  equations.  The  number  of  equations  is  equal  to 

(M  -  bmin  +  1)(N  -  rm;n  +  1)  +  (M  -  bmin  +1)  +  (N  -  rmin  +1).  The  equations  for  the 
probabilities  of  the  absorbing  boundary  states,  Po,n  and  Pmjo,  decouple  from  those 
for  the  transient  states  (those  whose  probabilities  tend  to  zero  for  long  times), 
however.  Therefore,  Equation  2-5  may  be  treated  by  solving  only  the 
(M  -  bmin  +  1)(N  -  rm;n  +  1)  equations  for  the  transient  states  and  computing  the 
probabilities  of  the  absorbing  boundary  states  by  integration. 
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In  any  case,  Equation  2-5  offers  no  difficulties  to  numerical  solution.  Moreover, 
introducing  the  nondimensional  time  parameter  x  =  krt  allows  us  to  write  Equa¬ 
tion  2-5  in  a  form  that  involves  only  the  nondimensional  parameter  Kill-Rate 
Ratio  =  KRR  =  k  =  kb/kr: 


Pm,„  =  “(mK  +  n)Pm.n  +  mKPrn,„+l  +  nPm+l,n> 

vb„„n  <  m  <  M  and  rmin  <  n  <  N; 

P0,n  =  nPbmin,n  »  hnin  -  11  -  N  [Eq.  2-6] 

Pm,0  =  mKPm.nnin,  t>min  <  HI  <  M 
Pm,„(°)  =  0  eXCePt  PM,n(°)  =  1 


The  state  probabilities  in  Equation  2-6  are  functions  of  x,  not  t,  and  so  are  differ¬ 
ent  functions  from  those  giving  the  state  probabilities  as  functions  of  t.  Also  the 
derivatives  in  Equation  2-6  are  with  respect  to  x,  not  t.  In  the  interest  of  unclut¬ 
tered  notation,  we  will  not  introduce  new  symbols  for  these  new  functions;  rather, 
in  Equation  2-6  and  in  all  following  parts  of  this  chapter,  the  probabilities  are 
functions  of  x,  and  dots  denote  differentiation  with  respect  to  x. 

Writing  the  initial  value  problem  for  the  M  vs.  N  probabilistic  engagement  in 
nondimensional  form  gives  useful  insights  into  the  nature  of  the  problem’s  solu¬ 
tions.  First,  because  the  only  parameter  appearing  in  Equation  2-6  is  k,  we  see  that 
the  probabilities  Pm>n,  which  in  general  depend  on  time  t  and  the  two  parameters  kr 
and  kb,  can  always  be  written  as  functions  of  the  nondimensional  time  x  and  just 
one  parameter,  the  kill-rate  ratio  k. 

Second,  the  fact  that  the  solutions  of  the  initial  value  problem  depend  only  on  k 
and  x  leads  to  the  conclusion  that  the  limiting  values  of  the  probabilities  as  time 
tends  to  infinity  (which  of  course  is  also  the  limit  as  x  tends  to  infinity)  are  func¬ 
tions  only  of  the  single  parameter  k. 

Figure  2-2  illustrates  these  ideas.  It  shows  the  evolution  of  certain  probabilities, 
obtained  by  numerical  integration  for  a  4  vs.  4  engagement  with  kb/kr  =  k  =  1 .46, 
and  in  which  the  combatants  fight  to  annihilation. 
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Figure  2-2.  Evolution  of  State  Probabilities 


— ♦—  P4,4  — P0,3  — A— P0,4  - P3,0  — P4,0 


To  find  the  value  of  one  of  the  probabilities  shown  in  Figure  2-2  at  a  given  time  t 
for  a  given  engagement,  one  would  find  the  value  of  x  corresponding  to  that  t  by 
computing  x  =  kat  and  reading  the  desired  value  from  the  figure.  For  example,  at 
t  =  20  time  units  in  an  engagement  in  which  the  mean  time  for  the  Red  aircraft  to 
make  a  kill  is  180  time  units,  the  value  of  x  is  20/180  ~  0.11;  thus  the  desired 
value  of  P4,4  is  0.4. 

We  will  be  more  interested  in  the  long-time  limiting  values  of  the  probabilities. 
Figure  2-2  indicates  that  for  x  >  2,  the  absorbing  boundary  states  have  essentially 
reached  their  long-time  limiting  values.  The  numerical  solution  confirms  that  for 
x  >  2,  the  total  probability  of  the  absorbing  boundary  states  is  larger  than  95  per¬ 
cent.  The  physical  time  required  to  reach  the  long-time  limit  will  of  course  vary 
with  the  value  of  the  dimensional  parameter  ka.  But  the  values  of  the  probabilities 
depend  only  on  the  single,  nondimensional  parameter  k. 

Relating  Loss  Data  to  Kill-Rate  Ratios 

The  engagement  models  underlying  SLAACM,  like  all  Lanchester-type  models, 
employ  kill  rates  for  the  opposing  forces.  In  detenninistic  Lanchester  models,  the 
parameters  are  deterministic  kill  rates,  i.e.,  kills  per  firer  per  unit  time.  In  stochas¬ 
tic  Lanchester  models  like  SLAACM,  the  parameters  are  the  reciprocals  of  mean 
times  between  kills  by  a  single  firer.  For  simplicity,  we  still  refer  to  parameters  of 
stochastic  Lanchester  models  as  “kill  rates.” 

Observed  engagement  data,  from  simulations  or  from  observations  of  actual  com¬ 
bat,  are  rarely,  if  ever,  rate  data.  Rather,  they  are  typically  loss  ratios  for  specific 
M  vs.  N  engagements.  To  apply  empirical  loss  ratio  data  to  rate-based  models,  we 
need  a  method  to  convert  from  loss  ratio  data  to  kill-rate  ratio  data.  To  do  this,  we 
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make  the  significant  assumption  that  the  long-time  limits  of  SLAACM’s  engage¬ 
ment  models  are  appropriate  models  for  the  observations. 

Thus,  with  the  basic  M  vs.  N  engagement  model,  we  assume  that  the  observed 
engagements  end  only  with  one  side’s  reaching  its  specific  breakaway  condi¬ 
tion — a  certain  number  of  losses,  or  annihilation — rather  than  with  a  clock-time- 
dependent  condition  such  as  fuel  exhaustion,  or  an  event-dependent  condition 
such  as  exhaustion  of  munitions.  More  advanced  engaged  models,  described  in 
Appendix  B,  do  account  for  missile  use. 

As  noted  above,  the  nondimensional  kill-rate  ratio  parameter,  k,  is  the  single 
parameter  of  interest  in  the  basic  engagement  model.  To  infer  a  value  of  k  from 
loss  ratio  data  for  an  M  vs.  N  engagement,  we  find  the  value  of  k  for  which  the 
long-time  limiting  value  of  the  ratio  of  expected  Blue  loss  to  expected  Red  loss  in 
the  M  vs.  N  probabilistic  model  that  models  the  engagement  is  equal  to  the  ob¬ 
served  loss  ratio.  In  this  way,  we  infer  a  nondimensional  quantity,  k,  from  a 
nondimensional  quantity,  the  observed  loss  ratio. 

Suppose,  for  example,  that  the  available  data  are  for  a  2  vs.  2  engagement  to  anni¬ 
hilation  and  that  they  show  a  loss  ratio  of  1  Blue  to  10  Reds.  Calculations  show 
that,  for  the  basic  2  vs.  2  probabilistic  engagement  model,  a  kill-rate  ratio  k  of 
8.28  makes  the  ratio  of  expected  Blue  loss  to  expected  Red  loss  equal  to  1 : 10. 2 

Although  we  usually  use  a  finite-step  iterative  method  to  compute  long-time  lim¬ 
its  directly,  it  is  of  course  possible  to  find  long-time  limits  with  time-based 
methods,  such  as  numerical  integration.  In  time-based  Markov  calculations,  the 
time  required  to  reach  the  long-time  limit  will  depend  on  the  value  chosen  for  kr, 
where  kb  is  equal  to  k  times  kr.  For  example,  analysis  shows  that  the  kill-rate  ratio 
of  8.28  results  in  a  1:10  loss  ratio  in  approximately  0.22  time  units  for  kr=  1,  in 
roughly  2.2  time  units  iterations  for  kr=  0. 1,  and  in  about  22  time  units  for 
kr  =  0.01.  The  0.2,  2,  and  22  time  units  represent  adequate  model  time  to  bring  the 
solution  to  the  long-time  limit,  for  their  corresponding  kill  rates.  Thus,  in  a  time- 
unit-based  model,  selection  of  the  individual  rates  is  arbitrary  as  long  as  the  ratio 
k  is  preserved  and  the  calculations  continue  to  a  time  on  the  order  of  twice  the  in¬ 
verse  of  the  sum  of  the  parameters  kb  and  kr.  (Appendix  A  describes  the  general 
analytic  method  that  produces  these  numerical  results.) 

Another  major  assumption  in  our  engagement  modeling  is  that  the  kill-rate  ratio  is 
independent  of  the  Blue:Red  aircraft  ratio.  We  assume  that  kill-rate  ratios  derived 
from  data  on  4  vs.  8  engagements  can  be  applied  directly  to  4  vs.  4  and  other  en¬ 
gagement  combinations.  A  corollary  assumption  is  that  kill-rate  ratios  are 
independent  of  the  state  of  an  engagement.  These  are  significant  and  powerful 


2  The  kill-rate  ratio  reflects  the  performance  of  individual  Blue  vs.  Red  pairs  and  is  numeri¬ 
cally  equal  to  the  loss  ratio  for  a  1  vs.  1  engagement.  A  10:1  kill-rate  ratio  would  be  required  to 
produce  a  1:10  loss  ratio  for  a  1  vs.  1  engagement,  but  only  an  8.28:1  kill-rate  ratio  is  necessary 
for  the  same  loss  ratio  in  a  2  vs.  2  engagement,  because  all  available  aircraft  are  contributing  to 
each  kill. 
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simplifying  assumptions,  which  we  are  always  eager  to  test  with  new  combat 
data. 

Available  data  from  TAC  BRAWLER  engagement  simulations  give  expected  val¬ 
ues  of  Red  losses  and  Blue  losses.  These  data,  some  of  them  classified,  are 
available  for  several  Red  and  Blue  aircraft  types  of  interest.  Government  subject 
matter  experts  have  scaled  these  results  to  model  additional  aircraft  combinations. 
A  SLAACM  utility  uses  bisection  iteration  and  the  long-time  limit  algorithms 
(discussed  below  and  in  Appendix  A)  to  calculate  the  values  of  k  from  the  loss 
ratios  for  each  of  the  Red  vs.  Blue  combinations.  The  calculations  are  repeated 
until  successive  iterations  agree  to  four  significant  figures. 

Relating  Probabilities  of  Detection  and  Probabilities  of  Kill 
Data  to  Kill-Rate  Ratios 

Air  combat  performance  data  are  often  expressed  in  terms  of  probabilities  of  de¬ 
tection,  Pd,  and  probabilities  of  kill,  Pk.  These  can  be  used  directly  to  calculate 
loss  ratios  and  KRRs.  Figure  2-3  diagrams  a  model  that  can  be  used  for  this  pur¬ 
pose. 


Figure  2-3.  Pd,  Pk  Engagement  Diagram 


The  model  described  is  a  1  Blue  vs.  1  Red  engagement  in  which  each  aircraft  can 
detect  the  opponent  and  fire  one  missile.  Each  side  has  an  initial  probability  of 
detection  plus  a  separate  probability  of  detecting  the  opponent  after  the  opponent 
has  fired  his  missile.  Outcomes  include  Blue  kills  Red,  Red  kills  Blue,  each  kill 

3  Loss  ratios  equal  kill-rate  ratios  for  1  vs.  1  engagements. 


2-7 


the  other,  neither  detect,  and  neither  kill.  This  model  can  be  easily  programmed  in 
Microsoft  Excel  to  develop  KRR  values  for  SLAACM. 

Table  2-1  contains  the  state  and  state  transition  infonnation  for  the  model. 

Table  2-1.  Pd,  Pk  Transition  Probability  Table 


State  transition 

From 

To 

Transition  probability 

Initial  detection 

Blue  detects,  Red  does  not 

(0,0, 0,0) 

(1,0, 0,0) 

Pbd  *  Qrd 

Red  detects,  Blue  does  not 

(0,0, 0,0) 

(0,0, 1,0) 

Prd  *  Qbd 

Both  Red  and  Blue  detect 

(0,0, 0,0) 

(1,0, 0,0) 

0.5Pbd  *  Prd 

(0,0, 1,0) 

0.5Pbd*Prd 

Neither  Blue  nor  Red  detect 

(0,0, 0,0) 

No  Eng’t 

Qbd  *  Qrd 

Consequences:  Blue  detects,  Red  does  not 

Blue  kills  Red 

(1,0, 0,0) 

BW 

Pkb 

Blue  misses,  Red  detects  and  kills  Blue 

(1,0, 0,0) 

RW 

Qkb  *  Prr  *Pkr 

Blue  misses,  Red  detects  and  misses 

(1,0, 0,0) 

NL 

Qkb  *  Prr  *  Qkr 

Blue  misses,  Red  does  not  detect 

(1,0, 0,0) 

NL 

Qkb  *  Qrr 

Consequences:  Red  detects,  Blue  does  not 

Red  kills  Blue 

(0,0, 1,0) 

RW 

Pkr 

Red  misses,  Blue  detects  and  kills  Red 

(0,0, 1,0) 

BW 

Qkr  *  Prb  *  Pkb 

Red  misses,  Blue  detects  and  misses 

(0,0, 1,0) 

NL 

Qkr  *  Prb  *  Qkb 

Red  misses,  Blue  does  not  detect 

(0,0, 1,0) 

NL 

Qkr  *  Qrb 

Consequences:  Both  Blue  and  Red  detect 

Blue  kills,  Red  misses 

(1,0, 1,0) 

BW 

Pkb  *  Qkr 

Red  kills,  Blue  misses 

(1,0, 1,0) 

RW 

Pkr*  Qkb 

Neither  Blue  nor  Red  kill 

(1,0, 1,0) 

NL 

Qkb  *  Qkb 

The  following  definitions  apply  to  Table  2-1: 

BW  =  Blue  win  state 

RW  =  Red  win  state 

NL  =  No-loss  state 

Pbw  =  Probability  of  Blue  win 

Prw  =  Probability  of  Red  win 

Pnl  =  Probability  of  no  loss 

Px  =  Probability  of  event  x 

Qx  =  (1-Px)  Probability  of  failure  of  event  x 

Pdb  =  Probability  that  Blue  detects  Red 
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Pdr  =  Probability  the  Red  detects  Blue 
Pkb  =  Probability  that  Blue  kills  Red 
Pkr  =  Probability  that  Red  kills  Blue 

Pbr  =  Probability  that  Blue  detects  Red  after  Red  fires  and  misses 
Prr  =  Probability  that  Red  detects  Blue  after  Blue  fires  and  misses  Red. 

The  model  of  Figure  2-1  and  Table  2-1  avoids  the  unrealistic  outcome  of  simulta¬ 
neous  destruction  of  both  adversaries  by  immediately  transitioning  the  “both 
detect”  case  to  either  the  Blue  detect  and  shoot  or  the  Red  detect  and  shoot  case, 
with  probability  one-half  for  each. 

Now  let  us  generate  KRR  values  from  the  results  of  the  model.  Please  recall  the 
equations  for  a  1  vs.  1  stochastic  Lanchester  equation, 

Pn  =-(kb  +  kr)Pu 

Pio=kbPn  Eq.  2-7] 

P0i  =  krPn 


which,  with  the  initial  condition  of  Pio(O)  =  P(oi(0)  =  0;  Pn(0)  =  1,  leads  to 
Pn  =  exp(kb  +  kr)t 

p,o  =  .  ukb,  t1  -  exp(kb  +  kr)t]  [Eq.  2-8] 

kb  +  kr 

P01  =  — — —  [l  -  exp(kb  +  kr)t] 

01  kb  +  kr  J 

These  equations  imply 


Pio  = 


kb 


kb  +  kr 


[l-Pn]and  P01  = 


kr 


kb  +  kr 


[i-p.,] 


[Eq.  2-9] 


Thus  for  the  stochastic  Lanchester  model,  at  time  t  the  probability  of  no  loss  is 
P ii,  the  probability  of  a  Blue  win  is  Pio,  and  the  probability  of  a  Red  win  is  Poi. 

For  the  Pd~Pk  model,  the  corresponding  probabilities  are  Pnl,  Pbw,  and  Prw-4  To 
infer  KRR  =  kb/kr,  we  may  evaluate  Pu  as  equal  to  PNl,  and  then  solve 


P  =  P 

1  10  iBW 


kb 


kb  +  kr 


M„]= 


KRR 
1  +  KRR 


Mu] 


[Eq.  2-10] 


for  KRR. 


4  One  may  choose  to  divide  the  Pd  -  Pk  probabilities  PBW,  Prw-  and  Pnl  by  1  -  QdrQdb,  tak¬ 
ing  the  position  that  the  stochastic  Lanchester  model  presumes  contact  between  the  adversaries. 
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Table  2-2  gives  two  examples  of  Pd-Pk  data  and  the  resulting  values  of  kill-rate 
ratio. 


Table  2-2.  Parameter  Values  and  Kill-Rate  Ratios 


Parameter 

Value,  Case  1 

Value,  Case  2 

Pbd 

0.8 

0.95 

Prd 

0.4 

0.1 

Pbk 

0.7 

0.7 

Prk 

0.6 

0.3 

Pbr 

0.9 

0.9 

Prr 

0.8 

0.8 

KRR 

2.2 

27.0 

Calculation  Method  for  the  Long-Time  Limit 

If  the  engagement  continues  to  completion — that  is,  until  the  probabilities  that  the 
system  is  not  in  one  of  the  absorbing  boundary  states  (0,  n)  or  (m,  0)  are  negligi¬ 
bly  small — then  the  probabilities  of  these  states  are  functions  only  of  k  and  the 
initial  values  M  and  N.  These  conditions  allow  significant  simplification  in  the 
calculation  of  state  probabilities.  One  can  obtain  values  of  the  long-time  limits  of 
the  probabilities  of  absorbing  boundary  states  by  a  simple,  finite-step  iteration. 
This  speeds  up  the  calculations  tremendously,  compared  to  time-based  Markov 
methods,  and  is  a  key  to  SLAACM  utility. 

An  Example 

As  an  example  of  the  use  of  the  initial  value  problem  (Equation  2-5),  we  consider 
4  Blue  aircraft  engaging  4  Red  aircraft.  Both  sides  break  away  on  sustaining  2 
losses,  that  is,  bmm  =  rmin  =  3. 

Figure  2-4  illustrates  the  states  in  the  engagement  and  the  transitions  between 
them.  In  that  figure,  the  absorbing  boundary  states  are  labeled  as  having  0  aircraft 
for  the  side  that  has  broken  away. 
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Figure  2-4.  Engagement  Diagram 

Red  Win  States 


3kr 


Blue  Win  States 


For  this  case,  the  general  initial  boundary  value  problem  (Equation  2-6  )  takes  this 
specific  form: 


P4/4  =  -(4k  +  4)P4  4 
P4  3  =  -(4k  +  3)P4  3  +  4kP4  4 
P3,4  =  -(3k  +  4)P43  +4Pm 
-^3,3  =  -(3k  +  3)P3,3  +3P4/3  +3kP34 

P4,o=4kP43  [Eq.  2-11] 

^3,0  =  3kP3  3 
P0,3  =  ^3,3 

Po,4  =  4P3,4 

Pm,n  (0)  =  0,  except  P4  4  (0)  =  1 

where  k  =  kb/kr  and  the  (0)  states  correspond  to  the  break  condition  of  two  air¬ 
craft. 

The  Red  side  wins  the  engagement  if  all  the  Blues  are  gone,  that  is,  if  the  en¬ 
gagement  is  in  a  state  (0,  n),  rmjn  <  n  <  N.  The  Blue  side  wins  if  the  engagement  is 
in  a  state  (m,  0),  bmin  <  m  <  M. 

Note  that  the  first  four  equations  in  Equation  2-11  describe  the  evolution  of  the 
probabilities  of  transient  states,  whose  probabilities  will  tend  to  zero  at  long  times. 
The  remaining  four  equations  describe  the  variation  of  the  absorbing  boundary 
states,  or  outcome  states.  At  long  times,  the  system  will  be  found  in  one  of  these 
states. 
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The  equations  for  the  transient  states  decouple  from  those  of  the  outcome  states. 
That  is,  the  equations  involve  only  the  four  transient-state  probabilities.  The  evo¬ 
lution  of  the  transient  states  can  thus  be  detennined  independently  of  the  outcome 
states. 

The  outcome  states’  evolution  is,  by  contrast,  completely  determined  by  the  tran¬ 
sient  states  and  the  initial  conditions.  For  example,  given  P4,3(t),  the  entire  history 
of  P4,o(t)  follows  from 

} P4,0(x)dx  =  P4,0(t)  -  P4o(0)  =  P4,0(t)  =  4k| P4J(x)dx  [Eq.  2-12] 

0  0 

where  we  have  used  the  initial  condition  P4,o(0)  =  0. 

Letting  t  — >  oo  in  Equation  2-12  we  see  that  the  long-time  limiting  value  of  P4,o(t) 
is  given  by 


lim  P4  0(t)  =  4KjP4  3(x)dx.  [Eq.  2-13] 


We  now  show  how  to  exploit  characteristics  of  the  engagement  we  are  consider¬ 
ing,  shown  by  Figure  2-4  and  reflected  in  properties  of  Equation  2-1 1,  to  allow 
calculation  of  long-time  limiting  probabilities  of  outcome  states  with  a  simple, 
finite  iterative  scheme.  This  engagement  is  an  example  of  what  has  been  called  a 
“pure  death”  process:  the  states  always  transition  to  states  with  fewer  participants 
Reflecting  this,  Figure  2-4  is  an  “acyclic”  graph.  That  is,  the  graph  gives  no  path 
by  which  to  return  to  a  previously  occupied  state. 

A  consequence  of  the  acyclic  character  of  Figure  2-4  is  found  in  the  structure  of 
Equation  2-11.  Defining  the  vector  x  as  (P44,  P43,  P34,  P33),  we  see  that  the  four 
equations  for  the  transient  states  can  be  written  as 


x  =  Ax 


[Eq.  2-14] 


where  the  matrix  A  is  given  by 


(-(4k +  4) 


A  = 


v 


4k 

4 

0 


0  0  0  ' 

-(4k +  3)  0  0 

0  -(3k +  4)  0 

3  3k  -(3k  +  3)v 


[Eq.  2-15] 


The  acyclic  character  of  Figure  2-4  is  reflected  in  the  lower-triangular  character  of 
the  matrix  A.  This  lower-triangular  feature  allows  calculation  of  the  integrals 
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from  0  to  go  of  each  of  the  transient-state  probabilities  with  a  simple  iterative 
scheme.  Defining  the  vector  x  by 


x  =  J  x(t)  dt ,  [Eq.  2-16] 

o 


we  see,  on  integrating  Equation  2-14  from  t  =  0  to  t  =  oo  and  remembering  that 
P4,4(0)  =  1  and  all  other  transient-state  probabilities  are  zero  at  t  =  0,  that 

-  e,  =  Ax  [Eq.  2-17] 


where  ei  is  the  vector  (1,  0,  0,  0).  Thus  x  is  detennined  by  the  solution  of  a  sys¬ 
tem  of  linear  algebraic  equations.  Moreover,  by  virtue  of  the  lower- triangular 
character  of  A,  that  system  is  easily  solved  with  an  obvious  iterative  scheme.  On 
writing  out  Equation  2-17  in  detail,  we  find 


-1  =  -(4k  +  4)xj 
0  =  4kX[  -(4k  +  3)x2 
0  =  4xj  -  (3k  +  4)x3 
0  =  3x2  +  3x3  -  (3k  +  3)x4 


[Eq.  2-18] 


To  solve  those  equations,  one  may  obtain  x,  from  the  first  equation,  then  x2  and 
x3  follow  from  the  second  and  third  equations,  and  finally,  x4  follows  from  the 

last  equation.  The  lower-triangular  character  of  the  matrix  A  guarantees  that  the 
nth  equation  in  the  system  of  Equation  2-18  involves  xn  and  does  not  involve  X: 

for  any  j  >  n.  Another  way  to  view  the  benefits  of  the  lower-triangular  character 
of  the  matrix  A  is  to  observe  that  this  makes  Equation  2-17  have  the  fonn  of  a 
system  of  equations  after  the  forward  course  of  Gaussian  elimination  has  been 
performed,  so  that  only  the  rapid  back  course  remains  to  be  done. 

Why  were  we  justified  in  integrating  the  x  from  0  to  co?  Texts  on  differential 
equations5  show  that  solutions  of  equations  such  as  Equation  2-14  are  always  lin¬ 
ear  combinations  of  functions  exp(kjt),  where  k:l  is  an  eigenvalue  of  the  matrix  A, 
when  these  eigenvalues  are  distinct.  The  eigenvalues  of  a  lower-triangular  matrix 
are  the  elements  on  the  main  diagonal.  Inspection  of  Equation  2-15  shows  that 
these  elements  are  distinct,  and  negative.  Thus  solutions  of  Equation  2-11  de¬ 
crease  exponentially  with  time,  and  the  integrals  we  use  exist. 

In  Appendix  A,  we  show  that  the  approach  to  calculating  long-time  limiting  prob¬ 
abilities  illustrated  here  can  be  applied  to  many  engagement  models.  We  use  this 
method,  which  from  now  on  we  will  call  “the  method  of  Appendix  A,”  frequently 
in  the  present  SLAACM  and  in  planned  extensions  to  it. 


5  W.E.  Boyce  and  R.C.  DiPrima,  Elementary’  Differential  Equations ,  Third  Edition  (New 
York:  Wiley,  1976),  Sections  7-6  through  7-9. 
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For  SLAACM,  the  key  result  from  engagement  calculations  is  the  discrete  bivari¬ 
ate  loss  distribution.  This  gives  the  loss  statistics  that  are  propagated  day-by-day 
in  a  campaign.  Figure  2-5  shows  an  example  of  the  discrete  bivariate  loss  distribu¬ 
tion,  for  a  4  vs.  4  engagement  in  which  both  sides  fight  to  annihilation  (breakpoint 
=  0),  and  k  =  1.46.  In  the  figure,  blue  bars  indicate  Blue  win  states  (states  with 
4  Red  losses),  and  red  bars  indicate  Red  win  states  (states  with  4  Blue  losses).  The 
figure  shows  that,  as  should  be  the  case  when  k  >  1,  the  Blues  are  likely  to  win 
the  engagement  and  the  total  probability  of  the  blue  bars  is  distinctly  larger  than 
the  total  probability  of  the  red  bars.  There  is  significant  dispersion  in  both  Red 
and  Blue  losses. 


Figure  2-5.  Bivariate  Loss  Distribution 


Blue  Loss 


Figure  2-6  shows  the  marginal  distributions  of  Blue  and  of  Red  losses  and  the  to¬ 
tal  win  probabilities.  The  marginal  distributions  show  that  the  marginal 
probabilities  are  not  necessarily  monotone.  They  also  show  that,  in  this  case,  the 
probability  of  4  losses  is  the  largest  Blue  probability. 

Figure  2-6.  Marginal  Loss  Distributions  with  KRR=1.46 


Blue  and  Red  Losses 

4v4  Blue  :  Red,  kb/kr  =  1.46,  Fight  to  Death 
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Such  a  “toe-up”  at  the  greatest  possible  loss  may  seem  counterintuitive.  In  fact, 
that  feature  is  certain  to  happen  for  kill-rate  ratios  k  sufficiently  close  to  1.  To  see 
this,  consider  a  4  vs.  4  engagement  to  annihilation,  with  k  =  1 .  The  probability  of 
a  Blue  win  must  equal  the  probability  of  a  Red  win,  and  these  probabilities  must 
both  equal  0.5.  The  probability  of  a  Red  win  is  the  probability  of  4  Blue  losses,  so 
the  probability  of  4  Blue  losses  must  be  0.5. 

The  sum  of  the  marginal  probabilities  of  Blue  losses  must  be  1.  Therefore,  none 
of  the  marginal  probabilities  of  0,  1,  2,  or  3  Blue  losses  can  exceed  0.5.  If,  as  in 
fact  is  the  case,  all  the  marginal  probabilities  of  0,  1,  2,  and  3  Blue  losses  are  posi¬ 
tive,  each  of  them  must  be  less  than  0.5,  so  that,  in  this  case  of  equal  strength 
forces  (k  =  1),  the  marginal  probability  of  4  Blue  losses  must  be  the  largest  mar¬ 
ginal  loss  probability. 


The  marginal  loss  probabilities  are  continuous  functions  of  k  (they  are  in  fact 
rational  functions  of  k,  as  may  be  seen  by  applying  the  method  of  Appendix  A), 
so  that,  as  k  increases  from  1,  at  least  for  a  while,  the  most  probable  Blue  loss 
must  be  4.  Figure  2-7  shows  that  the  Blue  toe-up  is  negligible  for  k  =  5. 


Figure  2-7.  Marginal  Loss  Distribution 
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Scalability  and  Coping  with  Large- 
Dimensioned  Problems 


Blue  and  Red  Loss  Probabilities 
4\4  Blue  :  Red,  kb/kr  =  5,  Fight  to  Death 


The  examples  illustrate  that  initial  value  problems  (Equations  2-5  and  2-6)  are 
tractable,  and  that  their  solutions  give  results  that  military  planners  may  find  use¬ 
ful.  These  examples  also  show  why  computational  probabilistic  models  have  not 
been  widely  used:  Even  moderate  values  of  M  and  N  generate  a  lot  of  equations. 
For  example,  an  engagement  of  4  Defense  aircraft  with  8  Attack  aircraft,  not  at  all 
uncommon,  leads  to  a  system  of  44  differential  equations.  The  more  sophisticated 
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models  that  we  consider  in  Appendix  B  can  lead  to  systems  of  thousands  of  equa¬ 
tions.  Computational  probabilistic  models  have  significant  scalability  concerns. 

Until  fairly  recently,  it  was  not  practical  to  develop  models  to  deal  with  such  large 
state  spaces.  Advances  on  four  fronts  have,  however,  changed  that.  First,  the 
memory  and  speed  of  widely  available  PCs  now  match  or  exceed  those  of  the 
supercomputers  of  only  a  decade  or  so  ago.  Second,  recently  available  applica¬ 
tions  packages  make  it  possible  to  build  and  analyze  large-dimensioned 
probabilistic  models.  In  the  development  of  SLAACM,  and  for  offline  analyses, 
we  have  used  Markov  analysis  tools  developed  by  NASA  for  analysis  of  high- 
reliability  electronics  to  set  up  and  solve  the  large-dimensioned  initial  value  prob¬ 
lems  to  which  computational  probabilistic  models  often  lead.  Appendix  B 
discusses  how  we  have  applied  these  tools  to  air  combat  analysis.  Third,  the  high¬ 
speed  algorithm  discussed  above  allows  us  to  calculate  engagement  results  ex¬ 
tremely  fast.  Fourth,  commercial  packages  are  now  available  to  solve  the  many- 
variable  integer  programming  problems,  to  which  computational  probabilistic 
models  lead. 

Composite  Engagements 


In  the  SLAACM  air-defense  scenario,  engagements  are  always  composite  en¬ 
gagements  of  a  particular  kind,  in  which  a  flight  of  4  defending  fighters  engages 
an  attack  package  made  up  of  4  advanced  escorts,  4  close  escorts,  and  4  bombers. 
We  discuss  that  scenario  in  detail  here  and  then  introduce  the  additional  engage¬ 
ments  we  have  added  for  the  attack  scenarios. 

The  defenders  first  engage  the  advanced  escorts.  If  the  defenders  win  that  en¬ 
gagement,  the  surviving  defenders  engage  the  close  escorts.  If  the  defenders  again 
prevail,  defenders  who  survive  the  engagement  with  close  escorts  engage  the 
bombers.  Results  of  this  final  engagement  detennine  the  distribution  of  the  num¬ 
ber  of  bombers  “leaking”  through  the  defense.  Figure  2-8  illustrates  the  opposing 
forces  for  this  composite  engagement. 

Figure  2-8.  Opposing  Forces  in  SLAACM  Composite  Engagement 


) 

) 

Sixteen  final  outcome  or  absorbing  states  are  possible  when  both  sides  fight  to 
annihilation.  The  outcome  states  may  be  described  by  the  ordered  quadruple  (A, 
C,  B,  D),  where  A  is  the  number  of  surviving  attack  advanced  escorts,  C  is  the 
number  of  attack  close  escorts  that  survive,  B  is  the  number  of  surviving  attack 
bombers,  and  D  is  the  number  of  surviving  defense  fighters.  Table  2-3  lists  the 
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and  Bombers  ( 
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16  possible  absorbing  quadruples  for  the  SLAACM  composite  engagement,  and 
Table  2-4  shows  the  64  total  transition  and  absorbing  states  possible  for  the 
scenario. 


Table  2-3.  Absorbing  States,  SLAACM  Composite  Engagement 


A 

c 

B 

D 

Attack  AE 

Attack  CE 

Attack  B 

Defense 

4 

4 

4 

0 

3 

4 

4 

0 

2 

4 

4 

0 

1 

4 

4 

0 

0 

4 

4 

0 

0 

3 

4 

0 

0 

2 

4 

0 

0 

1 

4 

0 

0 

0 

4 

0 

0 

0 

3 

0 

0 

0 

2 

0 

0 

0 

1 

0 

0 

0 

0 

4 

0 

0 

0 

3 

0 

0 

0 

2 

0 

0 

0 

1 

Table  2-4.  Transition  and  Absorbing  States, 
SLAACM  Composite  Engagement 


A 

B 

C 

D 

Attack  AE 

Attack  CE 

Attack  B 

Defense 

4 

4 

4 

4 

4 

4 

4 

3 

4 

4 

4 

2 

4 

4 

4 

1 

4 

4 

4 

0 

3 

4 

4 

4 

3 

4 

4 

3 

3 

4 

4 

2 

3 

4 

4 

1 

3 

4 

4 

0 

2 

4 

4 

4 

2 

4 

4 

3 

2 

4 

4 

2 

2 

4 

4 

1 

2 

4 

4 

0 
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Table  2-4.  Transition  and  Absorbing  States, 
SLAACM  Composite  Engagement 


A 

B 

c 

D 

Attack  AE 

Attack  CE 

Attack  B 

Defense 

1 

4 

4 

4 

1 

4 

4 

3 

1 

4 

4 

2 

1 

4 

4 

1 

1 

4 

4 

0 

0 

4 

4 

4 

0 

4 

4 

3 

0 

4 

4 

2 

0 

4 

4 

1 

0 

4 

4 

0 

0 

3 

4 

4 

0 

3 

4 

3 

0 

3 

4 

2 

0 

3 

4 

1 

0 

3 

4 

0 

0 

2 

4 

4 

0 

2 

4 

3 

0 

2 

4 

2 

0 

2 

4 

1 

0 

2 

4 

0 

0 

1 

4 

4 

0 

1 

4 

3 

0 

1 

4 

2 

0 

1 

4 

1 

0 

1 
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0 
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4 
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0 
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4 
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0 
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4 

2 

0 

0 

4 

1 

0 

0 

4 

0 

0 

0 

3 

4 

0 

0 

3 

3 

0 

0 

3 

2 

0 

0 

3 

1 

0 

0 

3 

0 

0 
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2 

4 

0 
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2 

3 

0 

0 

2 

2 

0 

0 

2 

1 

0 

0 

2 
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Table  2-4.  Transition  and  Absorbing  States, 
SLAACM  Composite  Engagement 


A 

B 

C 

D 

Attack  AE 

Attack  CE 

Attack  B 

Defense 

0 

0 

1 

4 

0 

0 

1 

3 

0 

0 

1 

2 

0 

0 

1 

1 

0 

0 

1 

0 

0 

0 

0 

4 

0 

0 

0 

3 

0 

0 

0 

2 

0 

0 

0 

1 

The  first  engagement  is  4  defenders  vs.  4  advanced  escorts.  Eight  outcome  states 
are  possible  from  this  engagement:  four  Attack  win  states — (4,  4,  4,  0),  (3,  4,  4, 

0),  (2,  4,  4,  0),  and  (1,  4,  4,  0) — in  which  the  advanced  escorts  defeat  the  four  de¬ 
fenders  with  0,  1,  2,  or  3  losses.  Four  Defense  win  outcomes  are  possible: 

(0,  4,  4,  4),  (0,  4,  4,  3),  (0,  4,  4,  2),  and  (0,  4,  4,  1),  when  Defense  defeats  the  ad¬ 
vanced  escorts  with  0,  1,  2,  or  3  losses.  We  use  the  method  of  Appendix  A  to 
determine  the  probabilities  of  these  states.  If  the  advanced  escorts  win  this  en¬ 
gagement,  the  composite  engagement  ends. 

The  Defense  win  states  cause  the  engagement  to  continue  into  engagements  be¬ 
tween  1,  2,  3,  or  4  defenders  and  the  4  close  escorts.  Four  possible  outcome  states 
are  Attack  win  states:  (0,  4,  4,  0),  (0,  3,  4,  0),  (0,  2,  4,  0),  and  (0,  1,  4,  0),  corre¬ 
sponding  to  close  escort  losses  of  0,  1,  2,  or  3  aircraft. 

The  engagement  between  1  defender  and  the  4  close  escorts  has  one  Defense  win 
state:  (0,  0,  4,  1).  The  engagement  between  2  defenders  and  4  close  escorts  has 
two  Defense  win  states:  (0,  0,  4,  1)  and  (0,  0,  4,  2).  The  engagement  between  3 
defenders  and  4  close  escorts  has  three  Defense  win  states — (0,  0,  4,  1),  (0,  0,  4, 
2),  and  (0,  0,  4,  3) — and  the  engagement  between  4  defenders  and  4  close  escorts 
has  four  Defense  win  states:  (0,  0,  4,  1),  (0,  0,  4,  2),  (0,  0,  4,  3),  and  (0,  0,  4,  4). 
Again,  we  use  the  method  of  Appendix  A  to  evaluate  the  probabilities  of  all  these 
outcomes. 

Thus  the  engagements  with  close  escorts  lead  to  four  Attack  win  states  and  to  four 
cases  of  engagements  between  defenders  and  bombers:  1  vs.  4,  2  vs.  4,  3  vs.  4, 
and  4  vs.  4.  These  engagements  lead  to  the  Attack  win  states  (0,  0,  1,  0),  (0,  0,  2, 
0),  (0,  0,  3,  0),  and  (0,  0,  4,  0)  and  to  the  Defense  win  states  (0,  0,  0,  1),  (0,  0,  0, 

2),  (0,  0,  0,  3),  and  (0,  0,  0,  4).  Using  the  method  of  Appendix  A  to  evaluate  the 
probabilities  of  the  outcome  states  in  the  defenders  vs.  bomber  engagements  com¬ 
pletes  determination  of  the  probabilities  of  all  16  outcome  states  in  Table  2-3. 
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SLAACM  contains  individual  kill-rate  ratios  for  each  Defense- Attack  aircraft 
type-pair  in  the  Defense  and  Attack  inventories.  SLAACM  uses  these  to  calculate 
probabilities  for  each  of  the  16  possible  outcomes  for  all  the  feasible  4  vs.  4+4+4 
combinations.  In  a  campaign.  Attack  and  Defense  use  the  resulting  probabilities 
to  construct  optimal  attack  packages  and  optimal  defenses. 

Modifications  for  Attack  SLAACM 

New  Scenarios 

New  engagement  scenarios  added  to  SLAACM  are  intended  to  represent  potential 
U.S.  attack  package  options  using  both  normal  and  low-observable  aircraft.. 

These  are  shown  in  Table  2-5. 

Table  2-5.  New  SLAACM  Scenarios 


Attack 

Defense 

4  escorts  +  2  bombers 

8,  4,  or  2  defenders  +  1  SAM  with  3  missiles 

4  escorts  +  2  bombers  +  2  SEADs 

8,  4,  or  2  defenders  +  1  SAM  with  3  missiles 

2  escorts  +  1  bomber 

8,  4,  or  2  defenders  +  1  SAM  with  3  missiles 

2  escorts  +  1  bomber  +  1  SEAD 

8,  4,  or  2  defenders  +  1  SAM  with  3  missiles 

In  these  scenarios  SEAD  aircraft  engage  the  SAM  site  (the  SAM  has  unlimited 
missiles  for  this  duel).  The  SAM  site  does  not  engage  the  escort  fighters,  and  the 
defending  fighters  do  not  engage  the  SEAD  aircraft.  If  the  SEAD  aircraft  are  de¬ 
stroyed  by  the  SAM,  the  SAM  site  (with  up  to  three  missiles)  engages  the 
bombers.  In  scenarios  with  no  SEAD  aircraft,  the  SAM  site  always  engages  the 
bombers  with  three  missiles.  The  limitation  to  three  missiles  accounts  for  the  lim¬ 
ited  time  interval  on  which  the  bombers  are  on  their  bomb  run,  rather  than  for  any 
limit  to  the  number  of  missiles  at  the  SAM  site. 

New  Options 

Several  options  have  been  added  to  SLAACM  to  more  accurately  represent  en¬ 
gagements  and  campaigns  and  to  allow  the  analyst  to  change  parameters  that  were 
previously  hardcoded.  These  new  options  include  the  following: 

♦  Escaping  pure  bombers  complete  mission.  SLAACM  has  been  revised  so 
that,  in  the  case  that  all  escorts  are  killed,  the  pure  bombers  either  escape 
and  go  home  or  escape  and  complete  their  mission 

♦  Escaping  fighter-bombers  complete  mission.  Fighter-bombers  fight 
through  as  before,  and  optionally  go  home  or  complete  their  missions. 
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♦  User-definable  breakpoints.  The  analyst  can  now  select  loss  breakpoints 
for  both  offensive  and  defensive  forces.  These  breakpoints  are  used  in  the 
engagement  calculations  to  detennine  the  bivariate  loss  probabilities. 

♦  Attacker  and  Defender  loss  values.  The  analyst  now  has  the  option  to  as¬ 
sign  the  payoff  function  values  the  attacker  and  defender  place  on  their 
own  loss  and  on  killing  an  opponent.  As  before,  the  LAD  (local  area  de¬ 
fender)  buttons  associated  with  each  defender  aircraft  type,  can  be  used  to 
force  that  defender  type  to  ignore  his  own  losses. 

♦  Bomb  payload  value.  The  bomb  payload  values,  calculated  from  the 
bomber  tonnage  and  CEP  payloads,  are  hardcoded,  but  are  now  displayed 
on  the  “RedBombs”  worksheet,  along  with  the  Attacker  and  Defender 
loss  values  to  give  the  analyst  insight  into  the  Attacker  and  Defender 
payoff  values. 

SAM  and  SEAD  Features 

The  following  subsections  describe  the  engagement  models  used  in  the  models  of 
SLA  ACM’ s  new  features. 

SAM-SEAD  Engagements 

When  attack  packages  include  SEAD  aircraft,  2  SEADs  engage  the  SAM  site  de¬ 
fending  the  package’s  target.  This  engagement  takes  place  before  any  of  the 
package’s  bombers  reach  the  target.  If  the  SEADs  destroy  the  site,  then  the  site 
does  not  engage  the  bombers. 

Two  options  are  available  for  the  SEADs  breakaway  behavior:  the  SEADs  may 
engage  to  annihilation,  or  they  may  disengage  on  sustaining  1  loss. 

We  model  the  SAM-SEAD  engagement  as  a  stochastic  Lanchester  engagement, 
in  which  the  times  between  kills  by  the  two  adversaries  each  are  exponentially 
distributed.  We  denote  the  parameter  of  the  SAM  site’s  time-to-kill  distribution 
by  ks  and  call  that  of  the  SEAD’s  distribution  kf.  The  parameter  kf  is  equal  to  the 
mean  time  between  kills  by  the  SEADs;  ks  is  the  mean  time  between  kills  by  the 
site. 

We  denote  the  state  of  the  engagement  by  the  ordered  pair  (f,  s),  where  f  denotes 
the  number  of  fighters  present.  The  parameter  s  has  the  value  1  if  the  site  is  active, 
0  if  it  is  destroyed.  Figure  2-9  diagrams  the  engagement  when  the  SEADs  fight  to 
annihilation.  The  labels  on  the  arrows  indicating  state  transitions  are  the  rates  of 
the  transition  indicated. 
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Figure  2-9.  SAM-SEAD  Engagement  Diagram,  SEADs  Break  at  0 
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(0,1) 


Let  Pij(t)  denote  the  probability  that  the  system  is  in  state  (i,  j)  at  time  t.  The 
forward  Chapman-Kolmogorov  equations  for  the  engagement  can  be  written  by 
inspection  of  Figure  2-9: 


-  21 


20 


10 


-  01 


=  -(2kf+ks)P21 
=  2kfP21 

=  -(kf  +ks)Pn  +ksP21 


—  kfPn 
=  ksP„ 


[Eq.  2-19] 


The  system  begins  in  state  (2,  1),  so  that  at  t  =  0,  all  P;j  are  zero  except  P21.  This 
initial  condition  and  equations  2-19  determine  the  complete  evolution  of  P;j. 


Consistent  with  SLAACM’s  treatment  of  other  stochastic  engagements,  we  use 
the  long-time  limiting  values  P”  of  the  Pj,  for  the  outcome  probabilities  of  the 

SAM-SEAD  engagement.  The  method  of  Appendix  A  allows  us  to  write  these 
values  explicitly.  In  the  long-time  limit  P2i  and  Pn  are  zero,  and 


00  _  2p 

20  ~  ,  , 

2p  + 1 


P 

(2p  +  l)(p  +  l) 


(2p  +  l)(p  +  l) 


[Eq.  2-20] 
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We  use  the  symbol  p  for  the  ratio  kf/ks.  The  probability  that  the  site  is  up  isP0" , 
the  probability  that  the  site  is  down  with  0  SEADs  lost  is  Pf0 ,  and  the  probability 
that  the  site  is  down  with  1  SEAD  lost  is  P^  . 

Similar  considerations  show  that,  for  a  SAM-SEAD  engagement  in  which  the 
SEADs  break  on  sustaining  1  loss,  the  probability  that  the  site  is  down  i s  P2y0 , 
which  is  also  the  probability  that  0  SEADs  are  lost.  For  this  engagement,  the 
probability  that  the  site  is  up  and  1  SEAD  is  lost  is  1  -  P2^  =  l/(2p  +1). 

SLAACM  incorporates  the  results  of  this  section  to  treat  the  outcomes  of  SAM- 
SEAD  engagements. 

SAM-Bomber  Engagements 

We  are  interested  in  the  probability  Sm(j,k)  that  j  of  k  attacking  bombers  survive 
when  a  SAM  site  fires  m  shots  at  the  bombers.  We  assume  that  the  site’s  battle 
management  is  sufficiently  sophisticated  that  the  outcomes  of  all  the  shots  fired 
are  statistically  independent,  with  single-shot  kill  probability  p.  Then  when  m  <  k, 
obviously 


Sm(j,k)  =  B(k- j,m,p),  0  <  j  <  k  [Eq.  2-21] 

where  B(i,  j,  p),  the  binomial  probability  distribution  function,  represents  the  prob¬ 
ability  of  i  successes  in  j  trials,  when  the  single-trial  probability  of  success  is  p. 

When  m  >  k,  Sm  (j,k)  may  be  reckoned  with  the  following  considerations.  Let  us 
begin  with  a  simple  example,  when  m  =  4  and  k  =  2.  Introducing  a  system  state 
variable  (i,  j),  where  i  is  the  number  of  shots  remaining  and  j  the  number  of  aircraft 
remaining,  the  system  evolves  as  shown  in  Figure  2-10.  In  that  figure,  q  =  1  -  p. 

Figure  2-10.  Four  SAM  Shots  vs.  Two  Bombers 
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The  diagram  of  Figure  2-10  proceeds  in  the  same  way  as  the  diagram  for  the  bi¬ 
nomial  distribution,  until  two  shots  have  been  fired.  For  three  shots,  since  the  state 
(2,  0)  is  an  absorbing  boundary  state  (the  engagement  ends  when  both  targets  are 
destroyed),  the  left-going  arrow  from  that  state  to  the  state  (1,  0)  is  absent.  Thus 
the  probability  of  state  (1,0)  is  not  B(2,  3,  p),  but,  rather,  B(2,  3,  p)  -  qB(2,  2,  p). 
Similarly,  the  probability  of  state  (0,  0)  is  B(2,  4,  p)  -  qB(2,  3,  p). 

The  probabilities  of  states  (0,  j),  1  <  j  <  2  may  be  read  from  the  bottom  row  of  the 
diagram  of  Figure  2-10.  Since  each  of  these  states  has  all  the  predecessors  of  the 
binomial  distribution,  they  are  binomial,  specifically,  B(2  -  j,  4,  p). 

Thus  for  the  example  of  Figure  2-10, 

S4  (j,  2)  =  B(2  -  j,4,p),  1  <  j  <  2  [Eq.  2-22] 


and 


S4  (0,2)  =  X  [B(2,4  -  n,  p)  -  qB(2,4  -  n  - 1, p)]  +  B(2,2,  p) .  [Eq.  2-23] 

n=0 


We  may  evaluate  Sm(j,k)  in  general  with  an  extension  of  the  diagram  of  Figure  2-10. 
The  extended  diagram  begins  in  state  (m,  k),  as  shown  in  Figure  2-11. 

Figure  2-11.  Extended  Diagram 


(m-k-1  ,  0) 


(0  ,  k) 


(0  ,  k-1) 


(0  ,1)  (0 , 0) 
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The  diagram  of  Figure  2-11  proceeds  in  the  same  way  as  the  diagram  for  the  bi¬ 
nomial  distribution,  until  k  shots  have  been  fired.  For  k  +  1  shots,  since  the  state 
(m  -  k,  0)  is  an  absorbing  boundary  state  (the  engagement  ends  when  all  k  targets 
are  destroyed),  the  left-going  arrow  from  that  state  to  the  state  (m  -  k  -  1,  0)  is 
absent.  Thus  the  probability  of  state  (m  -  k  -  1,  0)  is  not  B(k,  k  +  1,  p),  but, 
rather,  B(k,  k  +  1,  p)  -  qB(k,  k,  p).  As  more  shots  are  fired,  the  probability  of  state 
(m  -  k  -  n,  0)  is  similarly  seen  to  be  B(k,  k  +  n,  p)  -  qB(k,  k  +  n  -  1,  p),  until  all 
m  shots  have  been  fired,  bringing  us  to  the  bottom  row  of  the  diagram.  The  prob¬ 
ability  of  state  (0,  0)  is  B(k,  m,  p)  -  q  B(k,  m  -  1,  p). 

The  probability  of  ending  states  (0,  j),  read  from  the  bottom  row  of  the  diagram  of 
Figure  2-11,  is  B(k  -  j,  m,  p),  1  <  j  <  k.  Consequently 

Sm  (j,  k)  =  B(k  -  j,  m,  p),  1  <  j  <  k  [Eq.  2-24] 


and 


m-k-l 

Sm(0,k)  =  ^[B(k,m-n,p)-qB(k,m-n 


n=0 


l,p)]  +  B(k,k,p) . 


[Eq.  2-25] 


Now  let  us  consider  the  expected  value  of  the  number  of  bombers  destroyed  by 
SAMs,  for  the  cases  involving  SAMs.  We  consider  attack  package  4  as  a  repre¬ 
sentative  case.  In  this  package,  2  SEAD  aircraft  precede  4  escorts  accompanying 
2  bombers.  Each  package’s  target  is  defended  by  a  SAM  site. 

The  SEADs  engage  the  site.  The  defending  fighters  first  engage  the  escorts;  sur¬ 
vivors  engage  the  bombers.  If  SLAACM  II’s  user  has  chosen  options  pennitting 
victorious  fighter-bombers,  or  escaping  bombers,  to  continue  to  the  package’s  tar¬ 
get,  and  if  the  site  survives  the  SEADs’  attack,  bombers  reaching  the  target  are 
attacked  by  the  SAM  site. 

When  bombers  do  reach  the  target,  either  one  or  two  arrive.  A  viable  site  then 
fires  a  number  m  of  shots  at  the  bombers;  m  is  a  user  option  in  SLAACM  II. 

Thus  the  expected  number  Sk  of  bombers  destroyed  by  the  site  is 

Sk  =  PB(1)  x  Sm  (0,1)  +  PB(2  )  x  (2Sm  (0,2)  +  Sm  (1,2))  [Eq.  2-26] 

where  PB(j)  denotes  the  probability  that  j  bombers  reach  the  site.  Calculation  of 
Sk  thus  requires  PB(j),  and  the  ability  to  calculate  Sm(j,k)  .  Calculating  Sm(j,k)  is 
readily  done  using  Equations  2-24  and  2-25.  To  evaluate  PB(j),  one  may  consider 
the  activity  network  diagram  of  Figure  2-12. 

Please  note  that  the  sequence  of  events  diagrammed  in  Figure  2-12  is  appropriate 
for  calculating  the  probabilities  of  the  number  of  bombers  reaching  the  site,  but  it 
is  not  necessarily  appropriate  for  other  probability  calculations. 
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Figure  2-12.  Activity  Network  Diagram 


SEADs  defeat  site 


Referring  to  Figure  2-12,  we  see  that,  when  escaping  bombers  or  victorious 
fighter-bombers  do  not  continue  to  the  target,  the  only  way  in  which  bombers 
reach  the  target  is  when  the  SEADs  have  not  defeated  the  SAM  site,  and  the  es¬ 
corts  have  defeated  the  defenders.  Thus,  in  this  case, 


4 

PB(2)  =  P(site  up)  ^  Ares(d  min- 1,  j) 

j=e  min 


PB(1)  =  0 


[Eq.  2-27] 


where  Ares(i,  j)  denotes  the  probability  that  i  defenders  and  j  escorts  survive  the 
defender-escort  engagement,  dmin  denotes  the  minimum  number  of  defenders 
willing  to  continue  the  engagement,  and  emin  denotes  the  minimum  number  of 
escorts  willing  to  continue.  Let  PS(j)  be  the  probability  that  j  SEADs  survive  the 
SEAD-site  engagement;  then  when  the  SEADs  fight  to  annihilation,  P(site  up)  = 
PS(0).  If  the  SEADs  break  at  1  loss,  then  P(site  up)  =  PS(1). 
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If  escaping  bombers  and/or  victorious  fighter-bombers  do  continue  to  their  tar¬ 
gets,  then  when  bombers  are  fighter-bombers,  PB(j)  is  given  by 


PB(2)  =  P(site  up) 
PB(1)  =  P(siteup) 


ndef 


^  Ares(d  min- 1,  j)  +  ^  Ares(k,  e  min-  l)Bres(kb,2) 

j=emin  k=d  min 

ndef 

^Ares(k,  emin-  l)Bres(kb,l) 


k=d  min 


[Eq.  2-28] 


In  Equation  2-28,  Bres(i,  j)  denotes  the  probability  that  i  defenders  and  j  bombers 
survive  the  defender-bomber  engagement,  and  kb  =  dmin  -  1  when  the  bombers 
are  fighter-bombers,  while  kb  =  k  when  the  bombers  are  escaping  bombers. 

Calculations  of  the  expected  number  of  bombers  killed  by  sites  proceed  similarly 
for  attack  packages  5,  6,  and  7,  which  are  modifications  of  attack  package  4. 

Escaping  Bomber  Characteristics 

A  class  of  engagement  models  in  which  bombers  never  kill  defenders,  but  may 
escape  from  defenders,  is  discussed  in  “T  AC  AIR  Analytical  Support,  Interim 
Mathematical  Report  1,”  published  in  January  2008.  For  completeness,  we  in¬ 
clude  the  basic  discussion  of  the  “escaping  bomber”  engagement  models  here. 

Typically,  bomber  aircraft,  other  than  fighter-bomber  aircraft  used  as  bombers, 
are  quite  weak  compared  to  defending  fighter  aircraft.  Specifically,  the  mean  time 
between  kills  of  fighters  by  bombers  is  much  longer  than  the  mean  time  between 
kills  of  bombers  by  fighters.  This  is  a  realistic  model  of  what  would  happen  if  the 
fighters  actually  engaged  the  bombers. 

But  the  bombers  might  be  able  to  escape  the  fighters,  instead  of  engaging  with 
them.  Consequently,  a  better  model  of  engagements  between  bombers  and  de¬ 
fenders  should  not  allow  the  bombers  to  kill  a  defender,  but  should  allow  for  the 
possibility  that  bombers  escape  from  the  engagement.  The  following  paragraphs 
develop  a  model  of  that  situation  and  explain  how  the  new  model  has  been  inte¬ 
grated  into  an  extended  version  of  SLAACM. 

We  consider  an  engagement  of  i  defending  fighters  with  a  flight  of  4  bombers,  in 
which  the  bombers  do  not  attempt  to  destroy  the  fighters,  but,  rather,  attempt  to 
evade  them.  Consistent  with  SLAACM’s  general  modeling  approach,  we  assume 
that  the  times  between  kills  made  by  the  fighter  aircraft  are  exponentially  distrib¬ 
uted,  identically  and  independently,  with  mean  time  1/kb.  We  further  assume  that 
the  times  between  escapes  made  by  the  bomber  aircraft  are  identically  and  inde¬ 
pendently  exponentially  distributed,  with  mean  time  l/kx. 

Let  the  state  of  the  engagement  be  described  by  the  triple  (i,  j,  k),  where  i  denotes 
the  number  of  fighters  present,  j  the  number  of  bombers  present,  and  k  denotes  the 
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number  of  bombers  that  have  escaped  from  the  engagement.  Then  the  engage¬ 
ment’s  state  transition  diagram  is  as  shown  in  Figure  2-13. 

Figure  2-13.  State  Transition  Diagram 


(i,  o,  0) 


(i.  0,  1) 
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The  state  transition  diagram  of  Figure  2-13  implies  the  following  set  of  15  evolu¬ 
tion  equations  (forward  Chapman-Kolmogoroff  equations)  for  the  probabilities 
Pi,j,k(t)  that  the  system  is  in  state  (i,  j,  k)  at  time  t: 

Ki,k  =  “A  +  Jk,  +  ikbpij+i,k  +  U  +  l)k,pi.j+ik- 1.  4  -  7  -  °’  o  <  k  <  4  -  j  [Eq.  2-29] 
where  Pl k  (t)  =  0  if  j  >  4  or  k  <  0. 

Initially,  the  system  is  in  state  (i,  4,  0),  so  that 

A.0(°)  =  1.  A,a-(°)  =  0,j±4,k±0.  [Eq.  2-30] 
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It  is  convenient  to  introduce  the  nondimensional  time  z  =  kxt .  For  the  independent 
variable  x,  Equations  2-29  and  2-30  are 

Kj*  =  -<iK  +  +  ilcPij+hk  +  U  +  !)^,y+u-i  >  4  >  7  >  0,  0<k<4-j  [Eq.  2-31] 

and 


n4,0(°)  =  l’  Pij,k  (0)  =  0, 7  ^  4,  A:  ^  0 . 


[Eq.  2-32] 


In  Equation  2-3 1  the  prime  denotes  differentiation  with  respect  to  x, 
and  k  =  kb/  kx .  The  parameters  kb  and  kx  appear  in  Equation  2-3 1  only  in  the  sin¬ 
gle  parameter  k  .  (Strictly  speaking,  the  probabilities  as  functions  of  x  are  distinct 
from  the  probabilities  as  functions  of  t,  but  for  clarity  and  simplicity,  we  will  not 
introduce  distinct  notation  for  them.) 

The  initial  value  problem  of  Equations  2-31  and  2-32  may  be  solved  analytically 
or  numerically  by  standard  methods  to  give  the  time  evolution  of  the  state  prob¬ 
abilities.  Alternatively,  one  may  consider  the  long-time  limit  of  the  state 
probabilities.  This  limit  is  relevant  if  engagements  are  likely  to  last  for  times 
much  longer  than  the  mean  time  for  a  bomber  to  escape. 

The  long-time  limiting  probabilities  may  be  calculated  in  finitely  many  steps.  In¬ 
tegrating  Equations  2-3 1  with  respect  to  x  from  0  to  infinity  and  using  the  initial 
conditions  of  Equations  2-32  leads  to  a  system  of  10  linear  algebraic  equations  for 
the  integrals  from  0  to  infinity  of  the  transient-state  probabilities,  and  to  a  system 
of  5  expressions  for  the  long-time  limiting  values  of  the  probabilities  of  the  ab¬ 
sorbing  boundary  states,  that  is,  of  states  (i,  0,  k),  0  <  k  <  4,  as  linear 
combinations  of  integrals  of  the  transient  state  probabilities. 


The  10  equations  for  the  integrals  (j>i ,  k  of  the  transient-state  probabilities  are 


0/,4,O 
0i,j,i k 


1 

i/c +  4 

— U.’  is./<3,  o<i<4-y 

IK  +  7  IK  +  7 


[Eq.  2-33] 
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Equations  2-33  may  be  solved  recursively,  speeding  up  the  solution  process  con¬ 
siderably.  With  the  solutions,  the  long-time  limiting  probabilities  P™.  k  for  the 
absorbing  boundary  states  are  given  by 

^.0.0  =  /V4,,o 

PiM  =  +  4  1,0 

Kai  =  iKftxi  +  4u  [Eq.  2-34] 

Pi, 0,3  =  1,3  +  4  1,2 

PUO. 4  -  4l,3 

Equations  2-33  and  2-34  give  means  of  extending  SLAACM  to  model  bomber- 
defender  engagements  in  which  the  bombers  may  escape.  Like  the  equations  for 
the  M  vs.  N  stochastic  Lanchester  equations  of  SLAACM’ s  present  engagements, 
“escaping-bombers”  engagements  that  go  to  the  long-time  limit  have  only  one 
parameter,  k  .  These  models  are  available  in  current  versions  of  SLAACM. 
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Chapter  3 

How  SLAACM  Works: 

Bomber  Effectiveness  Parameter 


The  basic  campaign  modeled  in  SLAACM  is  a  bomber  attack  against  defended 
targets.  The  attacker’s  bombers  have  different  levels  of  effectiveness,  based  on 
the  attack  package.  The  offensive  effectiveness  of  a  bomber  depends  on  two  fac¬ 
tors:  the  payload  of  bombs  carried,  and  the  accuracy  with  which  those  weapons 
can  be  delivered  to  targets.  We  developed  a  measure  of  bomber  effectiveness 
based  on  these  factors. 

We  considered  first  the  relation  between  blast  overpressure  p  and  the  properties  of 
a  structure  designed  to  resist  that  overpressure.  We  used  the  simple  model  of  rein¬ 
forcing  structure  shown  in  Figure  3-1. 

Figure  3-1.  Reinforcing  Structure  Model 


— H  a  I" — 


In  this  simple  model,  two  columns  reinforce  the  structure.  Their  thicknesses  are  A, 
their  heights  are  h,  and  the  spacing  between  their  centers  is  d.  Their  dimension 
into  the  page  is  w. 


The  reinforcing  columns  act  as  cantilevers  against  lateral  loads.  A  pressure  p  act- 

pwh2 

ing  on  the  left  face  of  this  simple  structure  will  induce  axial  loads -  in  each 

2d 

column,  tension  in  the  left  and  compression  in  the  right.  The  pressure  will  also 
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induce  shear  loads 


pwh 

2 


in  each  column.  Since  both  loads  are  proportional  to  p, 


the  maximum  stress  induced  in  the  columns  will  also  be  proportional  to  p.  If  the 
columns’  dimensions  are  adjusted  so  that  a  material  of  given  maximum  allowable 
stress  can  support  the  loads  induced  by  the  overpressure,  the  area  of  the  columns’ 
bases  wA  will  be  proportional  to  the  maximum  stress  and,  therefore,  also  to  p. 


Thus  the  volume  of  the  reinforcing  structure,  2wAh,  will  be  proportional  to  p.  As¬ 
suming  that  the  cost  of  the  reinforcing  structure  is  proportional  to  its  volume  leads 
to  the  conclusion  that  the  cost  of  reinforcement  will  be  proportional  to  overpres¬ 
sure  p. 


Now  let  us  consider  how  blast  overpressure  varies  with  explosive  yield  E.  Dy¬ 
namic  overpressure  p  is  proportional  to  the  square  of  the  air  velocity  v  immedi¬ 
ately  behind  the  blast  wave  produced  by  the  explosion.  Dimensional  analysis 
shows  how  v  varies  with  E. 1  In  addition  to  E,  the  parameters  affecting  v  are  the 
undisturbed  air  density  p  and  the  time  t  required  for  the  blast  wave  to  reach  the 
locations  of  interest. 


According  to  the  principles  of  dimensional  analysis,  v  can  be  expressed  as  the 
product  of  a  term  E“p^t7  and  a  dimensionless  function  of  all  possible  dimensionless 
combinations  of  E,  p,  and  t.  If  E“ppty  is  equal  to  a  velocity,  that  term  must  have  the 
dimension  L/T,  where  L  represents  dimension  length  and  T  dimension  time.  Now, 
energy  E  has  dimension  ML2/T2,  where  M  denotes  dimension  mass;  p  has  dimen¬ 
sion  M/L3,  and,  of  course,  t  has  dimension  T.  Thus 


( 


V 


ML 


2  A 


[Eq.  3-1] 


Equation  3-1  implies  that 

a  +  P  =  0 

2a-3p  =  l  [Eq.  3-2] 

-2a  +  y  =  -1 


The  unique  solution  to  Equation  3-2  is  a  =  1/5,  P  =  -1/5,  y  =  -3/5.  Similar  ar¬ 
guments  to  those  leading  to  Equations  3-1  and  3-2  show  that  there  are  no  non- 
dimensional  combinations  of  E,  p,  and  t,  other  than  a  constant.  Therefore 


i 


V  oc 


vP; 


[Eq.  3-3] 


1  L.  Sedov,  Similarity’  and  Dimensional  Methods  in  Mechanics  (New  York:  Academic  Press, 
1959). 
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and  the  dynamic  overpressure  p  will  satisfy 

2  3  6 

p  oc  pv2  cc  E5p5t  5 .  [Eq.  3-4] 


Similar  dimensional  considerations  lead  to  the  conclusion  that  the  radius  r  of  the 
blast  wave  satisfies  the  proportionality 


r  oc 


E 

vp; 


i 

m  2 

t5, 


[Eq.  3-5] 


From  proportionalities  (Equations  3-4  and  3-5),  it  follows  that  the  radius  rL  at 
which  the  overpressure  falls  to  Pl  satisfies 


rL  00 


^E" 

Pv 


1 

3  _1 

Pl~3. 


[Eq.  3-6] 


Thus  the  lethal  radius  of  a  blast  wave  increases  as  the  cube  root  of  the  blast’s  en¬ 
ergy  and  decreases  as  the  cube  root  of  the  lethal  overpressure. 

Since  in  our  model  the  value  of  a  structure  is  proportional  to  the  overpressure  re¬ 
quired  to  destroy  it,  proportionality  (Equation  3-6)  implies  that  the  lethal  radius  of 
a  blast  with  energy  E  varies  directly  with  the  cube  root  of  E,  and  inversely  with 
the  cube  root  of  the  value  of  the  structures  it  affects.2 


We  now  use  the  information  about  the  variation  of  lethal  radius  with  explosion 
energy  and  structure  value  to  see  how  bombs’  single  shot  kill  probabilities 
(sspk’s)  vary  with  explosion  energy  and  value  of  structures  affected. 

The  sspk  of  a  bomb  with  lethal  radius  tl,  delivered  with  circular  error  probable 
(CEP),  is 


sspk  =  1 


r£ 

^ 1  )  CEP2 

UJ 


[Eq.  3-7] 


Introducing  the  variation  of  tl  with  explosive  energy  E  and  value  V  obtained 
above,  we  have 


sspk  =  1 


kE2/3V~2/3 

/  1  )  CEP2 

UJ 


[Eq.  3-8] 


2  These  conclusions,  obtained  with  simple  dimensional  analyses,  agree  with  those  from  the 
exact  Taylor-Sedov  solution  [Sedov,  1.  c.  ante]  for  values  immediately  behind  the  shock. 
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where  k  is  a  constant. 


Now, 


kE2/3y-2/3 


=  kE2/3V~2/3 

R  ref  vref 


V^ref  J 


V, 


ref 


V 


f  CEPref  V 
l  CEP  J 


Defining  the  parameter  ki  by 


k  =kE2/3V-2/3 

K1  -K£Lref  Vref  • 


[Eq.  3-9] 


[Eq.  3-10] 


we  have 


sspk  =  1 


f— 7i 

'Vref> 

UefJ  > 

{2J 


CEPref]2 
CEP  J 


[Eq.  3-11] 


Let  Eref  be  the  energy  of  a  500  lb  bomb.  Let  Vref  be  the  value  of  a  target  hardened 
to  withstand  a  500  lb  bomb.  We  will  consider  500  lb,  1,000  lb,  and  2,000  lb 
bombs  and  assume  that  bombs  are  always  matched  to  targets,  so  that  a  bomber  on 
a  mission  against  500  lb  targets  carries  500  lb  bombs,  and  one  on  a  mission 
against  1,000  lb  targets  carries  1,000  lb  bombs,  and  that  one  going  against  2,000 
lb  targets  carries  2,000  lb  bombs. 

With  those  assumptions,  the  product  of  the  term  involving  E  and  the  tenn  involv¬ 
ing  V  in  Equation  3-11  is  always  1,  and 


sspk  =  1 


( 1  \ 

Jv 


CEPref 

CEP 


[Eq.  3-12] 


We  arbitrarily  set  the  constant  ki  to  1;  we  only  want  a  systematic  way  to  compute 
the  value  to  Red  of  a  bomber  with  given  payload  P  that  is  capable  of  delivering 
bombs  with  relative  accuracy  CEP/CEPref. 

With  this  assumption,  the  value  of  the  expected  number  of  targets  killed  by  such  a 
bomber  carrying  500  lb  bombs  is 


V,= 


P 

500 


sspk0V5oo 


[Eq.  3-13] 
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where 


sspk0  =  1 


f 


n 


CEPreff 
CEP  J 


U ) 


[Eq.  3-14] 


and  where  V500  is  the  value  associated  with  a  target  hardened  to  withstand  500  lb 
bombs. 

When  the  bomber  is  loaded  with  1,000  lb  bombs,  the  value  of  the  expected  num¬ 
ber  of  targets  killed  is 


sspk0(2)225Vs0l),  [Eq.  3-15] 

since  the  value  of  a  target  goes  as  the  2/5  power  of  the  explosive  yield  required  to 
destroy  it. 

Finally,  if  the  bomber  is  loaded  with  2,000  lb  bombs,  the  value  of  the  expected 
number  of  targets  killed  is 

sspk0(4)2/5V!0C1.  [Eq.  3-16] 

Using  this  infonnation,  SLAACM  assigns  the  average  of  Vi,  V2,  and  V3,  divided 
by  V500,  as  the  value  of  the  bomber. 

The  value  for  the  CEP/CEPref  ratio  is  set  by  the  analyst  for  each  bomber  type  on 
the  “ExRatios”  worksheet.  We  use  the  ratio  of  1  if  the  bomber  cannot  deliver 
guided  (smart)  bombs  and  a  lower  value,  such  as  0.1  to  0.4,  if  it  can.  Bomber  de¬ 
structive  power  is  used  in  the  optimizer  payoff  functions,  but  it  has  been  useful  to 
report  the  bombs  delivered  in  terms  of  tonnage  in  the  output  tables  and  graphs. 

The  user  can  display  the  tonnage  of  smart  bombs  delivered  by  setting  set  a 
CEP/CEPref  ratio  display  value  on  the  “RedBombs”  worksheet. 

As  mentioned  in  Chapter  2,  the  destructive  power  of  the  bomber  payloads  are  now 
displayed  on  the  “RedBombs”  worksheet  to  give  the  analyst  insight  into  the  pay¬ 
off  values  used  in  the  optimizer. 


v2  = 


2000 
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Chapter  4 

How  SLAACM  Works: 

Campaign  Logic  and  Optimization 


This  chapter  discusses  the  attacker  and  defender  payoff  logic  and  optimization 
algorithms  used  in  SLAACM. 

Specifying  Attacker  and  Defender  Payoff 
Functions 


SLAACM  determines  optimal  attacker  attacks  and  optimal  defender  responses 
using  payoff  functions  for  each  side.  This  section  describes  in  detail  the  calcula¬ 
tions  used  for  the  basic  air  defense  scenario.  Other  scenarios  use  similar  logic. 

Attacker’s  Payoff  Function 

Attacker  values  an  attacker  package  by  considering  the  benefit  of  getting  bombers 
through  and  of  downing  defender  aircraft.  The  attacker  may  also  consider  penal¬ 
ties  for  losing  its  own  aircraft,  assigning  different  values  to  losing  a  fighter  and  to 
losing  a  bomber.  Specifically,  for  the  basic  air  defense  scenario,  a  package’s  value 
to  the  attacker  is  determined  in  this  way:  An  attacker  package  is  defined  by  the 
triple  (a,  c,  b),  where  a  is  the  advanced  escort  type,  c  the  close  escort  type,  and  b 
the  bomber  type  (the  bombers  may  be  fighter-bomber  aircraft).  In  Attack 
SLAACM,  attack  packages  can  differ  in  composition  from  the  basic  air  defense 
(a,  c,  b)  triple.  Nevertheless,  the  attacker’s  expected  payoff  is  calculated  using  the 
same  method.  The  outcome  statistics  for  an  engagement  between  the  package  and 
the  “planning  defender”  aircraft  (the  user  identifies  a  defender  type  for  the  at¬ 
tacker’s  planning)  are  determined,  and,  from  these,  the  expected  number  Bs  of 
surviving  bombers,  the  expected  number  RFk  of  lost  fighters,  and  the  expected 
number  Dk  of  downed  defenders  are  computed.  Then  VR,  the  attacker’s  value  for 
the  package,  is  given  by 

VR  =  clBs  +c2Dk  -c3RFk  -c4(4-Bs).  [Eq.4-1] 

The  positive  parameter  Ci  is  currently  set  at  the  value  of  a  bomber  computed  by  the 
method  described  in  Chapter  3.  Values  of  C2  and  C3  are  inputs  from  the  user.  When 
the  attacker  believes  the  defenders  to  be  loss-averse,  the  value  of  C2  is  typically  set 
at  20,  reflecting  the  fact  that  the  attacker  believes  the  defender  is  likely  to  abandon 
the  campaign  if  their  side  incurs  many  losses.  For  an  attacker  not  overly  concerned 
about  losses,  the  value  of  C3  is  typically  set  at  zero;  also,  if  the  attacker  faces  de¬ 
fenders  that  are  strong  in  comparison  with  his  aircraft  (values  of  k  significantly 
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greater  than  one),  even  values  of  C3  that  are  small  compared  with  one  may  cause  the 
attacker  not  to  dispatch  any  packages  at  all.  Currently,  C4  is  set  equal  to  C3. 

Defender’s  Payoff  Function 

The  defender  values  engaging  a  given  attack  package  by  considering  the  benefit 
of  downing  the  attacker’s  fighters  and  bombers,  and  the  penalty  of  losing  its  own 
aircraft.  Specifically,  the  value  to  the  defender  of  engaging  an  attacker  package 
with  fighters  of  a  given  type  is  calculated  by  computing  the  outcome  statistics  for 
that  composite  engagement  (described  in  Chapter  3)  and,  from  them,  the  expected 
value  Bk  of  bombers  killed,  the  expected  value  RFk  of  attack  fighters  killed,  and 
the  expected  value  Dk  of  defenders  killed.  Then  the  value  Vb  of  the  engagement  to 
the  defender  is  evaluated  as 

VB=d1Bk-d2Dk+d3RFk.  [Eq.  4-2] 

Currently,  di  is  set  to  the  attacker’s  value  for  the  bomber  type  in  the  package,  and 
d2  and  d3  are  user  inputs.  (A  prime  reason  for  allowing  user  assignment  of  the 
payoff  values  is  that  they  are  likely  to  be  different  for  friendly  defense  and 
friendly  attack  campaigns.) 

Determining  Attacks  and  Defenses 


A  key  part  of  SLAACM’s  operation  is  determining  the  engagements  that  occur 
each  day.  This  is  done  by  finding  optimal  sets  of  attacker  packages  dispatched  by 
the  attacker  and  by  determining  the  defender’s  response  in  two  ways:  the  response 
of  smart  defenders  that  know  the  makeup  of  the  attacker  packages  before  inter¬ 
cepting  them  and  make  optimal  defenses,  and  the  response  of  not-smart  defenders 
that  encounter  attack  packages  randomly.  This  section  explains  these  calculations. 

Determining  the  Attacker  Packages 

As  the  aggressor,  the  attacker  always  plans  its  attacks.  This  leads  to  an  integer 
programming  (IP)  problem,  which  SLAACM  treats  in  one  of  two  ways:  a  classic 
“greedy”  heuristic,  or  an  exact  solution  using  the  commercial  general  optimiza¬ 
tion  package  LINGO™.  This  section  describes  those  calculations. 

Attacker’s  Integer  Programming  Problem 

Each  day,  the  attacker  has  a  set  of  fighter,  bomber,  and  fighter-bomber  aircraft  to 
deploy.  In  some  scenarios,  the  attacker  also  has  SEAD  aircraft.  In  SLAACM,  the 
attacker  chooses  the  set  of  attacker  packages  to  maximize  the  total  of  the  expected 
payoff  of  the  packages  dispatched  (evaluating  payoff  as  described  above  in  “At¬ 
tacker’s  Payoff  Function”),  constrained  by  the  available  forces. 
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Here  we  formalize  the  IP  problems  performed  in  the  LINGO™  model  within 
SLAACM.  The  main  idea  behind  using  an  IP  tool  is  to  solve  the  following  opti¬ 
mization  problem:  the  attacker  wishes  to  come  with  optimal  attacker  packages 
based  on  a  detennined  payoff  function  (this  is  the  first  optimization).  Then,  the 
defender  aircraft  that  are  sophisticated  enough  to  have  a  priori  knowledge  of  what 
attacker  packages  are  coming  may  optimize  their  defending  packages  accordingly 
(the  second  optimization). 

A  basic  defense  scenario  attacker  package  consists  of  4  advanced  escorts,  4  close- 
in  escorts,  and  4  bombers.  Scenarios  for  Attack  SLAACM  also  include  SEADs 
and  varying  numbers  of  escorts  and  bombers. 

Since  some  fighters  may  be  bombers,  for  the  basic  air  defense  SLAACM  sce¬ 
nario,  we  could  have  a  set  of  T  possible  attacker  packages,  where,  given  R  types 
of  fighters  and  bombers,  we  could  have  |T|  =  R3 .  In  reality,  however,  most 
mathematically  feasible  combinations  of  packages  are  unrealistic  for  warfare  (one 
would  not,  for  example,  dispatch  a  package  of  heavy  bombers  escorting  fighters). 
Since  the  number  of  reasonable  attacker  packages  is  much  smaller  than  the  com¬ 
plete  enumeration,  we  assume  those  undesirable  packages  are  removed  in  advance 
of  the  optimization;  then,  we  limit  ourselves  to  those  remaining,  reasonable  possi¬ 
bilities  in  the  integer  program.  The  degree  to  which  one  wishes  to  limit  the 
number  of  potential  packages  depends  primarily  on  the  size  of  the  integer  pro¬ 
gram,  which  is  a  function  of  the  size  of  R  . 

In  Attack  SLAACM,  an  attack  package  may  consist  of  SEAD,  escort,  and  bomber 
aircraft.  If  the  package  has  S  types  of  SEAD  aircraft,  E  types  of  escort  aircraft, 
and  B  types  of  bomber  and  fighter-bomber  aircraft,  then  |t|  =  SEB  .  Of  course, 
here  too  not  all  combinations  are  feasible  militarily. 

We  define  a  set  of  attacker  packages  as  IcT.We  define  variable  r,,  as  the  num¬ 
ber  of  aircraft  of  type  j  used  in  attacker  package  i .  To  calculate  its  payoff,  the 

attacker  plans  its  attack  assuming  each  package  will  be  confronted  by  one  specific 
defender  aircraft  type,  usually  the  most  numerous  of  the  defender  aircraft.  Based 
on  the  payoff  functions  described  previously,  each  attacker  package  has  a  certain 
payoff,  denoted  as  p; .  We  denote  the  number  of  attack  aircraft  of  type  j  as  rij . 

At  this  point,  the  attacker  can  solve  for  its  optimal  attack  strategy,  denoted  by  the 
following  integer  optimization  problem: 

i 

max^piXi 

1=1 

I 

s-t-  ZryX>  -nj  VJeJ 

i=l 

Xi  e  {0,1,2,...}. 
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We  solve  for  the  values  of  xt ,  which  represent  the  optimal  number  of  packages  of 
type  i  sent.  This  problem  is  a  variant  of  the  classic  knapsack  problem.  The  solu¬ 
tion  to  this  integer  programming  problem  provides  the  optimal  set  of  attacker 
packages  dispatched  by  the  attacker. 

Heuristic  Solution 

The  exact  integer  programming  problem  can  be  solved  using  an  IP  software  tool 
such  as  LINGO.  In  SLAACM,  the  user  can  specify  that  the  exact  solution  tech¬ 
nique  be  used,  provided  the  user’s  PC  has  a  LINGO  license.  Although  utilizing  an 
optimization  package  is  the  ideal  circumstance — it  can  guarantee  optimal  solu¬ 
tions  if  run  to  completion — many  analysts  will  not  have  LINGO  available  on  their 
PCs.  For  this  reason,  SLAACM  includes  a  greedy  heuristic,  using  the  attacker  and 
defender  payoff  functions,  that  attempts  to  approximate  the  optimum  solution. 

The  heuristic  has  been  shown  to  yield  the  optimum  solution  when  the  attacker  and 
defender  aircraft  types  and  numbers  are  limited.  However,  as  the  choices  grow, 
and  particularly  when  the  attacker  has  options  of  using  aircraft  as  fighters  or 
fighter-bombers,  the  packages  selected  by  the  heuristic  will  deviate  significantly 
from  an  IP  solution.  This  difference  will  be  compounded  when  defenders  are 
smart  and  also  optimize.  (When  defender  aircraft  are  dominant,  the  different 
package  choices  may  produce  only  small  differences  in  the  outcome  statistics.)1 

The  greedy  heuristic  calculates  what  might  be  considered  the  “obvious”  answer.  It 
is  considered  greedy  because  it  takes,  without  more  holistic  considerations,  pack¬ 
ages  with  the  largest  payoffs  first.  This  approach  can  be  suboptimal  because  it 
does  not  consider  (or  mathematically  eliminate  from  consideration)  all  combina¬ 
tions  of  package  types. 

As  we  noted  earlier,  the  payoff  for  attacker  package  i  is  denoted  as  p;.  The  heu¬ 
ristic  first  orders  the  packages  from  greatest  to  least.  Without  loss  of  generality, 
we  denote  the  largest  payoff  as  package  pi,  the  second  largest  as  P2,  and  so  on. 
Thus,  we  have  Pj  <p2<...<Pi  .  Again,  we  define  rij  as  the  number  of  attack  air¬ 
craft  of  type  j  .  The  heuristic  can  then  be  described  by  the  following  procedure. 

Step  0:  Set  Hi  =  m  for  all  j  . 

Initialize  x}  =  0  for  alii . 

Set  i  =  l. 


1  For  further  information  on  integer  programming  and  greedy  heuristics,  the  interested  reader 
should  consider  one  of  the  many  textbooks  on  integer  programming  and  optimization,  such  as 
L.  Wolsey  and  G.  Nemhauser,  Integer  and  Combinatorial  Optimization  (New  York:  Wiley,  1988). 
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Step  1 :  If  n  j  - 1  >  0  for  all  j  such  that  r;j  >  0,  then  { 

n  |  =  n  j  - 1  for  all  j  such  that  ry  >  0 

xi  =X,  +1 

Go  to  Step  1 

}  Else  { 

Go  to  Step  2  } 

Step  2:  i  =  i  +  l. 

If  Pi  <  0  ,  quit. 

If  i  =  I  + 1 ,  quit. 

Else,  Go  to  Step  1 . 

In  essence,  this  heuristic  keeps  selecting  the  package  with  the  greatest  payoff  until  it 
exhausts  one  of  the  package’s  aircraft  types  (i.e.,  “hits  a  constraint”).  It  then  looks  for 
the  package  with  the  next  greatest  payoff  (assuming  all  its  aircraft  types  are  still  avail¬ 
able).  It  keeps  selecting  that  package  until  one  of  its  aircraft  types  is  exhausted,  and 
so  forth.  It  continues  this  procedure  until  there  are  no  more  packages  with  positive 
payoffs  for  which  the  advanced  escort,  close-in  escort,  or  bomber  aircraft  types  have 
not  been  exhausted. 

It  is  important  to  reemphasize  that  this  approach  will  not  necessarily  yield  an  op¬ 
timal  solution.  Indeed,  for  sample  campaigns  modeled  in  SLAACM,  the  above 
heuristic  typically  does  not  yield  the  same  solution  for  the  attacker  as  the  IP  solu¬ 
tion  solved  to  optimality  in  LINGO.  For  the  purposes  of  SLAACM,  it  appears  that 
it  is  rare  for  the  heuristic  solutions  to  differ  by  more  than  a  few  percent  from  op¬ 
timality.  Nevertheless,  theoretical  bounds  on  such  a  heuristic  for  a  general  IP 
problem  can  be  rather  far  from  optimality,  and  those  rare  instances  of  significant 
difference  may  be  militarily  significant. 

Exact  Solution 

The  attacker’s  IP  problem  may  be  solved  exactly  using  an  integer  programming 
solver.  LINGO  is  a  commercially  available  general  linear,  nonlinear,  and  integer 
programming  solver  tool.2  The  version  embedded  in  SLAACM  is  Extended 
LINGO™  9.0,  which  does  not  limit  the  number  of  variables  (integer  or  otherwise) 
or  constraints.  Of  course,  for  a  large  enough  problem,  even  sophisticated  IP 
solvers  will  have  trouble  producing  the  optimal  solution  in  reasonable  amounts  of 
time.  The  example  problems  run  in  SLAACM,  but  rarely  take  more  than  a  couple 
of  seconds  on  a  modem  desktop  PC. 

2  Lindo  Systems,  Inc.,  http://www.lindo.com/. 


4-5 


The  IP  model  below  is  written  in  LINGO’s  modeling  language  and  implemented 
in  SLAACM  for  the  attacker’s  optimization.  The  data  and  variables  are  defined 
elsewhere  in  SLAACM.  For  sake  of  clarity,  we  omit  portions  of  the  data,  variable 
definitions,  and  input/output  declarations. 

MODEL: 

SETS: 

RPACKAGE/1  .rposize/:  rpo,  X; 

RED  AC/1  .rac/:  rob; 

RMATRIX(RPACKAGE,REDAC):  rcnst; 

ENDSETS 

!  INTEGER  PROGRAMMING  PROBLEM; 

MAX  =  VALUE; 

VALUE  =  @SUM(RPACKAGE:  rpo*X); 

@FOR(  REDAC(J): 

@SUM(RPACKAGE(I):  rcnst(l,J)*X(l))  <=  rob(J) ); 

@SUM(RPACKAGE(I):  X(l))  <=  (redLim/12); 

@FOR(  R PACKAGE:  @GIN(X)); 

END 

Determining  the  Defender’s  Response 

Certain  defender  aircraft  are  identified  by  the  user  as  smart.  They  are  assumed  to 
know  the  composition  of  the  attacker  packages  before  intercepting  them  and  to 
coordinate  their  defense  to  maximize  its  total  value  (value  computations  are  de¬ 
scribed  above  in  the  section  “Defender’s  Payoff  Function”).  The  remaining 
defender  aircraft  encounter  attacker  packages  randomly.  The  following  subsec¬ 
tions  describe  how  SLAACM  determines  the  set  of  engagements  made  by  smart 
defenders  and  by  the  remaining  defender  aircraft. 

Defender’s  Integer  Programming  Problem 

Given  the  set  of  attacker  packages,  the  smart  defenders  determine  which  types  to 
engage  by  maximizing  the  total  defender  payoff,  subject  to  constraints  imposed 
by  their  order  of  battle.  As  stated  previously,  the  defender  responds  with  optimal 
defenses  of  the  attacker  if  a  given  aircraft  is  smart  enough  to  have  a  priori  knowl¬ 
edge  of  the  attacker  package  composition.  We  assume  the  attacker  has  obtained 
its  solution  before  the  defender  calculates  its  optimal  response.  We  denote  the  op¬ 
timal  solution  for  the  attacker — that  is,  x  =  (xj  ,x2,...,Xj)  =  (a1,a2,...,aI) — as 
constraints  for  the  defender’s  optimization.  In  other  words,  at  is  the  number  of 
attacker  packages  of  type  i . 

The  defender  has  his  own  payoff  function  for  each  type  of  potential  attacker  pack¬ 
age.  We  denote  pik  as  the  payoff  to  the  defender  for  intercepting  an  attacker 
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package  i  with  the  defender  n  -ship  package  of  aircraft  type  k  (in  the  traditional 
air  defense  scenario,  n  =  4 ).  Again,  the  components  of  the  payoff  function  have 
been  described  previously;  they  include  a  loss-aversion  factor  for  the  defender, 
along  with  positive  payout  for  expected  attacker  kills  and  bomb  damage  pre¬ 
vented.  We  denote  the  set  of  defender  aircraft  types  with  a  priori  knowledge  of 
attacker  packages  as  K'  c  K .  We  denote  the  number  of  defender  aircraft  of  type 
k  e  K'  as  mk  .  We  then  formulate  and  solve  the  following  IP  problem  for  the  de¬ 
fender.  In  short,  we  want  to  solve  for  the  optimal  number  of  defender  n  -ship 
package  of  type  k  that  intercepts  attacker  package  i ,  denoted  by  yik  . 

I  K' 

max  XZp.ry.k 

i=l  J'=l 

I 

s-t-Znyik-mk  VkeK' 

i=l 


Zy^ai  vieI 

k=l 

yikE  {0,1,2,...}  [Eq.  4-3] 


The  IP  model  written  in  LINGO’s  modeling  language  is  embedded  into 
SLAACM.  LINGO  solves  the  IP  problem  and  returns  the  solution  to  variables 
defined  in  SLAACM.  At  this  point,  the  remaining  not-smart  defender  aircraft,  de¬ 
noted  by  K'cK,  where  K'  uK'=K  and  K'  n  K"  =  0 ,  engage  the  attacker 
packages  that  were  not  selected  by  the  smart  defenders.  That  is,  the  defender  air¬ 
craft  without  a  priori  knowledge  engage  randomly  the  attacker  packages  that 
remain  after  the  above  IP  problem  is  solved;  in  other  words,  those  defender  air¬ 
craft  engage  a  set  of  attacker  packages  (a[,a'2 ,  where 
f  k  h 


max 


a,  -Zy.k’° 


Those  encounters  are  described  in  a  later  section. 


V  k=l  J 


Heuristic  Solution 


As  stated  previously,  SLAACM  allows  the  user  to  specify  whether  to  use  LINGO 
to  solve  the  IP  formulation  or  to  use  a  greedy  heuristic.  If  SLAACM  uses  the  heu¬ 
ristic  (or  solves  the  IP  problem  in  LINGO),  it  must  do  so  for  both  the  attacker’s 
and  defender’s  optimization. 

The  greedy  heuristic  for  the  defender  works  similarly  to  the  heuristic  for  the  at¬ 
tacker.  Following  the  notation  from  the  IP  formulation,  we  denote  pik  as  the 

payoff  to  the  defender  for  intercepting  an  attacker  package  i  with  the  defender  n  - 
ship  package  of  aircraft  type  k  .  We  then  order  the  pik  from  greatest  to  least  over 
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all  i  and  k  .  We  denote  that  new  set  of  payoffs  as  p,  >  p2  >...  >  pB  ,  and  note 
their  associated  i  and  k  . 

It  is  worth  noting  that  defenders  that  are  not  members  of  the  local  air  defense 
(LAD)  are  presumed  to  be  loss  averse,  and  it  is  not  uncommon  for  some  pik  <  0 ; 
in  other  words,  intercepting  certain  attacker  packages  may  have  negative  payoffs 
for  the  defender.  This  result  occurs  because  loss-averse  defenders  suffer  a  large 
penalty  for  their  expected  losses.  Neither  the  heuristic  nor  the  LINGO  model  uses 
any  aircraft  whose  payoff  values  result  in  pik  <  0 ,  since  selecting  those  combina¬ 
tions  would  of  course  reduce  the  defender’s  objective  function. 

Again,  we  denote  the  number  of  attacker  packages  of  type  i  sent  as  aj .  As  be¬ 
fore,  we  denote  the  number  of  defender  aircraft  of  type  k  e  K'  as  mk  .  The  set  of 
decision  variables  (number  of  defender  n  -ship  packages  of  type  k  that  intercept 
attacker  package  i )  is  denoted  by  yik  .  As  with  the  attacker,  we  wish  to  solve  a 
multidimensional  knapsack  problem. 

We  can  then  describe  the  optimization  heuristic  we  employ  by  the  following  algo¬ 
rithm: 


Step  0:  Set  mk  =  mk  for  all  k  . 

Set  a)  =  a;  for  all  i . 

Initialize  yik  =0  for  all  i,k. 

Begin  with  b  =  1 . 

Step  1:  Let  i,k  be  the  values  for  p~  such  that  p,.  =  pb 
If  mk  -  n  >  0  and  a;  >  0 ,  Then  { 
mk  =ffik  -n 
a,  =  a,  -  1 

yik  =  yik  + 1 

Go  to  Step  1 

}  Else  { 

Go  to  Step  2} 

Step  2 :  b  =  b  + 1 . 

If  pb  <  0,  quit. 

If  b  =  B  +  1 ,  quit. 

Else,  Go  to  Step  1 . 
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In  essence,  this  heuristic  finds  the  solution  with  the  largest  payoff  to  the  defender; 
notes  the  associated  attacker  and  defender  package  types  (i.e.,  i  and  k  ,  respec¬ 
tively);  checks  to  see  if  n  defender  aircraft  of  that  type  are  available;  checks  to 
see  if  the  attacker  sent  any  packages  of  that  type;  and  chooses  those  defender  air¬ 
craft.  It  sends  the  lesser  of  the  number  of  defender  aircraft  of  that  type  still 
available  or  the  number  of  attacker  packages  of  that  type  that  were  sent.  Then,  the 
heuristic  looks  for  the  next  highest  payoff  and  repeats  the  above  process  until  all 
smart  defender  aircraft  with  positive  payoffs  are  assigned  or  until  there  are  no 
more  attacker  packages  to  which  to  assign  a  smart  defender  4-ship  package. 

This  process  is  similar  to  the  one  for  the  attacker  heuristic,  but  it  needs  to  be  re¬ 
stated  that  this  method  does  not  guarantee  optimality  for  the  defender,  because  it 
does  not  consider  the  optimization  of  the  entire  defender  fleet  over  all  possibili¬ 
ties.  Instead,  it  looks  for  a  greedy  local  solution.  That  said,  both  the  defender’s 
and  the  attacker’s  heuristics  are  computationally  easy  to  implement.  Indeed,  with¬ 
out  the  ability  to  solve  moderately  complex  integer  programs,  these  techniques 
can  form  good  approximations.  However,  neither  heuristic  is  able  to  provide  tight 
bounds  on  its  perfonnance,  at  least  for  a  general  knapsack  problem.  Therefore,  the 
preferred  method,  given  the  choice  of  the  two  methods,  is  to  solve  both  optimiza¬ 
tions  exactly  using  an  IP  solver. 

Exact  Solution 

The  defender’s  IP  problem  may  be  solved  exactly  using  LINGO.  In  general,  the 
defender’s  optimization  is  a  little  more  complex  than  the  attacker’s,  because  it 
requires  assigning  defender  packages  to  attacker  packages.  While  the  attacker 
simply  creates  optimized  attacker  packages,  the  defender  is  effectively  optimizing 
over  two  indices  (which  defender  packages  should  be  sent,  and  which  attack 
packages  should  be  intercepted).  In  most  cases,  however,  the  number  of  variables 
in  the  defender’s  optimization  is  not  larger  than  the  number  of  decision  variables 
in  the  attacker’s  optimization,  because  the  attacker  sends  only  a  small  fraction  of 
its  potential  package  types.  In  both  cases,  LINGO  obtains  the  optimal  solution  to 
the  IP  problem  in  a  few  seconds  on  a  standard,  modern  desktop  PC  for  most  sam¬ 
ple  problems  in  SLAACM. 

The  IP  model  written  below  in  LINGO’s  modeling  language  solves  the  defender’s 
optimization  problem.  It  is  similar  in  syntax  to  the  attacker’s  LINGO  model. 
Again,  the  data  and  the  variable  names  are  defined  elsewhere  in  SLAACM.  For 
the  sake  of  clarity,  we  omit  most  data  and  variable  definitions,  as  well  as  in¬ 
put/output  calls. 

MODEL: 

SETS: 

PPACKAGE/1  .bsize/:  sboob; 

REDAC/1  .rsize/:  rsent; 

PMATRIX(REDAC,  PPACKAGE):  bpo,  X; 

ENDSETS 
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!  INTEGER  PROGRAMMING  PROBLEM; 

I  Objective  Function; 

MAX  =  VALUE; 

VALUE  =  @SUM(PMATRIX:  bpo*X); 

I  subject  to; 

@FOR(  PPACKAGE(J):  @SUM(REDAC(I):  X(l,J)*ndef)  <=  sboob(J) ); 

@FOR(  REDAC(I):  @SUM(PPACKAGE(J):  X(I,J))  <=  rsent(l) ); 

@FOR(PMATRIX(l,J):  @GIN(X(I,J))); 

END 

Random  Engagements  by  Not-Smart  Defender  Aircraft 

The  optimal  engagements  made  by  smart  defender  aircraft  may  leave  some  at¬ 
tacker  packages  unengaged.  For  determining  the  packages  that  need  to  be 
intercepted  by  not-smart  defenders,  we  assume,  optimistically,  that  an  attacker 
package  intercepted  by  a  smart  attacker  flight  does  not  continue  the  engagement. 
This  assumption  is  reasonable  when  the  smart  defenders  are  much  stronger  than 
any  of  the  attacker’s  fighters  and  bombers. 

Other  SLAACM  versions  have  considered  the  opposite  case,  which  may  be  im¬ 
plemented  in  SLAACM  with  simple  reprogramming,  in  which  the  attacker  forces 
that  survive  the  smart  defense  reorganize  into  optimal  attacker  packages  and  con¬ 
tinue  the  assault. 

The  defender  aircraft  not  identified  as  smart  encounter  attacker  packages  ran¬ 
domly.  They  do  not  conduct  optimized  defensive  encounters.  There  is  still  an 
element  of  battle  management  for  these  aircraft,  however,  and  that  is  the  effi¬ 
ciency  with  which  they  can  locate  and  engage  the  attacker  packages. 

With  perfect  battle  management,  if  the  number  of  defending  flights  is  at  least  as 
large  as  the  number  of  attacker  packages,  every  attacker  package  would  be  inter¬ 
cepted.  The  intercepting  aircraft  types  would  be  random,  however,  rather  than 
optimal.  If  the  number  of  attacker  packages  is  larger  than  the  number  of  defending 
flights,  every  defending  flight  would  engage  an  attacker  package. 

With  less-than-perfect  battle  management,  not  every  attacker  package  would  be 
intercepted  in  the  former  case,  and  not  every  defending  flight  would  be  engaged 
in  the  latter  case.  SLAACM ’s  battle  management  feature  allows  the  user  to  enter  a 
“goodness”  parameter  that  characterizes  the  effectiveness  of  battle  management 
for  defender  aircraft  that  are  not  smart.  The  following  paragraphs  explain  how 
that  parameter  affects  SLAACM’s  calculations. 

Let  K  distinct  types  of  defender  flights  deal  with  an  attack  by  J  distinct  types  of 
attacker  packages.  The  defenders  are  assumed  not  to  know  the  makeup  of  indi¬ 
vidual  attacker  packages  before  interception,  so  that  the  type  of  attacker  package 
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engaged  by  a  given  defending  flight  is  the  result  of  random  selection  from  the  set 
of  attacker  packages. 

Let  m;  be  the  number  of  defender  flights  of  type  i,  and  let  nj  be  the  number  of  at¬ 
tacker  packages  of  type  j .  Then  the  total  number  of  attacker  packages,  N,  and  the 
total  number  of  defending  flights,  M,  are  given  respectively  by 

N  =  i>j;  M  =  £mi.  [Eq.  4-4] 

1  1 

First  we  consider  the  case  of  perfect  battle  management.  If  M  >  N,  then  every  at¬ 
tacker  package  will  be  intercepted.  Not  all  defending  flights  engage;  in  this  simple 

3 

analysis,  we  assume  that  the  fraction  N/M  of  each  defending  flight  type  engages. 

Then,  a  measure  of  Ey ,  the  central  tendency  of  the  number  Ey  of  mj  vs.  nj  en¬ 
gagements,  is 


N  nj  minj 

E:;  =  — rri:  —  = - 

1]  M  1  N  M 


[Eq.  4-5] 


For  this  simple  analysis,  we  take  E-  to  be  the  number  of  m,  vs.  nj  engagements. 

Summing  E  -  over  j  shows  that  the  fraction  of  each  defending  flight  type  engaged 

is  N/M,  as  it  should  be;  summing  over  i  shows  that  all  attacking  flights  of  each 
type  are  engaged. 

When  M  <  N,  every  defending  flight  engages,  but  not  all  attacker  packages  can  be 
engaged.  The  estimate  for  E-  analogous  to  the  one  given  in  Equation  4-5  is 


_  ri:  mm: 

E;;  =  ' 


1  N 


N 


[Eq.  4-6] 


Summing  this  estimate  for  E-  over  i  shows  that  the  fraction  M/N  of  each  attacker 

package  type  is  engaged;  summing  over  j  shows  that  every  defending  flight  of 
each  type  is  engaged. 

The  Battle  Management  Efficiency  Factor 

In  some  cases,  M  and  N  are  both  0(10“).  For  such  large  M  and  N,  providing  this 
“perfect”  defender’s  battle  management  could  exceed  the  capabilities  of  available 
systems.  We  account  for  the  limitations  of  the  defender’s  battle  management  with 
a  simple  adaptation  of  the  “perfect”  case.  We  make  the  adaptation  by  introducing 
“ghost”  attacker  packages  or  defender  flights.  A  defender  flight  that  engages  a 


3  We  recognize  that  this  simple  analysis  treats  a  complex  combinatoric  problem  crudely. 
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ghost  attacker  package  does  not  actually  engage;  an  attacker  package  engaged  by 
a  ghost  defender  flight  is  not  actually  intercepted. 

When  M  >  N,  we  introduce  mg  defending  flights  of  type  “ghost,”  and  proceed  as 
in  the  above  analysis  (where  the  defender  has  perfect  battle  management).  This 
gives  us  a  new  estimate  forE^  : 


_  m=n;  m.n: 

E  = _ i _ I _ = _ - _ Lb 

IJ  M  +  mg  M  1 


[Eq.  4-7] 


where 


bi  = 


M 


M  +  mg 

When  M  <  N,  we  introduce  ng  attacker  packages  of  type  “ghost,”  and  find 


[Eq.  4-8] 


_  mm:  m.ri: 

E  —  1  J  _  1  J 

ij  ~  N  +  n„  “  M 


[Eq.  4-9] 


where 


b2  =  — — — .  [Eq.  4-10] 

~  N  +  ng 

Noting  the  similarity  of  Equations  4-8  and  4-10,  we  introduced  into  SLAACM  the 
simple  model 


-  min> 

Eii=^rb  [Eq-4-11] 

and  allow  the  user  to  choose  the  “battle  management  efficiency  factor” 
b,0  <  b  <  1. 

Once  seen,  Equation  4-11  is  so  simple,  and  seemingly  obvious,  that  it  is  fair  to  ask 
why  we  did  not  simply  introduce  it  into  SLAACM  without  analysis.  The  reason  is 
that  we  wanted  to  understand  what  that  simple  choice  implied  about  how  the  bat¬ 
tle  proceeded.  The  analysis  given  here  supplies  that  understanding. 
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Chapter  5 

How  SLAACM  Works: 

SLAACM  Campaign  Calculations 


This  chapter  describes  the  calculation  protocols  used  by  SLAACM  for  campaign 
analysis. 

Determining  and  Propagating 
Loss  and  Destruction  Statistics 


The  calculations  described  in  Chapter  4,  under  “Determining  Attacks  and  De¬ 
fenses,”  generate  the  sets  of  engagements  that  take  place  in  each  day’s  fighting. 
There  are  then  specified  numbers  of  engagements  between  a  specific  attacker 
package  (with  its  specified  aircraft  types  for  the  advanced  escorts,  close  escorts, 
SEAD  aircraft,  and  bombers)  and  a  flight  of  2,  4,  or  8  defenders  of  a  specific  type, 
and  SAMs.  SLAACM  uses  the  method  described  in  Chapter  2,  under  “Composite 
Engagements,”  to  generate  probabilities  for  all  outcome  states  of  each  kind  of  en¬ 
gagement  that  occurs. 

Then,  to  evaluate  attacker  losses,  SLAACM  calculates  expected  values  of  the 
losses  for  advanced  escorts,  close  escorts,  SEADs,  and  bombers.  These  expected 
values  are  accumulated  for  each  attacker  aircraft  type,  for  engagements  between 
smart  defenders  and  the  attacker  packages,  and  for  engagements  between  non¬ 
smart  defenders  and  the  attacker  packages  not  engaged  by  smart  defender  aircraft. 

The  total  expected  losses  for  each  attacker  type  are  rounded  to  integers  and  sub¬ 
tracted  from  the  attacker  order  of  battle.  SLAACM  then  generates  the  attacker 
order  of  battle  for  the  following  day  using  the  surviving  attackers  plus  any  rein¬ 
forcements  scheduled  to  arrive  for  that  day.  Reinforcements  are  a  user  input  to 
SLAACM. 

Expected  values  of  surviving  bombers  are  multiplied  by  bomber  payloads  and 
rounded  to  the  nearest  integer,  to  give  tons-of-bombs  results.  Tons  delivered  by 
bombers  with  relative  circular  errors  probable  (CEPs)  less  then  a  user-definable 
level  are  accounted  for  separately  from  tons  delivered  by  other  bombers  and  iden¬ 
tified  as  “smart  bombs.” 

In  treating  defender  losses,  SLAACM  considers  measures  of  dispersion  as  well  as 
measures  of  central  tendency.  The  probability  distributions  for  outcome  states  of 
the  engagements  occurring  each  day  are  used  to  calculate  both  means  and  vari¬ 
ances  of  the  losses  of  each  type  of  defender  aircraft,  for  each  day.  Mean  losses, 
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rounded  to  the  nearest  integer,  are  subtracted  from  the  defender’s  order  of  battle 
and,  after  increments  from  any  scheduled  reinforcements,  are  posted  as  the  de¬ 
fender’s  order  of  battle  for  the  next  day. 

Variances  in  the  losses  of  each  type  of  defender  aircraft  are  accumulated  day-by¬ 
day,  and  the  total  variances  are  reported  to  give  a  measure  of  dispersion  in  de¬ 
fender  losses. 

This  procedure  neglects  a  potentially  important  source  of  dispersion,  the  effect  of 
propagating  losses  other  than  the  means  of  both  attacker  and  defender  losses.  The 
present  dispersion  calculations  are  used  for  checking  the  validity  of  using  only 
mean  values.  Propagating  even  a  limited  number  of  results,  such  as  maximum, 
mean,  and  minimum  losses,  leads  to  explosive  growth  of  computations  when  all 
combinations  are  considered.  There  may  value  to  propagating  only  the  strings  of 
worst  cases  and  best  cases  to  establish  bounds  for  the  mean  results,  especially 
when  forces  are  evenly  matched.  For  the  cases  in  which  defender  aircraft  are  al¬ 
ways  much  stronger  than  the  attacker  aircraft,  there  is  little  dispersion  in  attacker 
losses — nearly  all  attacker  aircraft  that  are  engaged  are  lost,  nearly  all  the  time — 
so  the  significant  dispersion  in  each  day  is  the  dispersion  in  defender  losses.  But, 
again  by  virtue  of  the  defender’s  much  greater  strength,  the  dispersions  never  lead 
to  significant  probabilities  of  the  defender  losing  substantial  fractions  of  his 
forces,  except  near  the  ends  of  campaigns  in  which  the  defender  is  losing.  The 
converse  is  true  for  strong  attackers. 

SLAACM’s  Sequence  of  Calculations 


This  section  describes  SLAACM’s  step-by-step  calculations.  It  explains  how 
SLAACM  uses  the  methods  described  above  to  model  air-to-air  campaigns. 

User  Inputs 

The  SLAACM  user  inputs  attacker  and  defender  orders  of  battle  for  the  first  day 
of  fighting  and  lists  any  scheduled  reinforcements. 1  The  following  subsections 
summarize  the  inputs,  focusing  on  their  meanings  and  use.  (Chapter  6  details  the 
methods  for  entering  data  into  SLAACM.) 

Inputs  for  the  Defender  Force 

The  user  may  identify  certain  defender  aircraft  as  smart.  Smart  defender  aircraft 
can  determine  the  composition  of  attacker  packages  before  intercepting  them,  and 
can  coordinate  their  defensive  response  among  themselves  to  provide  optimal 
defenses. 


1  “Day”  is  used  to  signify  a  set  of  attack  and  defense  sorties.  There  could  be  more  than  one  set 
of  sorties  on  a  single  calendar  day,  or  there  could  be  gaps  of  calendar  days  between  attacks. 
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Defender  aircraft  based  long  distances  from  the  battle  space  will  have  to  be  ra¬ 
tioned  and  rotated  to  maintain  continuous  defensive  coverage.  The  user  may 
select  a  combat  air  patrol  (CAP)  factor  for  any  defender  aircraft  type,  so  that  only 
a  specified  fraction  of  the  inventory  of  that  defender  type  is  available  to  meet  an 
attack. 

The  user  may  identify  certain  defender  aircraft  as  local  air  defenders  (LADs). 
SLAACM  treats  LADs  differently  from  other  defender  types  when  considering 
defender  dispatch  options.  Normally,  the  defender  values  the  life  of  its  forces  and 
will  not  dispatch  a  defender  type  if  engaging  a  attacker  package  gives  a  negative 
value  of  the  defender’s  payoff  function.  Defender  types  designated  as  LADs, 
however,  will  be  dispatched  in  all  cases. 

Inputs  for  the  Attacker  Force 

The  user  inputs  the  bomb  payload  of  each  attacker  aircraft  and  specifies  whether 
the  aircraft  is  to  be  used  as  a  fighter  only,  or  as  a  bomber  only,  or  as  either  a 
fighter  or  a  bomber.  The  user  also  inputs  the  ratio  of  the  CEP  with  which  bomber 
aircraft  can  deliver  bombs,  to  the  reference  CEP.  As  discussed  in  Chapter  3,  the 
default  ratio  is  1.0  for  dumb  bombs  and  0.1  to  0.4  for  smart  bombs.  If  delivery  of 
cruise  missiles  is  of  interest,  the  user  inputs  the  number  of  cruise  missiles  that  can 
potentially  be  carried  by  each  aircraft. 

The  attacker’s  order  of  battle  may  be  so  large  that  airfield  capacity  and  battle 
management  capability  may  preclude  using  all  forces.  The  user  may  specify  a 
maximum  number  of  attacker  aircraft  that  can  be  dispatched  in  a  given  attack. 

The  attacker  optimizes  its  dispatch  choices  with  the  assumption  that  each  package 
encounters  a  specified  defender  fighter  type.  The  user  identifies  this  type. 

Input  for  Relative  Capabilities  of  Attacker  and  Defender  Aircraft 

Currently,  SLAACM  models  engagements  with  the  M  vs.  N  probabilistic  en¬ 
gagement  model,  and  it  assumes  that  outcome  probabilities  have  their  long-time 
limiting  values  as  described  in  Chapter  2  under  “Calculations  for  the  Long-Time 
Limit.”  Thus  only  one  aircraft  parameter — the  ratio  k  of  the  parameter  of  the  de¬ 
fender  time-between-kill  distribution  to  that  parameter  for  the  attacker  time- 
between-kill  distribution — affects  outcomes. 

Values  of  k  parameters  may  not  be  readily  available,  and  SLAACM  includes  a 
utility  to  calculate  k  from  the  defender-to-attacker  “loss  ratio,”  that  is,  the  gener¬ 
ally  available  output  from  simulations  or  combat.  (Some  simulations  provide 
exchange  ratios,  which  are  ratios  of  defender  relative  losses  to  attacker  relative 
losses.  Multiplying  an  exchange  ratio  by  the  ratio  of  the  initial  number  of  defend¬ 
ers  to  the  initial  number  of  attackers  gives  the  loss  ratio.)  The  user  inputs  the  loss 
ratios  and  the  source  scenario,  such  as  4  defenders  vs.  8  attackers  with  defenders 
breaking  away  after  2  losses.  SLAACM  then  determines  the  value  of  k  that  gives 
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the  observed  loss  ratio  and  uses  it  in  subsequent  calculations.  SLAACM  uses  bi¬ 
section  to  find  k  values  correct  to  four  significant  figures. 

With  all  these  inputs  completed,  the  user  enters  the  number  of  days  to  run  the 
campaign,  specifies  the  battle  management  efficiency  parameter,  as  described  in 
Chapter  4  under  “The  Battle  Management  Efficiency  Factor,”  and  instructs 
SLAACM  to  use  either  heuristic  optimization  or  exact  optimization  with  LINGO. 
SLAACM  then  performs  the  sequence  of  calculations  described  below. 

Steps  in  SLAACM’s  Calculations 

SLAACM  first  calculates  values  of  the  attacker  payoff  function,  as  described  in 
Chapter  4  under  “Attacker’s  Payoff  Function,”  for  all  possible  choices  of  ad¬ 
vanced  escorts,  close  escorts,  and  bombers  that  can  make  up  an  attacker  assault 
package.  There  will  often  be  hundreds  of  such  choices.  SLAACM  stores  the  pay¬ 
off  values  for  use  in  later  optimizations. 

SLAACM  then  begins  the  day-by-day  calculations  of  outcomes  in  the  campaign. 
Using  either  the  heuristic  or  the  exact  solution,  as  specified  by  the  user,  SLAACM 
determines  the  set  of  attacker  packages  that  maximizes  the  total  attacker  payoff 
function,  subject  to  constraints  imposed  by  the  attacker’s  available  forces  and  to 
the  maximum  number  of  attacker  aircraft  that  can  be  dispatched  in  1  day  as  speci¬ 
fied  by  the  user.  These  calculations  are  described  in  Chapter  4  under 
“Determining  Attacks  and  Defenses.” 

SLAACM  considers  this  set  of  attacker  packages  to  determine  the  optimal  defense 
by  the  smart  defender  aircraft.  After  calculating  the  defender  payoff,  as  described  in 
Chapter  4  under  “Defender’s  Payoff  Function,”  for  each  possible  combination  of 
smart  defender  defense  with  an  attacker  package  type  actually  dispatched, 
SLAACM  determines  the  defender  dispatch  option  that  maximizes  the  total  de¬ 
fender  payoff,  constrained  by  the  available  defender  forces  and  the  CAP  values 
input  by  the  user.  SLAACM  uses  either  the  heuristic  or  the  exact  solution  of  the  re¬ 
sulting  IP  problem,  as  directed  by  the  user. 

As  described  above  under  “Determining  and  Propagating  Loss  and  Destruction 
Statistics,”  SLAACM  then  evaluates  the  losses  to  all  forces  from  engagements 
between  the  attacker  packages  and  the  smart  defenders.  Then  SLAACM  calcu¬ 
lates  the  results  of  engagements  between  the  attacker  packages  that  were  not 
engaged  by  smart  defender  aircraft.  Only  those  attacker  packages  not  engaged  by 
smart  defenders  are  available  for  engagements  with  the  not-smart  defender  air¬ 
craft. 

SLAACM  calculates  the  defender  payoff  for  each  combination  of  not-smart  de¬ 
fender  aircraft  and  a  attacker  package  type  that  was  not  intercepted  by  a  flight  of 
smart  defenders.  Then,  following  the  procedure  described  in  Chapter  4  under 
“Random  Engagements  by  Not-Smart  Defender  Aircraft,”  and  considering  the 
user’s  input  of  battle  management  efficiency  factor,  SLAACM  determines  the 
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numbers  of  engagements  between  the  several  types  of  not-smart  defenders  and  the 
attacker  package  types  available  to  them.  As  described  above  under  “Determining 
and  Propagating  Loss  Statistics,”  SLAACM  updates  the  losses  of  all  aircraft 
types,  calculates  the  tons  of  smart  and  not-smart  bombs  delivered,  and,  in  view  of 
any  scheduled  reinforcements,  updates  attacker  and  defender  orders  of  battle  for 
the  next  day  of  the  campaign. 

After  calculating  outcomes  for  the  number  of  days  specified  by  the  user, 
SLAACM  produces  daily  records  of  attacker  and  of  defender  orders  of  battle,  of 
attacker  losses,  of  defender  losses  and  their  variances  due  to  engagements,  of  tons 
of  smart  and  of  not-smart  bombs  delivered,  and  of  cruise  missiles  delivered.  As 
described  in  Chapters  6  and  7,  SLAACM  also  produces  graphs  of  attacker  and  of 
defender  drawdowns,  the  quantities  and  types  of  bombs  delivered,  and  informa¬ 
tion  on  “who  shot  whom.” 
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Chapter  6 

Users  Guide  to  Classic  SLAACM 


Classic  SLAACM  is  housed  in  a  Microsoft  Excel  workbook  with  extensive  Visual 
Basic  for  Applications  code.  The  current  version  of  Classic  (Defender)  SLAACM 
was  updated  to  use  object-oriented  programming  principles.  In  addition,  it  has  an 
option  to  conduct  integer  programming  by  linking  to  an  external  optimization  pro¬ 
gram  called  LINGO. 1  Using  Excel  to  house  SLAACM  has  several  advantages. 
Most  analysts  are  familiar  with  the  Excel  environment  and  will  be  able  to  navi¬ 
gate  it  easily.  Other  than  acquiring  the  SLAACM  workbook,  users  do  not  need  to 
install  additional  software  unless  they  want  to  run  integer  programming  optimiza¬ 
tions  in  LINGO.  Finally,  users  can  easily  copy  and  paste  output  values  and  charts 
into  their  preferred  presentation  formats  for  displaying  results  to  others. 

This  chapter  describes  how  to  operate  Classic  SLAACM.  It  begins  with  a  descrip¬ 
tion  of  the  model  inputs.  It  then  describes  how  to  run  the  model  once  all  of  the 
input  parameters  are  set  up.  Next,  the  chapter  describes  the  outputs  reporting 
SLAACM’s  results.  The  last  section  addresses  how  the  outputs  are  analyzed  to 
determine  the  effects  of  the  various  engagement  scenarios. 


Inputs 


SLAACM’s  input  parameters  are  defined  in  five  worksheets: 

♦  BlueSupply,  which  is  used  to  specify  the  quantities  of  Blue  aircraft  avail¬ 
able  on  each  day  of  the  campaign. 

♦  RedSupply,  which  is  used  to  specify  the  quantities  of  Red  aircraft  avail¬ 
able  on  each  day  of  the  campaign. 

♦  BlueOOB,  which  is  used  to  define  Blue’s  order  of  battle  and  other  perti¬ 
nent  parameters.  BlueOOB  is  the  main  sheet  from  which  SLAACM  runs 
are  conducted.  (SLAACM  also  has  a  RedOOB  worksheet,  but  it  is  purely 
an  output  worksheet  and  does  not  receive  any  user  inputs.) 

♦  ExRatios,  which  is  used  to  specify  engagement  parameters  such  as  loss 
ratios  and  breakpoints. 

♦  CM,  which  is  used  to  specify  the  number  of  cruise  missiles  that  a  Red  air¬ 
craft  could  carry. 

We  discuss  each  sheet  below  and  describe  the  proper  way  to  complete  each  one. 

1  For  more  information  on  LINGO,  see  http://www.lindo.com/products/lingo/. 
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BlueSupply  and  RedSupply  Worksheets 

For  each  force,  Blue  and  Red,  the  user  will  need  to  set  up  the  initial  quantities  of 
aircraft  available  and  the  reinforcements  available  on  each  day  of  the  campaign. 
The  user  will  need  to  type  those  quantities  into  the  appropriate  worksheets,  called 
BlueSupply  and  RedSupply. 

Figure  6-1  is  an  example  of  the  BlueSupply  worksheet.  As  shown  in  the  figure, 
the  worksheet  has  a  row  for  each  Blue  aircraft  type  and  a  column  for  each  day  of 
the  campaign.  Initial  aircraft  quantities  are  entered  into  the  column  labeled  0.  Re¬ 
inforcements  may  be  provided  for  subsequent  days  of  the  campaign.  In  the  figure, 
we  see  that  reinforcements  are  entered  for  Day  2  and  Day  3. 2 

Figure  6-1.  Classic  SLAACM:  BlueSupply  Worksheet 


The  RedSupply  worksheet  is  designed  in  the  same  way  as  the  BlueSupply  sheet. 
Figure  6-2  is  an  example.  In  the  figure,  we  see  that  only  initial  aircraft  quantities 
are  provided  and  no  reinforcements  are  specified. 


2  Note  that  a  day  is  used  in  the  model  as  a  convenient  designator  to  represent  an  individual  at¬ 
tack/defense  sortie.  In  an  actual  campaign,  there  may  be  more  than  one  sortie  on  a  particular  cal¬ 
endar  day,  or,  conversely,  sorties  may  be  distributed  among  several  calendar  days. 
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Figure  6-2.  Classic  SLAACM:  Red  Supply  Worksheet 


The  BlueSupply  and  RedSupply  worksheets  can  accommodate  up  to  20  aircraft 
types  each.  If  new  rows  are  needed  to  accommodate  additional  aircraft  types,  they 
can  added  using  the  “Add  a  New  Row”  button. 


Caution! 

Any  time  a  row  is  added  or  deleted,  it  is  imperative  that  the  buttons  on  the 
BlueSupply  and  RedSupply  worksheets  are  used  to  add  or  delete  the  rows. 

Do  not  simply  insert  or  delete  rows  and  columns. 


Several  of  SLAACM’s  worksheets  reference  the  aircraft  types.  To  keep  the  work¬ 
sheets  synchronized,  the  “Add  a  New  Row”  and  “Delete  a  Row”  buttons  must  be 
used.  Using  typical  worksheet  approaches  to  inserting  and  deleting  rows  will 
likely  result  in  model  errors  or  misinterpretation  of  results. 

BlueOOB  Worksheet 

The  worksheet  named  BlueOOB  (which  stands  for  Blue  order  of  battle),  is  the 
main  worksheet  for  input  parameters  for  the  Blue  aircraft.  It  is  also  the  worksheet 
used  to  launch  campaign  analyses.  BlueOOB  also  serves  as  an  output  worksheet, 
as  does  the  RedOOB  worksheet;  the  output  portion  of  BlueOOB  will  be  discussed 
in  the  output  section  of  this  chapter.  Figure  6-3  shows  the  BlueOOB  worksheet. 
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Figure  6-3.  Classic  SLAACM:  BlueOOB  Worksheet 


Caution! 

Rows  or  columns  should  never  be  deleted  from  this  worksheet! 


The  aircraft  listed  in  the  BlueOOB  worksheet  are  linked  to  those  that  were  shown 
in  the  BlueSupply  worksheet  (Figure  6-1).  The  Day  0  quantities  are  also  linked  to 
those  in  the  BlueSupply  worksheet.  The  user  who  attempts  to  click  the  aircraft 
names  or  Day  0  quantities  on  the  BlueOOB  sheet  will  receive  the  warning  shown 
in  Figure  6-4. 

Figure  6-4.  Warning:  Do  Not  Alter  Aircraft  of  BlueOOB  Worksheet 


The  first  four  columns  on  the  BlueOOB  worksheet  define  aircraft-specific  input 
parameters  for  Blue: 

♦  “Smart?  ”  column.  The  user  may  use  this  column  to  identify  specific  Blue 
aircraft  as  smart.  Smart  Blue  aircraft  can  determine  the  composition  of 
Red  attack  packages  before  intercepting  them  and  can  coordinate  their 
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defensive  response  among  themselves  to  provide  optimal  defenses.  This 
coordination  and  optimization  is  discussed  in  Chapter  5.  Multiple  Blue  air¬ 
craft  types  may  be  designated  as  smart.  In  Figure  6-3,  the  first  two  aircraft 
are  selected  to  be  smart. 

♦  “LAD?  ”  column.  The  user  may  use  this  column  to  identity  certain  Blue 
aircraft  as  local  air  defenders.  SLAACM  treats  LADs  differently  from 
other  Blue  types  when  considering  Blue  dispatch  options.  In  general,  a 
non-LAD  Blue  has  a  penalty  for  its  own  losses  that  is  significantly  greater 
than  a  LAD  Blue,  which  presumably  is  less  loss  averse  due  to  defending 
its  own  territory.  Some  non-LAD  Blues  may  have  a  negative  payoff  func¬ 
tion  for  engaging  certain  Red  packages  dispatched.  Blue  will  not  dispatch 
that  type.  LADs,  however,  are  dispatched  in  all  cases.  Multiple  Blue  air¬ 
craft  types  may  be  designated  as  LADs.  In  Figure  6-3,  the  eighth  aircraft  is 
designated  as  such.  Note  that  SLAACM  also  has  a  separate  aircraft  type 
labeled  BlueLAD,  which  is  intended  to  represent  an  indigenous  “Green” 
defender.  In  earlier  versions  of  SLAACM,  this  Blue  LAD  was  the  only 
aircraft  in  the  model  that  placed  no  value  on  its  own  survival.  The  current 
version  allows  all  aircraft  to  be  optionally  “fearless.”  This  facilitates  the 
analysis  of  supplying  different  types  of  aircraft  to  the  indigenous  forces. 

♦  “CAP  Factor”  column.  This  factor  addresses  the  situations  in  which  Blue 
fighters  are  stationed  at  remote  airfields  and  must  divide  their  forces  to 
maintain  continuous  combat  air  patrol  coverage  over  the  battle  space.  The 
user  may  select  a  CAP  factor  for  any  Blue  aircraft  type,  so  that  only  a 
specified  fraction  of  the  inventory  of  that  Blue  type  is  available.  CAP  fac¬ 
tor  values  are  chosen  from  a  drop-down  box  for  each  aircraft  type.  Values 
range  from  1  to  10.  A  CAP  factor  of  1  indicates  that  no  CAP  limitation  is 
placed  on  that  aircraft  type.  In  Figure  6-3,  several  aircraft  have  CAP  fac¬ 
tors  greater  than  1.  For  example,  the  first  aircraft  has  been  assigned  a  CAP 
factor  of  4.  That  aircraft,  BlueFl,  also  begins  with  an  initial  quantity  of 
24  aircraft.  Therefore,  Blue_Fl  has  only  24/4  =  6  available  aircraft  on  the 
first  day,  which  means  that  only  one  4-ship  package  of  BlueFls  can  be 
deployed.  In  addition,  if  Blue  Fl  losses  on  subsequent  campaign  days 
were  to  exceed  8  aircraft,  Blue  could  no  longer  dispatch  Blue  Fls,  since 
15/4  <  4.  Thus,  Blue  Fl  would  be  unable  to  maintain  a  4-ship  package  at 
all  times.  This  assumption  is  fairly  conservative:  if  there  are  not  enough 
aircraft  to  fulfill  the  combat  air  patrol  at  all  times,  those  flights  are  not  dis¬ 
patched. 

♦  “ Index  of  Blue  a/c  for  Red’s  planning ”  column.  Only  one  aircraft  type  can 
be  selected  for  this  parameter.  In  Figure  6-3,  the  first  aircraft  type  is  des¬ 
ignated  as  Red’s  “planning  aircraft.”  This  aircraft  type  is  used  by  Red  to 
compute  its  expected  payoff  function  and  optimize  its  dispatch  choices. 
Experience  has  shown  that  Red  optimization  is  relatively  insensitive  to  the 
selection  of  the  Blue  type,  so  nonnal  practice  is  to  select  the  most  numer¬ 
ous  Blue  type  for  Red  planning. 
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In  addition  to  the  aircraft-specific  parameters,  several  other  input  parameters  must 
be  selected  on  the  BlueOOB  worksheet. 


At  the  top  of  the  worksheet,  the  number  of  days  in  the  campaign  is  specified.  In 
Figure  6-3,  the  number  of  campaign  days  is  8.  The  longer  the  campaign,  the 
longer  it  will  take  SLAACM  to  run,  since  each  campaign  day  requires  the  full 
breadth  of  SLAACM  computations.  Initially,  the  user  may  need  to  start  out  with  a 
longer  campaign  and  then  decrease  its  length  as  a  steady  state  in  aircraft  losses  is 
reached.  Eight  is  a  good  initial  number  of  days  to  use  since  some  campaigns  will 
complete  earlier  and  some  later,  depending  on  the  mixes  and  capabilities  of  Blues 
and  Reds  specified. 

Beneath  the  number  of  campaign  days  is  a  parameter  labeled  “BM  r|.”  This 
parameter  is  the  battle  management  efficiency  factor.  BM  r)  is  a  value  greater  than 
0  and  less  than  or  equal  to  1 .  Simply  stated,  a  BM  r)  less  than  1  indicates  the  per¬ 
centage  of  time  that  a  Blue  defender  will  be  unable  to  find  a  Red  attack  package. 
In  Figure  6-3,  the  efficiency  factor  is  shown  as  0.85,  which  means  that  15  percent 
of  the  not-smart  Blue  defenders  will  chase  “ghost”  attackers.  No  losses  will  result, 
but  in  those  cases,  Blue  will  not  be  confronting  an  actual  attacker,  so  its  efficiency 
is  reduced.  Since  smart  Blues  know  Red’s  attack  package  compositions  in  ad¬ 
vance,  we  assume  they  have  perfect  management. 

Next  to  the  campaign  days  and  the  efficiency  factor  is  a  drop-down  box  labeled 
“LINGO?”  If  “No”  is  selected,  as  shown  in  Figure  6-3,  SLAACM  will  use  its  in¬ 
ternal  heuristic  optimization  to  calculate  engagement  results.  If  “Yes”  is  selected, 
LINGO  will  be  used  to  perfonn  integer  optimization  (a  licensed  copy  of  LINGO 
software  is  required  to  be  installed  on  the  PC  for  this  feature  to  work).  Both  the 
heuristic  and  the  exact  integer  programming  solution  are  described  in  detail  in 
Chapter  5.  If  LINGO  is  used,  there  are  two  LINGO  input  worksheets:  “RedOpt” 
and  “BlueOpt.”  However,  those  sheets  are  already  complete  and  should  not  be 
altered  by  the  user. 

Below  the  parameters  just  described  are  four  value  boxes: 

♦  Red  value  of  Blue  loss 

♦  Blue  value  of  Blue  loss 

♦  Red  value  of  Red  loss 

♦  Blue  value  of  Red  loss. 

These  parameters  establish  the  value  of  killing  an  opponent,  as  well  as  the  penalty 
a  side  takes  for  losing  one  of  its  own.  These  values,  which  are  applied  to  fighter 
aircraft  only,  are  used  in  the  optimization  objective  function.  The  Red  value  of 
Blue  loss  represents  the  benefit  to  Red  of  destroying  a  Blue  aircraft,  while  Blue 
value  of  Blue  loss  is  the  penalty  to  Blue  for  losing  one  of  its  own  aircraft.  Values 
can  be  between  zero  and  100. 
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In  the  case  shown  in  Figure  6-3,  we  see  that  the  value  of  a  Blue  loss  to  Red  is  20, 
while  the  penalty  to  Blue  for  losing  its  own  aircraft  is  40.  Therefore,  Blue’s  aver¬ 
sion  to  a  loss  is  twice  the  value  of  Red’s  benefit  for  obtaining  the  Blue  loss.  Con¬ 
versely,  we  see  in  this  example  that  Red  is  not  penalized  for  a  Red  fighter  loss 
(the  value  is  zero)  and  Blue  barely  achieves  a  benefit  for  a  Red  fighter  loss  (the 
value  is  one). 

Another  parameter,  shown  on  the  far  right  side  of  Figure  6-3,  is  labeled  “Escaping 
bombers  complete  mission?”  The  drop-down  box  allows  for  either  a  yes  or  no  se¬ 
lection.  If  yes  is  chosen,  then  bombers  that  escape  attack  go  on  to  complete  their 
mission  and  deliver  bombs  to  their  target.  If  no  is  chosen,  then  escaping  bombers 
do  not  complete  their  missions. 

The  final  input  selection  on  BlueOOB  is  a  drop-down  box  labeled  “Alerts.”  The 
default  setting  for  this  feature  is  “Off.”  However,  if  “On”  is  selected,  message 
boxes  will  become  visible  as  the  analysis  runs. 

ExRatios  Worksheet 

The  ExRatios  worksheet  contains  a  row  for  each  Red  aircraft  that  has  been  desig¬ 
nated  on  the  RedSupply  worksheet.  The  worksheet  has  several  parameter  col¬ 
umns,  discussed  below,  and  a  column  for  each  Blue  aircraft  labeled  “F  v.  Bluex” 
where  “F”  indicates  that  Red  is  a  fighter  and  “x”  corresponds  to  the  label  provided 
on  BlueSupply.  Columns  are  repeated  for  each  Blue  aircraft,  except  they  are  la¬ 
beled  “B  v.  Blue  x”  where  “B”  indicates  that  Red  is  a  bomber  and  “x”  is  as  just 
described.  The  default  setting  in  SLAACM  assumes  that  loss  ratios  are  for  en¬ 
gagements  between  4  Blues  and  8  Reds,  with  the  Reds  and  Blues  fighting  to  anni¬ 
hilation.  However,  those  assumptions  can  be  changed  on  this  sheet  under  “Red 
Start,”  “Blue  Start,”  “Red  Quit,”  and  “Blue  Quit.”  Red  Start  and  Blue  Start  corre¬ 
spond  to  the  engagement  type,  e.g.,  8  and  4,  respectively.  Red  Quit  and  Blue  Quit 
denote  the  breaking  points,  e.g.,  0  and  0  means  that  both  sides  will  fight  to  annihi¬ 
lation.  The  user  will  want  to  specify  values  that  correspond  to  the  conditions  un¬ 
der  which  the  loss  ratios  were  obtained. 

To  keep  the  rows  and  columns  synchronized  with  the  BlueSupply  and  RedSupply 
worksheets,  it  is  imperative  that  no  rows  or  columns  on  ExRatios  be  inserted  or 
deleted. 


Caution! 

As  stated  earlier,  use  the  add  and  delete  buttons  on  the  BlueSupply 
and  RedSupply  worksheets  to  add  or  delete  rows. 
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Figure  6-5  is  a  screenshot  of  the  ExRatios  worksheet. 


Figure  6-5.  Classic  SLAACM:  ExRatios  Worksheet 
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Only  red-colored  values  are  inputs  on  this  sheet. 

Fighter-Bombers  (Usage  =  0)  engage  Blue  A/C.  even  when  going  as  Bombers.  Values  posted  on  this  sheet  represent  loss-ratios  for  those 
engagements. 

Bombers-Only  (Usage  =  -1 )  do  not  engage  Blue  A/C,  and  attempt  to  escape.  Values  posted  on  this  sheet  represent  the  ratio  of  (mean  time  for 
Red  to  escape  /  (mean  time  for  Blue  to  kill). 
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The  first  column  on  the  ExRatios  worksheet  is  labeled  “Usage.”  For  each  Red  air¬ 
craft  type,  a  value  of  1,  -1,  or  0  needs  to  be  entered.  A  value  of  1  indicates  a 
fighter.  When  bombers  are  designated  by  the  Usage  =  -1  code,  they  are  assumed 
to  try  to  escape  rather  than  engage  the  Blue  fighters.  Because  the  bombers  are  es¬ 
caping  and  not  engaging  in  combat,  the  values  for  fighter  vs.  bomber  encounters 
in  the  ExRatios  spreadsheet  no  longer  represent  loss  ratios.  Rather,  they  represent 
the  ratio  of  the  mean  time  for  Red  to  escape  to  the  mean  time  for  Blue  to  kill. 
Bombers  with  the  Usage  =  0  code  are  assumed  to  be  fighter  bombers  and  will  en¬ 
gage;  their  encounter  values  in  the  ExRatios  spreadsheet  represent  kill-rate  ratios. 
(A  note  reminding  the  user  of  this  difference  has  been  added  to  the  ExRatio 
spreadsheet  and  is  shown  in  Figure  6-5.) 

The  second  column  indicates  the  name  of  the  Red  aircraft  type;  it  is  linked  to  the 
RedSupply  worksheet.  The  third  column,  “Payload,  lb,”  is  used  to  enter  the  num¬ 
ber  of  pounds  of  bombs  that  the  associated  Red  aircraft  can  carry. 

The  fourth  column  is  labeled  “cepr,”  which  stands  for  the  circular  error  probable 
ratio  (CEP/CEPref).  The  cepr  is  a  value  greater  than  0  and  less  than  or  equal  to  1, 
where  1  indicates  that  the  bomber  cannot  deliver  smart  bombs  (guided  munitions) 
and  0.2  (or  less)  indicates  that  it  can.  The  smart  bomb  threshold  of  0.2  can  be  ad¬ 
justed  if  necessary  on  the  worksheet  named  “RedBombs.”  (We  do  not  include  the 
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smart  bomb  threshold  on  the  ExRatios  input  worksheet  because  we  recommend 
that  the  default  value  of  0.2  be  used.) 

The  remaining  “F  v.  Bluex”  and  “B  v.  Bluex”  columns  are  filled  with  loss  ratio 
values  for  Usage  values  of  1,  kill-rate  ratios  for  Usage  values  of  0,  and  the  ratio  of 
the  mean  time  for  Red  to  escape  to  the  mean  time  for  Blue  to  kill.  Those  values 
are  generally  determined  from  experience  or  simulation  and  are  beyond  the  scope 
of  this  users  guide.  It  is  expected  that  users  will  have  access  to  relevant  loss  ratio 
data.  Once  the  loss  ratios  are  completed,  the  large  button  at  the  bottom  of  the 
screen  labeled  “Get  Kill  Rate  Ratios”  needs  to  be  clicked.  That  button  will  run  a 
macro  that  uses  the  loss  ratios  to  compute  kill-rate  ratios.  The  kill-rate  ratios  are 
automatically  entered  onto  the  worksheet  labeled  “Rdata.”  As  an  added  safeguard, 
every  time  SLAACM  is  run,  the  kill-rate  ratio  macro  is  run  to  ensure  that  the  most 
recent  values  are  entered  onto  Rdata. 

SLAACM  will  not  run  if  the  ExRatios  worksheet  contains  any  blank  values.  If  a 
loss  ratio  is  left  blank  on  ExRatios,  the  error  message  shown  in  Figure  6-6  will 
appear. 


Figure  6-6.  Classic  SLAACM:  Loss  Ratio  Error 


Microsoft  L'.c 


Invalid  Loss  Ratio 


Then  the  user  will  be  taken  to  the  ExRatios  worksheet.  Any  blank  values  will  be 
highlighted  as  shown  in  Figure  6-7. 


Figure  6-7.  Classic  SLAACM:  Blank  Loss  Ratio  FI igh lighted 
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The  final  user  inputs  are  on  the  worksheet  named  “CM,”  which  stands  for  cruise 
missiles.  On  this  worksheet,  the  user  can  specify  the  number  of  cruise  missiles  a 
Red  aircraft  could  carry.  Figure  6-8  shows  that  Column  C  is  used  to  input  the 
cruise  missile  capacity  of  each  aircraft. 

Figure  6-8.  Classic  SLAACM:  CM  Worksheet 


No  tradeoff  between  cruise  missiles  and  bomb  payload  is  required  because  the 
cruise  missile  computation  is  a  secondary  analysis  that  is  meant  to  provide  “what- 
if’  information  to  the  user.  In  this  computation,  the  number  of  Red  aircraft  or  each 
type  that  get  past  Blue  defenses  each  day  is  multiplied  by  the  number  of  cruise 
missiles  specified  on  the  CM  worksheet.  This  analysis  merely  gives  an  idea  of  the 
potential  number  of  cruise  missiles  that  Red  could  launch  if  its  aircraft  were 
armed  to  do  so. 

Running  SLAACM 


Once  all  of  the  inputs  have  been  entered,  the  user  needs  to  return  to  the  BlueOOB 
worksheet.  To  launch  the  analysis,  the  user  clicks  the  “Run  Campaign”  button. 

If  the  heuristic  is  being  used,  no  other  action  is  required  on  the  part  of  the  user. 

If  LINGO  is  being  used,  then  a  series  of  message  boxes  will  appear  to  let  the  user 
know  what  day  is  being  considered  and  which  optimization  is  taking  place.  The 
user  simply  needs  to  click  “OK”  to  proceed.  The  optimization  in  LINGO  usually 
takes  only  a  couple  of  seconds,  and  the  screen  often  clears  automatically.  If  it 
does  not  clear  by  itself,  the  user  may  need  to  click  “Close,”  provided  LINGO  has 
solved  for  the  globally  optimal  solution  (noted  by  “Global  Opt”  in  the  “State”). 
This  case  is  shown  in  Figure  6-9. 
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Figure  6-9.  Classic  SLAACM:  LINGO  Solution  for  Red  Optimization 


If  the  “State”  shows  “Feasible”  as  its  status,  then  the  user  should  wait  until  the 
global  optimal  solution  is  reached.  If  the  solver  is  running  for  a  long  time,  and  the 
difference  between  the  “Best  Obj”  and  “Obj  Bound”  is  small  (e.g.,  less  than  0.5 
percent),  then  the  user  may  click  “Interrupt  Solver”;  the  best  feasible  solution  at 
that  point  will  be  entered  into  SLAACM.  This  event  occurs  rarely  in  SLAACM, 
since  LINGO  is  usually  able  to  solve  these  integer  programming  problems  very 
quickly. 

Figure  6-10  shows  LINGO  before  it  has  reached  optimality;  it  has  reached  a  fea¬ 
sible — but  not  necessarily  optimal — solution.  In  this  example,  the  difference 
between  the  best  objective  and  the  objective  bound  is  less  than  1/1,712  =  0.05 
percent,  so  one  could  certainly  interrupt  the  solver  at  this  point. 

Figure  6-10.  Classic  SLAACM:  Feasible  LINGO  Solution 
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Outputs 


Once  the  analysis  has  been  completed,  the  user  has  several  output  results  to  re¬ 
view.  SLAACM  produces  daily  records  of  Red  and  Blue  orders  of  battle,  of  Red 
and  Blue  losses,  of  their  variances  due  to  engagements,  of  tons  of  smart  and  not- 
smart  bombs  delivered,  and  of  cruise  missiles  delivered. 

Six  worksheets  report  SLAACM ’s  output  results: 

♦  BlueOOB 

♦  RedOOB 

♦  RedBombs 

♦  Losses 

♦  MBCharts 

♦  CM. 

We  discuss  the  outputs — charts  and  worksheets  tables — below.  By  reviewing  the 
charts  and  tables,  the  user  can  gain  an  understanding  of  specific  campaign  out¬ 
comes. 


Charts 


The  main  results  are  presented  graphically.  From  the  BlueOOB  worksheet,  the 
user  can  click  the  “View  Charts”  button,  or  simply  proceed  to  the  worksheet  tab 
labeled  “MB  Charts.”  The  MB  Charts  worksheet  contains  eight  charts. 

Figure  6-1  la  shows  the  quad  chart  “dashboard”  of  the  four  main  SLAACM  re¬ 
sults. 
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Figure  6-1  la.  Classic  SLAACM:  Output  Charts 


The  “Red  Forces”  chart  in  the  upper  left  corner  of  Figure  6-1  la  shows  the  quan¬ 
tity  of  Red  aircraft  on  each  day  of  the  campaign.  In  the  example,  most  Red  air¬ 
craft  are  annihilated  by  Day  4.  The  aircraft  that  have  some  quantity  remaining,  but 
which  reach  steady-state  quantities,  are  most  likely  not  being  sent.  For  example, 
Red  F 1  and  Red  F2  are  weak  fighters,  so  they  are  never  sent;  instead,  Red  is 
choosing  to  escort  its  bombers  with  more  capable  aircraft.  A  quick  check  of  the 
“Tons  of  Bombs  Dropped”  chart  in  the  lower  right  comer  shows  that  it  looks  like 
the  bomb  levels  were  decreasing  by  Day  5,  indicating  that  fewer  Red  bombers  are 
getting  through  Blue  defenses  beyond  that  day;  however,  after  Day  5  we  see  a 
steady  increase  in  bombs  dropped.  We  can  see  that  no  additional  Blue  fighters  are 
sent  after  Day  5  and  a  low  level  of  Red  bombers  continue  to  reach  their  target. 
That  is  because  we  had  set  up  Blue  to  be  highly  averse  to  losing  aircraft.  If  the 
Blue  value  of  Blue  Loss  parameter  is  set  to  10  (half  of  the  Red  value),  then  the 
Blue  aircraft  are  sent  and  the  campaign  is  over  by  Day  5. 

The  “Blue  Forces”  chart  in  the  upper  right  corner  of  Figure  6-1  la  shows  the  daily 
quantities  of  Blue  aircraft.  We  see  that  Blue  sustains  some  losses,  although  not 
nearly  as  many  as  Red;  therefore,  we  can  assume  that  Blue’s  forces  are  superior  to 
Red’s.  In  addition,  in  the  first  few  days,  we  see  increases  in  some  quantities  of 
Blue  aircraft.  That  is  due  to  Blue  receiving  reinforcements.  (See  the  BlueSupply 
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worksheet  in  Figure  6-1.)  Blue  LADs  are  broken  out  separately  from  the  Blue 
forces  and  are  shown  in  the  chart  in  the  lower  left  comer  of  Figure  6-1  la. 

The  Red  Forces,  Blue  Forces,  and  LAD  charts  show  the  quantity  of  each  aircraft 
type  that  is  available  on  each  day  of  the  campaign.  The  charts  do  not  indicate  the 
number  of  aircraft  sent  into  battle. 

Additional  charts  on  the  worksheet  are  shown  in  Figure  6-1  lb.  Those  charts  show 
number  of  kills  by  Blue  and  Red  aircraft.  Colors  and  dimensions  are  used  to  dis¬ 
play  comparisons  of  kills  to  show  how  specific  aircraft  perform. 

Figure  6-1 1b.  Classic  SLAACM:  Output  Charts  (continued) 


Figure  6-12  shows  a  graph  of  the  potential  number  of  cruise  missile  launches.  As 
discussed  in  the  previous  section,  this  result  assumes  that  the  Red  bombers  that 
got  through  Blue  defenses  were  armed  with  cruise  missiles  in  lieu  of  bombs  and 
were  able  to  launch  their  cruise  missiles. 
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Figure  6-12.  Classic  SLAACM:  Cruise  Missile  Launches 
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Worksheet  Tables 

Because  SLAACM  is  housed  in  an  Excel  workbook,  the  user  can  browse  the  data 
used  to  build  the  charts.  All  of  the  charts  in  Figure  6-11  are  built  from  data  in 
worksheet  tables: 

♦  The  Red  Forces  chart  uses  the  data  on  the  “RedOOB”  worksheet. 

♦  The  Blue  Forces  and  LAD  charts  use  data  from  the  “BlueOOB”  work¬ 
sheet.  (As  we  discussed  earlier,  the  BlueOOB  worksheet  is  both  an  input 
and  an  output.  After  running  SLAACM,  available  aircraft  quantities  are 
loaded  into  the  BlueOOB  worksheet  for  each  day  of  the  campaign.) 

♦  The  Tons  of  Bombs  chart  uses  data  from  the  “RedBombs”  worksheet. 

♦  The  Cruise  Missile  chart  is  built  from  data  on  the  “CM”  worksheet,  where 
the  chart  resides. 

It  can  be  helpful  to  look  at  aircraft  losses  to  get  a  better  idea  of  how  the  campaign 
progressed.  SLAACM  has  a  worksheet  named  “Losses”  that  shows  the  loss  data. 
Figure  6-13  shows  the  Losses  worksheet.  Blue  and  Red  losses  are  shown  by  air¬ 
craft  type.  In  addition,  the  worksheet  shows  the  standard  deviation  of  the  losses 
for  Blue  aircraft. 
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Figure  6-13.  Classic  SLA  ACM:  Losses 


Figure  6-14  shows  the  pounds  and  tons  of  bombs  dropped  by  Red  aircraft  each 
day  of  the  campaign.  In  addition,  the  bomb  types  are  broken  out  by  smart  and  not- 
smart  types.  These  are  the  data  used  to  build  the  “Tons  of  Bombs  Dropped”  chart 
in  the  lower  right  comer  of  Figure  6-1  la.  Upon  inspection  of  the  data  in 
Figure  6-14,  we  see  that  the  campaign  actually  ends  on  Day  5,  not  Day  4  as  it  had 
appeared  on  the  chart;  one  can  see  that  2  tons  of  bombs  were  dropped  on  Day  5. 

The  RedBombs  worksheet  also  contains  a  table  that  shows  bomber  payloads, 
cepr,  and  bomber  values.  The  values  for  bomber  effectiveness  used  in  the  optimi¬ 
zations  are  calculated  based  on  payload  tonnage  and  accuracy  as  in  previous  ver¬ 
sions  of  SLAACM.  Now,  however,  those  calculated  values — along  with  the 
fighter  values — are  displayed  on  the  RedBombs  worksheet  to  support  “eyeball” 
analysis  by  the  user.  The  RedBombs  spreadsheet  with  this  display  is  shown  in 
Figure  6-14. 
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Figure  6-14.  Classic  SLAACM:  Red  Bombs 


User  Note 

The  output  worksheets  in  SLAACM  are  dynamic  and  are  overwritten 
during  each  model  run.  Users  who  want  to  save  run  data  should  copy 
and  paste  the  values  into  an  external  file. 


Analyses  of  Sample  Runs 


In  this  section,  we  conduct  a  comparative  analysis  of  the  heuristic  and  LINGO 
using  two  SLAACM  runs  to  illustrate  some  of  the  subtle  differences  that  occur 
under  reasonably  similar  scenarios. 

First,  we  run  SLAACM  using  the  heuristic,  with  the  inputs  shown  in  Figures  6-1 
through  6-5;  we  then  repeat  the  run  using  LINGO.  As  shown  in  Table  6-1,  the  re¬ 
sults  are  identical;  the  heuristic  performs  the  same  as  the  LINGO  optimization. 
This  case  is  not  uncommon,  and  it  gives  us  confidence  that  the  heuristic  can  be  a 
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good  approximation  for  the  optimal  solution  (the  heuristic  and  the  exact  integer 
programming  formulations  are  covered  in  detail  in  Chapter  2). 


Table  6-1.  Classic  SLAACM:  Results  of  Initial  SLAACM  Run 


Item 

Heuristic 

LINGO 

Tons  dropped 

Smart  bombs 

J  1,241 

|  1,241 

Dumb  bombs 

174 

174 

Aircraft  losses 

Blue 

19 

19 

Red 

397 

397 

To  illustrate  a  situation  in  which  the  heuristic  and  LINGO  solutions  differ,  con¬ 
sider  the  same  scenario,  except  that  all  Blue  aircraft  are  smart.  This  change  re¬ 
quires  checking  the  boxes  on  the  BlueOOB  worksheet  (Figure  6-3).  The  results 
are  shown  in  Table  6-2.  For  both  optimizers,  the  tonnage  of  smart  bombs  dropped 
is  reduced  by  almost  50  percent,  demonstrating  the  impact  of  smart  Blues  prefer¬ 
entially  engaging  the  highest-value  Red  packages.  Table  6-2  also  shows  signifi¬ 
cant  differences  in  the  heuristic  results  and  LINGO’s  exact  IP  results. 

Table  6-2.  Classic  SLAACM:  Results  of  Setting  All  Blues  to  Be  Smart 


Item 

Heuristic 

LINGO 

Tons  dropped 

Smart  bombs 

640 

602 

Dumb  bombs 

289 

289 

Aircraft  losses 

Blue 

19 

17 

Red 

402 

404 

LINGO’s  optimal  solution  resulted  in  38  fewer  tons  of  smart  bombs  being 
dropped,  2  fewer  Blue  aircraft  losses,  and  2  additional  Red  aircraft  losses.  In 
short,  the  IP  solution  using  LINGO  resulted  in  Blue  reducing  the  amount  of  dam¬ 
age  while  sustaining  fewer  losses.  This  example  illustrates  what  we  have  seen  in 
other  results,  as  well  as  what  one  might  expect.  The  heuristic  and  the  exact  solu¬ 
tion  to  the  IP  problem  in  LINGO  generally  produce  qualitatively  the  same  results. 
However,  if  one  is  looking  for  very  precise  results,  the  heuristic  may  fall  short  of 
true  optimality  by  several  percentage  points. 
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Summary 


Classic  SLAACM  is  a  fast,  flexible,  robust  model  that  contains  the  key  input  pa¬ 
rameters  necessary  to  define  the  characteristics  of  a  realistic  air-to-air  campaign. 
Model  results  in  both  tabular  and  graphical  fonnat  clearly  display  the  impacts  of 
parameter  choices  and  provide  insight  into  campaign  scenarios.  The  Excel  work¬ 
book  implementation  allows  results  to  be  easily  copied  and  pasted  to  other  appli¬ 
cations  for  reporting  and  presentation. 

The  examples  show  differences  in  results  using  the  heuristic  optimizer  and  solv¬ 
ing  the  exact  IP  formulation  with  LINGO.  In  the  many  practical  cases  we  have 
analyzed,  it  has  been  rare  for  the  heuristic  results  to  differ  from  the  LINGO  results 
by  more  than  10  percent,  but  careful  analyses  may  require  the  use  of  exact  solu¬ 
tion  methods,  especially  for  actual  campaign  planning. 
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Chapter  7 

Users  Guide  to  Attack  SLAACM 


Attack  SLAACM  is  housed  in  a  Microsoft  Excel  workbook  with  extensive  Visual 
Basic  for  Applications  code.  The  current  version  of  Attack  SLAACM  was  de¬ 
signed  using  object-oriented  programming  principles.  In  addition,  it  has  an  option 
to  conduct  integer  programming  by  linking  to  an  external  optimization  program 
called  LINGO.1  Using  Excel  to  house  SLAACM  has  several  advantages.  Most 
analysts  are  familiar  with  the  Excel  environment  and  will  be  able  to  navigate  it 
easily.  Other  than  acquiring  the  SLAACM  workbook,  users  do  not  need  to  install 
additional  software  unless  they  want  to  run  integer  programming  optimizations  in 
LINGO.  Finally,  users  can  easily  copy  and  paste  output  values  and  charts  into 
their  preferred  presentation  formats  for  displaying  results  to  others. 

In  previous  SLAACM  documentation,  friendly  forces  are  identified  as  Blue  and 
Green,  and  enemy  forces  are  Red.  With  the  ability  of  friendly  forces  to  defend  or 
attack,  we  have  changed  the  nomenclature  in  this  report  to  “Attack”  and 
“Defense.”  The  previous  version  has  been  called  “Red  SLAACM,”  “Air  Defense 
SLAACM,”  and  “Classic  SLAACM,”  and  the  new  SLAACM  version  incorporat¬ 
ing  attack  has  been  called  “Blue  SLAACM”  and  “Attack  SLAACM.”  For  consis¬ 
tency,  in  this  chapter  we  use  the  terms  “Classic  SLAACM”  and  “Attack 
SLAACM”  to  differentiate  between  the  two  versions.  Despite  the  names,  the  new 
version  includes  both  defensive  and  offensive  scenarios  and  is  a  superset  of  the 
older  version. 

This  chapter  describes  how  to  operate  Attack  SLAACM.  It  begins  with  a  descrip¬ 
tion  of  the  model  inputs.  It  then  describes  how  to  run  the  model  once  all  of  the 
input  parameters  are  set  up.  Next,  the  chapter  describes  the  outputs  reporting  At¬ 
tack  SLAACM ’s  results.  The  last  section  addresses  how  the  outputs  are  analyzed 
to  determine  the  effects  of  the  various  engagement  scenarios. 


Inputs 


Attack  SLAACM ’s  input  parameters  are  defined  in  six  worksheets: 

♦  BlueSupply,  which  is  used  to  specify  the  quantities  of  Defense  aircraft 
available  on  each  day  of  the  campaign. 

♦  RedSupply,  which  is  used  to  specify  the  quantities  of  Attack  aircraft  avail¬ 
able  on  each  day  of  the  campaign. 


1  For  more  information  on  LINGO,  see  http://www.lindo.com/products/lingo/. 
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♦  BlueOOB,  which  is  used  to  define  the  order  of  battle  and  other  pertinent 
parameters.  BlueOOB  is  the  main  sheet  from  which  SLAACM  runs  are 
conducted.  (SLAACM  also  has  a  RedOOB  worksheet,  but  it  is  purely  an 
output  worksheet  and  does  not  receive  any  user  inputs.) 

♦  ExRatios,  which  is  used  to  specify  engagement  parameters  such  as  loss  ra¬ 
tios  and  breakpoints. 

♦  SiteData,  which  is  used  to  define  SEAD  vs.  SAM  and  SAM  vs.  bomber 
parameters. 

♦  CM,  which  is  used  to  specify  the  number  of  cruise  missiles  that  a  Red  air¬ 
craft  could  carry. 

We  discuss  each  sheet  below  and  describe  the  proper  way  to  complete  each  one. 

BlueSupply  and  RedSupply  Worksheets 

For  each  force,  Blue  (Defense)  and  Red  (Attack),  the  user  will  need  to  set  up  the 
initial  quantities  of  aircraft  available  and  the  reinforcements  available  on  each  day 
of  the  campaign.  The  user  will  need  to  type  those  quantities  into  the  appropriate 
worksheets,  called  BlueSupply  and  RedSupply. 

Figure  7-1  is  an  example  of  the  BlueSupply  worksheet. 

Figure  7-1.  Attack  SLAACM:  BlueSupply  Worksheet 
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As  shown  in  the  figure,  the  worksheet  has  a  row  for  each  Blue  (Defense)  aircraft 
type  and  a  column  for  each  day  of  the  campaign.  Initial  aircraft  quantities  are 
entered  into  the  column  labeled  0.  Reinforcements  may  be  provided  for  subse- 
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quent  days  of  the  campaign.  In  the  figure,  we  see  that  reinforcements  are  entered 
for  Day  2  and  Day  3. 2 


The  RedSupply  worksheet  is  designed  in  the  same  way  as  the  BlueSupply  sheet. 
Figure  7-2  is  an  example.  In  the  figure,  we  see  that  only  initial  aircraft  quantities 
are  provided  and  no  reinforcements  are  specified. 

Figure  7-2.  Attack  SLAACM:  RedSupply  Worksheet 
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The  BlueSupply  and  RedSupply  worksheets  can  accommodate  up  to  20  aircraft 
types  each.  If  new  rows  are  needed  to  accommodate  additional  aircraft  types,  they 
can  be  added  using  the  “Add  a  New  Row”  button. 


Caution! 

Any  time  a  row  is  added  or  deleted,  it  is  imperative  that  the  buttons  on  the  Blue¬ 
Supply  and  RedSupply  worksheets  are  used  to  add  or  delete  the  rows. 

Do  not  simply  insert  or  delete  rows  and  columns. 


Several  of  SLAACM’s  worksheets  reference  the  aircraft  types.  The  “Add  a  New 
Row”  and  “Delete  a  Row”  buttons  must  be  used  to  keep  the  worksheets  synchro¬ 
nized.  Using  typical  worksheet  approaches  to  inserting  and  deleting  rows  will 
likely  result  in  model  errors  or  misinterpretation  of  results. 


2  Note  that  a  day  is  used  in  the  model  as  a  convenient  designator  to  represent  an  individual 
attack/defense  sortie.  In  an  actual  campaign,  there  may  be  more  than  one  sortie  on  a  particular 
calendar  day,  or,  conversely,  sorties  may  be  distributed  among  several  calendar  days. 
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BlueOOB  Worksheet 


The  worksheet  named  BlueOOB  (which  stands  for  Blue  order  of  battle)  is  the 
main  worksheet  for  input  parameters  for  the  Blue  (Defense)  aircraft.  It  is  also  the 
worksheet  used  to  launch  campaign  analyses.  BlueOOB  also  serves  as  an  output 
worksheet,  as  does  the  RedOOB  worksheet;  the  output  portion  of  BlueOOB  will 
be  discussed  in  the  output  section  of  this  chapter.  Figure  7-3  shows  the  BlueOOB 
worksheet. 

Figure  7-3.  Attack  SLAACM:  BlueOOB  Worksheet 
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Caution! 

Rows  or  columns  should  never  be  deleted  from  this  worksheet! 


The  aircraft  listed  in  the  BlueOOB  worksheet  are  linked  to  those  shown  in  the 
BlueSupply  worksheet  (Figure  7-1).  The  Day  0  quantities  are  also  linked  to  those 
in  the  BlueSupply  worksheet.  The  user  who  attempts  to  click  the  aircraft  names  or 
Day  0  quantities  on  the  BlueOOB  sheet  will  receive  the  warning  shown  in 
Figure  7-4. 

Figure  7-4.  Warning:  Do  Not  Alter  Aircraft  of  BlueOOB  Worksheet 
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The  first  four  columns  on  the  BlueOOB  worksheet  define  aircraft-specific  input 
parameters  for  Blue: 

♦  “Smart?  ”  column.  The  user  may  use  this  column  to  identify  specific  Blue 
(defense)  aircraft  as  smart.  Smart  Blue  aircraft  can  determine  the  composi¬ 
tion  of  Red  attack  packages  before  intercepting  them  and  can  coordinate 
their  defensive  response  among  themselves  to  provide  optimal  defenses. 
This  coordination  and  optimization  is  discussed  in  Chapter  5.  Multiple 
Blue  aircraft  types  may  be  designated  as  smart.  In  Figure  7-3,  the  first  air¬ 
craft  is  selected  to  be  smart. 

♦  “LAD?  ”  column.  The  user  may  use  this  column  to  identify  certain  Blue 
(defense)  aircraft  as  local  air  defenders.  SLAACM  treats  LADs  differently 
from  other  Blue  types  when  considering  Blue  dispatch  options.  In  general, 
a  non-LAD  Blue  has  a  penalty  for  its  own  losses  that  is  significantly 
greater  than  a  LAD  Blue,  which  presumably  is  less  loss  averse  due  to  de¬ 
fending  its  own  territory.  Some  non-LAD  Blues  may  have  a  negative  pay¬ 
off  function  for  engaging  certain  Red  packages  dispatched.  Blue  will  not 
dispatch  that  type.  LADs,  however,  are  dispatched  in  all  cases.  Multiple 
Blue  aircraft  types  may  be  designated  as  LADs.  In  Figure  7-3,  the  third 
aircraft  is  designated  as  such. 

♦  “CAP  Factor”  column.  This  factor  addresses  the  situations  in  which  Blue 
(defense)  fighters  are  stationed  at  remote  airfields  and  must  divide  their 
forces  to  maintain  continuous  combat  air  patrol  coverage  over  the  battle 
space.  The  user  may  select  a  CAP  factor  for  any  Blue  aircraft  type,  so  that 
only  a  specified  fraction  of  the  inventory  of  that  Blue  type  is  available. 
CAP  factor  values  are  chosen  from  a  drop-down  box  for  each  aircraft  type. 
Values  range  from  1  to  10.  A  CAP  factor  of  1  indicates  that  no  CAP  limi¬ 
tation  is  placed  on  that  aircraft  type.  In  Figure  7-3,  all  aircraft  have  CAP 
factors  of  1 .  However,  if  an  aircraft  had  been  assigned  a  CAP  factor  of  4, 
that  aircraft  would  have  only  n/4  available  aircraft  on  the  first  day.  If  there 
are  not  enough  aircraft  to  fulfill  the  combat  air  patrol  at  all  times,  those 
flights  are  not  dispatched. 

♦  “ Index  of  Blue  a/c  for  Red’s  planning ”  column.  Only  one  aircraft  type  can 
be  selected  for  this  parameter.  In  Figure  7-3,  the  first  aircraft  type  is  des¬ 
ignated  as  Red’s  “planning  aircraft.”  This  aircraft  type  is  used  by  Red  to 
compute  its  expected  payoff  function  and  optimize  its  dispatch  choices. 
Experience  has  shown  that  Red  optimization  is  relatively  insensitive  to  the 
selection  of  the  Blue  type,  so  nonnal  practice  is  to  select  the  most  numer¬ 
ous  Blue  type  for  Red  planning. 

In  addition  to  the  aircraft-specific  parameters,  several  other  input  parameters  must 
be  selected  on  the  BlueOOB  worksheet. 
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At  the  top  of  the  worksheet,  the  number  of  days  in  the  campaign  is  specified.  In 
Figure  7-3,  the  number  of  campaign  days  is  4.  The  longer  the  campaign,  the 
longer  it  will  take  SLAACM  to  run,  since  each  campaign  day  requires  the  full 
breadth  of  SLAACM  computations.  Initially,  the  user  may  need  to  start  out  with  a 
longer  campaign  and  then  decrease  its  length  as  a  steady  state  in  aircraft  losses  is 
reached. 

Beneath  the  number  of  campaign  days  is  a  parameter  labeled  “Battle  Management 
r|.”  This  parameter  is  the  battle  management  efficiency  factor,  q  is  a  value  greater 
than  0  and  less  than  or  equal  to  1 .  Simply  stated,  q  less  than  1  indicates  the  per¬ 
centage  of  time  that  Blue  defender  will  be  unable  to  find  a  Red  attack  package.  In 
Figure  7-3,  the  efficiency  factor  is  shown  as  0.9,  which  means  that  10  percent  of 
the  not-smart  Blue  defenders  will  chase  “ghost”  attackers.  No  losses  will  result, 
but  in  those  cases,  Blue  will  not  be  confronting  an  actual  attacker,  so  its  efficiency 
is  reduced.  Since  smart  Blues  know  Red’s  attack  package  compositions  in  ad¬ 
vance,  we  assume  they  have  perfect  management. 

Next  to  the  campaign  days  and  the  efficiency  factor  is  a  drop-down  box  labeled 
“LINGO?”  If  “No”  is  selected,  as  shown  in  Figure  7-3,  SLAACM  will  use  its  in¬ 
ternal  heuristic  optimization  to  calculate  engagement  results.  If  “Yes”  is  selected, 
LINGO  will  be  used  to  perform  integer  optimization  (a  licensed  copy  of  LINGO 
software  is  required  to  be  installed  on  the  PC  for  this  feature  to  work).  Both  the 
heuristic  and  the  exact  integer  programming  solution  are  described  in  detail  in 
Chapter  5.  If  LINGO  is  used,  there  are  two  LINGO  input  worksheets:  “RedOpt” 
and  “BlueOpt.”  However,  those  sheets  are  already  complete  and  should  not  be 
altered  by  the  user. 

Below  the  parameters  just  described  are  four  value  boxes: 

♦  Red  value  of  Blue  loss 

♦  Blue  value  of  Blue  loss 

♦  Red  value  of  Red  loss 

♦  Blue  value  of  Red  loss. 

These  parameters  establish  the  value  of  killing  an  opponent,  as  well  as  the  penalty 
a  side  takes  for  losing  one  of  its  own.  These  values,  which  are  applied  to  fighter 
aircraft  only,  are  used  in  the  optimization  objective  function.  The  Red  value  of 
Blue  loss  represents  the  benefit  to  Red  of  destroying  a  Blue  aircraft,  while  Blue 
value  of  Blue  loss  is  the  penalty  to  Blue  for  losing  one  of  its  own  aircraft.  Values 
can  be  between  zero  and  100. 

In  the  case  shown  in  Figure  7-3,  we  see  that  the  value  of  a  Blue  loss  to  Red  is  20, 
while  the  penalty  to  Blue  for  losing  its  own  aircraft  is  only  1.  Therefore,  Blue’s 
aversion  to  a  loss  is  one -twentieth  the  value  of  Red’s  benefit  for  obtaining  the 
Blue  loss.  Conversely,  we  see  in  this  example  that  Red  is  not  penalized  for  a  Red 
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fighter  loss  (the  value  is  zero)  and  Blue  achieves  a  benefit  for  a  Red  fighter  loss 
equal  to  Red’s  benefit  for  killing  a  Blue  aircraft  (the  value  is  20). 

Further  to  the  right  are  two  drop-down  boxes  for  Attack  package  and  Defender 
package  designations.  New  engagement  scenarios  have  been  added  to  Attack 
SLAACM  to  represent  potential  U.S.  attack  package  options  using  both  nonnal 
and  low-observable  aircraft .  Table  7-1  lists  the  attack  package  options. 

Table  7-1  Attack  Package  Options 


Drop-down  value 

Attack  package 

2 

Basic  hostile  12-ship  (4  LE,  4  CE,  4  B)  attack  package 

3 

4  escort,  2  SEAD,  2  bomber  attack  package 

5 

4  escort,  2  bomber  attack  package 

6 

2  escort,  2  SEAD,  1  bomber  attack  package 

7 

2  escort,  1  bomber  attack  package 

The  defender  package  selections  are  enumerated  in  Table  7-2. 

Table  7-2.  Defender  Package  Options 


Drop-down  value 

Defender  package 

2 

2  defenders 

4 

4  defenders  (basic  friendly  defense  package) 

8 

8  defenders  (basic  hostile  defense  package) 

On  the  far  right  hand  side  of  Figure  7-3  are  two  yes/now  drop-down  boxes.  The 
first  is  labeled  “Victorious  F/B  complete  mission?”  The  second  is  labeled  “Escap¬ 
ing  bombers  complete  mission?”  If  yes  is  chosen,  then  fighter-bombers  and 
bombers,  respectively,  that  escape  attack  go  on  to  complete  their  mission  and 
deliver  bombs  to  their  target.  If  no  is  chosen,  then  escaping  fighter-bombers  and 
bombers  do  not  complete  their  bombing  missions. 

The  final  input  selection  on  BlueOOB  is  a  drop-down  box  labeled  “Alerts.”  The 
default  setting  for  this  feature  is  “Off.”  However,  if  “On”  is  selected,  message 
boxes  will  become  visible  as  the  analysis  runs. 
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ExRatios  Worksheet 


The  ExRatios  worksheet  contains  a  row  for  each  Red  attack  aircraft  that  has  been 
designated  on  the  RedSupply  worksheet.  The  worksheet  has  several  parameter  col¬ 
umns,  discussed  below,  and  a  column  for  each  Blue  defender  aircraft  labeled  “F  v. 
Bluex”  where  “F”  indicates  that  Red  is  a  fighter  and  “x”  corresponds  to  the  label 
provided  on  BlueSupply.  Columns  are  repeated  for  each  Blue  aircraft,  except  they 
are  labeled  “B  v.  Blue  x”  where  “B”  indicates  that  Red  is  a  bomber  and  “x”  is  as 
just  described.  The  default  setting  in  SLAACM  assumes  that  loss  ratios  are  for 
engagements  between  4  Blues  and  8  Reds,  with  the  Reds  and  Blues  fighting  to 
annihilation. 

Those  assumptions  can  be  changed  on  this  sheet  under  “Red  Start,”  “Blue  Start,” 
“Red  Quit,”  and  “Blue  Quit.”  Red  Start  and  Blue  Start  correspond  to  the  engage¬ 
ment  type,  e.g.,  8  and  4,  respectively.  Red  Quit  and  Blue  Quit  denote  the  breaking 
points;  for  example,  0  and  0  means  that  both  sides  will  fight  to  annihilation.  The 
user  will  want  to  specify  values  that  correspond  to  the  conditions  under  which  the 
loss  ratios  were  obtained. 

To  keep  the  rows  and  columns  synchronized  with  the  BlueSupply  and  RedSupply 
worksheets,  it  is  imperative  that  no  rows  or  columns  on  ExRatios  be  inserted  or 
deleted. 


Caution! 

As  stated  earlier,  use  the  add  and  delete  buttons  on  the  BlueSupply  and  RedSupply 
worksheets  to  add  or  delete  rows. 


Figure  7-5  is  a  screenshot  of  the  ExRatios  worksheet. 
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Figure  7-5.  Attack  SLAACM:  ExRatios  Worksheet 
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Only  red-colored  values  are  inputs  on  this  sheet 

Fighter-Bombers  (Usage  =  0)  engage  Blue  A/C,  even  when  going  as  Bombers  Values  posted  on  this  sheet  represent  loss-ratios  for  those 
engagements 

Bombers-Only  (Usage  =  -1 )  do  not  engage  Blue  A/C,  and  attempt  to  escape  Values  posted  on  this  sheet  represent  the  ratio  of  (mean  time  for 
Red  to  escape  /  (mean  time  for  Blue  to  kill). 


*  Usage  =  -1  =>  a/c  is  bomber  only;  0  =>  bomber  OR  fighter;  1  =>  fighter  only;  2  =  SEAD.  NB:  a/c  usable  as  bombers  must  h. 


Info  for  Calculating  Kill-Rate  Ratios 
Blue  Start  Red  Start  Blue  Quit  Red  Quit 


8 


Get  Kill  Rate 
Ratios 


Sel  SEAD  vs. 
SAM  values 


h  <  ►  m\  BlueSupply  /  BlueOOB  /  RedSupply  /  RedOOB  /  BlueKIs  /  RedKIs  \  ExRatios/  SteData  ,| 
Ready 


m 


The  first  column  on  the  ExRatios  worksheet  is  labeled  “Usage.”  For  each  Red  air¬ 
craft  type,  a  value  of  0,  1,  -1,  or  2  needs  to  be  entered.  A  value  of  0  indicates  a 
fighter-bomber.  A  value  of  1  indicates  a  fighter.  When  bombers  are  designated  by 
the  Usage  =  -1  code,  they  are  assumed  to  try  to  escape  rather  than  engage  the 
Blue  fighters.  Because  the  bombers  are  escaping  and  not  engaging  in  combat,  the 
values  for  fighter  vs.  bomber  encounters  in  the  ExRatios  spreadsheet  no  longer 
represent  loss  ratios.  Rather,  they  represent  the  ratio  of  the  mean  time  for  Red  to 
escape  to  the  mean  time  for  Blue  to  kill.  Bombers  with  the  Usage  =  0  code  are 
assumed  to  be  fighter  bombers  and  will  engage;  their  encounter  values  in  the 
ExRatios  spreadsheet  represent  loss  ratios.  A  value  of  2  indicates  a  SEAD  aircraft. 
(A  note  reminding  the  user  of  this  difference  is  on  the  ExRatio  spreadsheet  and  is 
shown  in  Figure  7-5.) 


The  second  column  indicates  the  name  of  the  Red  aircraft  type;  it  is  linked  to  the 
RedSupply  worksheet.  The  third  column,  “Payload,  lb,”  is  used  to  enter  the  num¬ 
ber  of  pounds  of  bombs  that  the  associated  Red  aircraft  can  carry. 
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The  fourth  column  is  labeled  “cepr,”  which  stands  for  the  circular  error  probable 
ratio  (CEP/CEPref).  The  cepr  is  a  value  greater  than  0  and  less  than  or  equal  to  1, 
where  1  indicates  that  the  bomber  cannot  deliver  smart  bombs  (guided  munitions) 
and  0.2  (or  less)  indicates  that  it  can.  The  smart  bomb  threshold  of  0.2  can  be  ad¬ 
justed  if  necessary  on  the  worksheet  named  “RedBombs.”  (We  do  not  include  the 
smart  bomb  threshold  on  the  ExRatios  input  worksheet  because  we  recommend 
that  the  default  value  of  0.2  be  used.) 

The  remaining  “F  v.  Bluex”  and  “B  v.  Bluex”  columns  are  filled  with  loss  ratio 
values  for  Usage  values  of  1,  kill-rate  ratios  for  Usage  values  of  0,  and  the  ratio  of 
the  mean  time  for  Red  to  escape  to  the  mean  time  for  Blue  to  kill.  Those  values 
are  generally  determined  from  experience  or  simulation,  and  are  beyond  the  scope 
of  this  users  guide.  It  is  expected  that  users  will  have  access  to  relevant  loss  ratio 
data.  Once  the  loss  ratios  are  completed,  the  large  button  at  the  bottom  of  the 
screen  labeled  “Get  Kill  Rate  Ratios”  needs  to  be  clicked.  That  button  will  run  a 
macro  that  uses  the  loss  ratios  to  compute  kill-rate  ratios.  The  kill-rate  ratios  are 
automatically  entered  onto  the  worksheet  labeled  “Rdata.”  As  an  added  safeguard, 
every  time  SLAACM  is  run,  the  kill-rate  ratio  macro  is  run  to  ensure  that  the  most 
recent  values  are  entered  onto  Rdata. 


SLAACM  will  not  run  if  the  ExRatios  worksheet  contains  any  blank  values.  If  a  loss 
ratio  is  left  blank  on  ExRatios,  the  error  message  shown  in  Figure  7-6  will  appear. 


Figure  7-6.  Attack  SLAACM:  Loss  Ratio  Error 


Then  the  user  will  be  taken  to  the  ExRatios  worksheet  and  any  blank  values  will 
be  highlighted  as  shown  in  Figure  7-7. 


Figure  7-7.  Attack  SLAACM:  Blank  Loss  Ratio  FI igh lighted 
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SiteData  Worksheet 

The  SiteData  worksheet  is  particular  to  the  Attack  SLAACM  model.  Using  this 
worksheet,  users  define  SEAD  and  SAM  parameters.  Figure  7-8  is  a  screenshot  of 
the  SiteData  worksheet. 

Figure  7-8.  Attack  SLAACM:  SiteData  Worksheet 


Aircraft  designated  as  SEAD  by  using  the  value  of  2  on  ExRatios  worksheet  can 
have  their  kill-rate  ratios  set  on  the  SiteData  worksheet.  In  Figure  7-8,  the  fourth 
listing  is  a  SEAD  aircraft.  On  the  SiteData  worksheet,  a  KRR  of  1 .2  has  been  en¬ 
tered,  indicating  that  this  SEAD  aircraft  has  a  KRR  of  1.2  against  a  SAM  site.  For 
the  SAM  site,  a  loss  can  be  either  destruction  or  shutdown  for  the  duration  of  the 
engagement.  For  the  SEAD,  a  loss  means  loss  of  the  aircraft.  For  aircraft  that  are 
not  designated  as  SEAD,  KRR  values  on  this  sheet  are  ignored,  so  by  convention 
we  set  them  to  zero. 

In  addition,  on  the  SiteData  worksheet,  the  single-shot  kill  probability  for  a  SAM 
missile  against  a  bomber  can  be  set  in  the  box  labeled  “Site  sspk.” 

Finally,  the  number  of  missiles  that  a  SAM  site  fires  at  the  set  of  bombers  reach¬ 
ing  their  target  is  set  using  the  box  labeled  “#  of  shots.” 

SEAD  aircraft  engage  the  SAM  site  (the  SAM  has  unlimited  missiles  for  this 
duel).  The  SAM  site  does  not  engage  the  escort  fighters,  and  the  defending  fight¬ 
ers  do  not  engage  the  SEAD  aircraft.  If  the  SEAD  aircraft  are  destroyed  by  the 
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SAM,  the  SAM  site  (with  up  to  three  missiles)  engages  the  bombers.  In  scenarios 
with  no  SEAD  aircraft,  the  SAM  site  always  engages  the  bombers  with  three  mis¬ 
siles.  The  limitation  to  three  missiles  accounts  for  the  limited  time  interval  on 
which  the  bombers  are  on  their  bomb  run,  rather  than  for  any  limit  to  the  number 
of  missiles  at  the  SAM  site. 

CM  Worksheet 

The  final  user  inputs  are  on  the  worksheet  named  “CM,”  which  stands  for  cruise 
missiles.  On  this  worksheet,  the  user  can  specify  the  number  of  cruise  missiles  a 
Red  aircraft  could  carry.  Figure  7-9  shows  that  Column  C  is  used  to  input  the 
cruise  missile  capacity  of  each  aircraft. 

Figure  7-9.  Attack  SLAACM:  CM  Worksheet 
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No  tradeoff  between  cruise  missiles  and  bomb  payload  is  required  because  the 
cruise  missile  computation  is  a  secondary  analysis  that  is  meant  to  provide  “what- 
if’  information  to  the  user.  In  this  computation,  the  number  of  Red  aircraft  of  each 
type  that  get  past  Blue  defenses  each  day  is  multiplied  by  the  number  of  cruise 
missiles  specified  on  the  CM  worksheet.  This  analysis  merely  gives  an  idea  of  the 
potential  number  of  cruise  missiles  that  Red  could  launch  if  its  aircraft  were 
armed  to  do  so. 
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Running  SLAACM 


Once  all  of  the  inputs  have  been  entered,  the  user  needs  to  return  to  the  BlueOOB 
worksheet.  To  launch  the  analysis,  the  user  clicks  the  “Run  Campaign”  button. 

If  the  heuristic  is  being  used,  no  other  action  is  required  on  the  part  of  the  user. 

If  LINGO  is  being  used,  then  a  series  of  message  boxes  will  appear  to  let  the  user 
know  what  day  is  being  considered  and  which  optimization  is  taking  place.  The 
user  simply  needs  to  click  “OK”  to  proceed.  The  optimization  in  LINGO  usually 
takes  only  a  couple  of  seconds,  and  the  screen  often  clears  automatically.  If  it 
does  not  clear  by  itself,  the  user  may  need  to  click  “Close,”  provided  LINGO  has 
solved  for  the  globally  optimal  solution  (noted  by  “Global  Opt”  in  the  “State”). 
This  case  is  shown  in  Figure  7-10. 

Figure  7-10.  Attack  SLAACM:  LINGO  Solution  for  Red  Optimization 


If  the  “State”  shows  “Feasible”  as  its  status,  then  the  user  should  wait  until  the 
global  optimal  is  reached.  If  the  solver  is  running  for  a  long  time,  and  the  differ¬ 
ence  between  the  “Best  Obj”  and  “Obj  Bound”  is  small  (e.g.,  less  than  0.5  per¬ 
cent),  then  the  user  may  click  “Interrupt  Solver”;  the  best  feasible  solution  at  that 
point  will  be  entered  into  SLAACM.  This  event  occurs  rarely  in  SLAACM,  since 
LINGO  is  usually  able  to  solve  these  integer  programming  problems  very  quickly. 

Figure  7-11  shows  LINGO  before  it  has  reached  optimality;  it  has  reached  a  fea¬ 
sible — but  not  necessarily  optimal — solution.  In  this  example,  the  difference  be¬ 
tween  the  best  objective  and  the  objective  bound  is  less  than  1/1,712  =  0.05 
percent,  so  one  could  certainly  interrupt  the  solver  at  this  point. 
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Figure  7-11.  Attack  SLAACM:  Feasible  LINGO  Solution 


Outputs 


Once  the  analysis  has  been  completed,  the  user  has  several  output  results  to  re¬ 
view.  Attack  SLAACM  produces  daily  records  of  Red  and  Blue  orders  of  battle, 
of  Red  and  Blue  losses,  of  their  variances  due  to  engagements,  of  tons  of  smart 
and  not-smart  bombs  delivered,  and  of  cruise  missiles  delivered. 

Seven  worksheets  report  Attack  SLAACM ’s  output  results: 

♦  BlueOOB 

♦  RedOOB 

♦  RedBombs 

♦  Losses 

♦  MBCharts 

♦  CM 

♦  sam. 

We  discuss  the  outputs — charts  and  worksheet  tables — below.  By  reviewing  the 
charts  and  tables,  the  user  can  gain  an  understanding  of  specific  campaign  out¬ 
comes. 
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Charts 


The  main  results  are  presented  graphically.  From  the  BlueOOB  worksheet,  the 
user  can  click  the  “View  Charts”  button,  or  simply  proceed  to  the  worksheet  tab 
labeled  “MB  Charts.”  There  are  eight  charts  on  the  MB  Charts  worksheet. 

Figure  7-12a  shows  the  “dashboard”  of  the  four  main  SLAACM  results. 

The  “Red  Forces”  chart  in  the  upper  left  corner  of  Figure  7-12a  shows  the  quan¬ 
tity  of  Red  Attack  aircraft  on  each  day  of  the  campaign.  In  the  example,  most  Red 
aircraft  are  annihilated  by  Day  5.  A  quick  check  of  the  “Tons  of  Bombs  Dropped” 
chart  in  the  lower  right  comer  shows  that  it  looks  like  the  bomb  levels  were  de¬ 
creasing  by  Day  6,  indicating  that  fewer  Red  bombers  were  getting  through  Blue 
defenses  beyond  that  day.  However,  the  “Blue  Forces”  chart  shows  a  steady  de¬ 
cline  in  Blue  Defenders,  which  means  that  some  level  of  attackers  continue  to  en¬ 
gage  with  them. 

Blue  LADs  are  broken  out  separately  from  the  Blue  Forces  and  are  shown  in  the 
chart  in  the  lower  left  comer  of  Figure  7-12a. 

The  Red  Forces,  Blue  Forces,  and  LAD  charts  show  the  quantity  of  each  aircraft 
type  that  is  available  on  each  day  of  the  campaign.  The  charts  do  not  indicate  the 
number  of  aircraft  sent  into  battle. 

The  attack  and  defender  packages  that  were  selected  on  BlueOOB  are  displayed 
on  the  right  side  of  Figure  7-12a. 
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Figure  7-1 2a.  Attack  SLAACM:  Output  Charts 
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Additional  charts  on  the  worksheet  are  shown  in  Figure  7- 12b.  Those  charts  show 
number  of  kills  by  Blue  and  Red  aircraft.  Colors  and  dimensions  are  used  to  dis¬ 
play  comparisons  of  kills  to  show  how  specific  aircraft  perform. 
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Figure  7-1 2b.  Attack  SLAACM:  Output  Charts  (continued) 
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Figure  7-13  shows  a  graph  of  the  potential  number  of  cruise  missile  launches.  As 
discussed  in  the  previous  section,  this  result  assumes  that  the  Red  bombers  that 
got  through  Blue  defenses  were  armed  with  cruise  missiles  in  lieu  of  bombs  and 
were  able  to  launch  their  cruise  missiles. 

Figure  7-13.  Attack  SLAACM:  Cruise  Missile  Launches 


The  “sam”  worksheet  shows  results  of  SAM  and  SEAD  activities.  Figure  7-14  is  a 
screenshot  of  this  worksheet.  The  “sam”  worksheet  shows,  by  day,  the  number  of 
SAMs  killed,  SEADs  killed,  and  bombers  killed  by  SAMs.  In  the  above  example, 
we  see  that  about  175  SAMs,  86  SEADs,  and  19  bombers  are  killed  in  this  sce¬ 
nario. 
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Figure  7-14.  Attack  SLAACM:  SAM,  SEAD,  and  Bomber  Results 


Worksheet  Tables 

Because  SLAACM  is  housed  in  an  Excel  workbook,  the  user  can  browse  the  data 
used  to  build  the  charts.  All  of  the  charts  in  Figure  7-11  are  built  from  data  in 
worksheet  tables: 

♦  The  Red  Forces  chart  uses  the  data  on  the  “RedOOB”  worksheet. 

♦  The  Blue  Forces  and  LAD  charts  use  data  from  the  “BlueOOB”  work¬ 
sheet.  (As  we  discussed  earlier,  the  BlueOOB  worksheet  is  both  an  input 
and  an  output.  After  running  SLAACM,  available  aircraft  quantities  are 
loaded  into  BlueOOB  worksheet  for  each  day  of  the  campaign.) 

♦  The  Tons  of  Bombs  chart  uses  data  from  the  “RedBombs”  worksheet. 

♦  The  Cruise  Missile  chart  is  built  from  data  on  the  “CM”  worksheet,  where 
the  chart  resides. 

♦  The  SAMs,  SEADs,  and  bombers  chart  is  built  from  data  on  the  “sam” 
worksheet. 

It  can  be  helpful  to  look  at  aircraft  losses  to  get  a  better  idea  of  how  the  campaign 
progressed.  SLAACM  has  a  worksheet  named  “Losses”  that  shows  the  loss  data. 
Figure  7-15  shows  the  Losses  worksheet.  Blue  and  Red  losses  are  shown  by  air¬ 
craft  type.  In  addition,  the  worksheet  shows  the  standard  deviation  of  the  losses 
for  Blue  aircraft. 
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Figure  7-15.  Attack  SLAACM:  Losses 


Figure  7-16  shows  the  pounds  and  tons  of  bombs  dropped  by  Red  aircraft  each 
day  of  the  campaign.  In  addition,  the  bomb  types  are  broken  out  by  smart  and  not- 
smart  types.  These  are  the  data  used  to  build  the  “Tons  of  Bombs  Dropped”  chart 
in  the  lower  right  comer  of  Figure  7-12a. 
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The  RedBombs  worksheet  also  contains  a  table  showing  bomber  payloads,  cepr, 
and  bomber  values.  The  values  for  bomber  effectiveness  used  in  the  optimizations 
are  calculated  based  on  payload  tonnage  and  accuracy,  as  in  previous  versions  of 
SLAACM.  Now,  however,  those  calculated  values — along  with  the  fighter  val¬ 
ues — are  displayed  on  the  RedBombs  worksheet  to  support  “eyeball”  analysis  by 
the  user.  The  RedBombs  spreadsheet  with  this  display  is  shown  in  Figure  7-16. 

Figure  7-16.  Attack  SLAACM:  Red  Bombs 
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The  output  worksheets  in  SLAACM  are  dynamic  and  are  overwritten  during  each 
model  run.  Users  who  want  to  save  run  data  should  copy  and  paste  the  values  into 

an  external  file. 
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Summary 


Attack  SLAACM  is  a  fast,  flexible,  robust  model  that  contains  the  key  input  pa¬ 
rameters  necessary  to  define  the  characteristics  of  a  realistic  air-to-air  campaign. 
Model  results  in  both  tabular  and  graphical  fonnat  clearly  display  the  impacts  of 
parameter  choices  and  provide  insight  into  campaign  scenarios.  The  Excel  work¬ 
book  implementation  allows  results  to  be  easily  copied  and  pasted  to  other  appli¬ 
cations  for  reporting  and  presentation. 
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Appendix  A 

Direct  Computation  of  Long-Time  Limiting 
Probabilities  of  Boundary  States 


The  structure  of  the  evolution  equations  for  probability  distributions  of  the 
engagements  we  consider  makes  it  possible  to  find  the  long-tenn  limits  of  the 
absorbing  boundary  states  directly,  in  finitely  many  steps,  without  formally 
solving  those  differential  equations.  In  some  interesting  cases,  the  finite  steps  may 
be  carried  out  by  straightforward  iteration.  This  appendix  explains  these  facts  and 
gives  some  examples. 

The  evolution  equations  for  the  state  probabilities  in  our  engagement  models 
decompose  naturally  into  a  set  of  equations  for  the  transient  states  and  a  set  for 
the  absorbing  boundary  states.  Equations  2-5  and  2-6  are  examples  of  this  for  the 
basic  M  vs.  N  probabilistic  engagement  model.  The  advanced  engagement  models 
of  Appendix  D  also  have  this  property,  except  for  the  low-observable  model, 
which  is  treated  completely  by  the  method  given  in  Appendix  D. 

If  we  make  the  probabilities  of  the  transient  states  the  components  of  a  vector  x, 
and  the  components  of  the  absorbing  boundary  states  into  the  components  of  a 
vector  y,  the  evolution  of  the  state  probabilities  can  be  described  by 


x  =  Ax  ,  x(0)  =  x0 
y  =  Bx 


[Eq.  A-l] 


To  illustrate  these  ideas  for  an  example  that,  unlike  the  one  of  Chapter  2,  has 
recurrent  states,  let  us  consider  treating  a  1  vs.  1  engagement  with  a  two-phase  kill 
model  with  break-lock.  We  use  the  state  description  (m,  i;  n,  j)  where  m  denotes 
the  number  of  Blue  aircraft,  i  the  number  of  Blues  tracking  opponents,  n  the 
number  of  Red  aircraft,  and  j  the  number  of  Reds  tracking  opponents.  Figure  A-l 
is  a  diagram  of  the  states  in  this  engagement  and  their  transitions. 
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Figure  A-1.  Diagram  of  1  vs.  1  Engagement 
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The  evolution  equations  are 

Pioio  =  — (kb1  +  kr1)P1010  +bxP1011  +rxP1110 
Pmo  =  — (bb2  +krj  +rx)P1110  +kbjP1010  +bxPnn 
Pion  =  — ' (^2  +kbj  +bx)P1011  +  krjP1010  +rxPnll 
Pun  =  — (kb2  +  kr2  +rx  +bx)Pmi  +kr1P1110  +kbjP1011 
Piooo  =  kb2(P1110  +Pnn) 

Pooio  =  kr2(Pion  +P 

mi) 


[Eq.  A-2] 


The  first  four  of  these  describe  the  transient  states,  and  the  last  two  describe  the 
absorbing  boundary  states.  Defining  xi  as  Pioio,  x 2  as  Pmo,  X3  as  P1011,  and  x4  as 
Pi  111,  and  defining  yi  as  Piooo,  yi  as  Pooio,  we  see  that,  in  this  case,  the  matrices  A 
and  B  of  Equation  A-1  are 


(kb!  +  kr, ) 

rx 

bx 

0  ) 

kbj 

-0*2  +krl  +rx) 

0 

bx 

[Eq.  A-3] 
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"0  kb 2  0  kb2" 

v0  0  kr2  kr2 


[Eq.  A-4] 
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Now,  the  solutions  of  systems  of  ordinary  differential  equations  of  the  fonn 
x  =  Ax  for  constant  coefficient  matrices  A  can  always  be  expressed  as  finite 
linear  combinations  of  generalized  exponential  functions  of  the  time,  that  is, 

functions  of  the  fonn  tJ  eXt ,  where  the  j  are  positive  integers  and  the  A  are  the 
eigenvalues  of  the  matrix  A.1  If  the  transients  are,  in  fact,  transient,  then  the  real 
parts  of  the  k  are  all  strictly  less  than  0.  It  follows  that  the  x; ,  defined  by 


=  {  xi(t)dt  [Eq.A-5] 

o 

exist.  Then,  on  integrating  the  differential  equations  of  Equation  A-l  from  0  to  oo, 
using  the  initial  condition  on  x  and  y,  and  remembering  that  the  x;(t)  tend  to  zero 
as  t  — »  oo,  we  find 


-  x0  =  Ax 

Ylim  =  B* 


[Eq.  A-6] 


where  ynm  denotes  the  limit  of  y  as  t  — »  oo.  It  follows  that 


Ylim 


-BA~ 


x 


o  • 


[Eq.  A-7] 


Equation  A-7  expresses  the  long-time  limiting  values  of  the  absorbing  boundary 
state  probabilities  as  the  result  of  finite  operations,  that  is,  matrix  inversion  and 
multiplication. 

Let  us  continue  to  illustrate  these  concepts  with  the  1  vs.  1  example  considered 
above.  For  a  numerical  example,  we  take  kri  =  1,  kr2  =  2,  rx  =  1,  kbi  =  3,  kbi  =  4, 
and  bx  =  2.  Then  the  matrices  A  and  B  take  the  values 


and 


'-4  1  2 

3-6  0 

1  0 
0  1  3 


0  ' 
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B  = 


4  0 
0  2 


4^ 

2. 


[Eq.  A-8] 


[Eq.  A-9] 


It  is  as  well  to  check  that  the  eigenvalues  of  A  are  distinct  and  negative.  Direct 
calculation  shows  that  this  is  in  fact  the  case:  the  eigenvalues  are  approximately 


1  W.E.  Boyce  and  R.C.  DiPrima,  Elementary’  Differential  Equations ,  Third  Edition  (New 
York:  Wiley,  1976),  Sections  7-6  through  7-9. 


A-3 


-2.44,  -6.44,  -6.56,  and  -10.56.  Assured  by  this  that  the  Xj  exist,  we  go  on  to 
find 


Ylim 


-BA-1x 


0  “ 


"0.827^ 

v0.173y 


[Eq.  A- 10] 


That  is,  the  probability  that  the  Blue  aircraft  defeats  the  Red  one  is  roughly 
83  percent. 

The  eigenvalues  of  the  plant  matrix  A  have  another  use:  they  tell  one  whether  or 
not  it  makes  sense  to  look  for  the  long-time  limiting  values  of  the  absorbing 
boundary- state  probabilities.  Those  probabilities  are  of  interest  if,  but  only  if,  an 
actual  engagement  can  continue  long  enough  for  the  system  to  be  in  boundary 
states  with  a  probability  near  one.  The  time  for  this  to  happen  is  the  time  for 
which  the  transient  states’  probabilities  are  all  much  less  than  one. 

When  the  eigenvalues  of  A  are  distinct  (and,  of  course,  have  negative  real  parts) 
so  that  the  Xi(t)  are  linear  combinations  of  the  functions  cxp(Xjt),  that  time  will  be 
no  greater  than  the  time  at  which  exp(A*t)  becomes  much  less  than  one,  where  A* 
is  the  eigenvalue  closest  to  zero.  The  time  to  make  that  happen  is  approximately 
2/A.*.  If  the  rate  parameters  kbi,  ...  of  our  example  are  in  inverse  minutes,  the 
mean  time  for  the  Blue  aircraft  to  acquire  a  target  is  20  seconds,  the  mean  time  for 
them  to  launch  a  successful  missile  is  15  seconds,  and  the  mean  time  for  them  to 
break  a  Red  lock  is  30  seconds,  while  the  mean  time  for  the  Reds  to  make  lock  is 
60  seconds,  the  mean  time  for  them  to  launch  a  successful  missile  is  30  seconds, 
and  the  mean  time  for  them  to  break  lock  is  60  seconds.  The  negative  eigenvalue 
of  the  plant  matrix  A  with  smallest  magnitude  is  -2.44  inverse  minutes,  which 
implies  that  the  time  for  the  transient  phase  of  the  1  vs.  1  engagement  to  be  over  is 
roughly  50  seconds.  Very  likely,  the  combatants  will  have  enough  fuel  to  fight 
that  long,  and  so  the  long-time  limit  is  meaningful  in  our  example. 

Although  we  have  found  numerical  evaluation  of  Equation  A-7  quite  helpful  in 
generating  insight  with  examples  of  modest  dimension,  determining  long-time 
limiting  probabilities  in  this  way  may  not  be  practical  for  systems  of  large 
dimension.  In  an  important  class  of  engagement  models,  however,  the  difficult 
task  of  computing  A'1  xo  can  be  done  with  a  straightforward  (although  possibly 
lengthy)  iterative  scheme.  That  class  is  those  engagement  models  whose  diagrams 
are  acyclic,  like  the  example  of  Chapter  2.  For  these  models,  the  plant  matrix  A  is 
lower  triangular,  and  this  provides  the  iterative  scheme.  In  these  cases,  one  has 

aiAi  =  “1 

a2iXi  +  a22x2  =  0  [Eq.  A-ll] 
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so  that 


Xj  =  -l/an 

x2  =  -Xj  a 21  /a22  [Eq.  A-12] 


With  the  x  determined,  the  long-time  limiting  values  of  the  absorbing  boundary 
states  follow  from  Equation  A-7. 

Determining  the  set  of  states  actually  occupied  in  a  given  engagement  can  be 
somewhat  tedious.  For  example,  finding  the  long-time  limit  of  a  missile-tracking 
engagement  model  for  4  vs.  4,  when  each  aircraft  has  6  missiles,  involves  more 
than  42,000  states.  Obviously,  it  is  not  practical  to  evaluate  these  cases  by  hand. 
In  these  cases,  we  generally  use  C++  code  to  find  both  the  set  of  states  occupied 
and  the  solution  (Equation  A-12). 
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Appendix  B 

Analyzing  Large  State-Space  Engagement 
Models  with  ASSIST  and  STEM 


Evolution  equations  like  Equations  2-5  and  2-6  of  Chapter  2  are  simple  in  principle. 
They  are  systems  of  linear  ordinary  differential  equations  with  constant  coefficients. 
The  solutions  of  the  initial  value  problems  formed  by  such  equations  and  specifica¬ 
tions  of  the  starting  probabilities  of  the  various  states  can  be  written  symboli¬ 
cally — and  perhaps  treated  practically — with  matrix  exponentiation,  and  they 
generally  offer  no  difficulty  to  numerical  solution  by  difference  schemes,  but  they 
can  be  of  large  dimension. 

A  basic  4  vs.  4  engagement  has  24  states  and  32  transitions;  adding  6  missiles  to 
the  Blue  aircraft  produces  1,632  states  and  9,087  transitions;  adding  6  missiles  to 
both  Red  and  Blue  produces  108,534  states  and  2,416,252  transitions.  Also,  it  is 
somewhat  awkward  to  detennine  the  long-time  limits  of  the  absorbing  boundary 
states  by  integrating  the  differential  equations  numerically.  This  appendix  dis¬ 
cusses  NASA-developed  tools  that  we  have  found  helpful  in  dealing  with  both 
these  issues. 

The  decision  to  apply  the  NASA  tools  to  the  problem  of  fighter  combat  was  ser¬ 
endipitous.  One  of  us  was  engaged  in  Markov  analysis  of  fighter  combat,  and  an¬ 
other  of  us  was  using  the  NASA  tools  for  safety-related  reliability  analysis.  The 
need  for  a  tool  to  efficiently  develop  relatively  complex  fighter  combat  models 
surfaced  during  informal  discussions.  Using  the  NASA  tools,  we  can  develop  and 
run  complex  models  for  Blue  versus  Red  combat,  including  tracking  of  Blue  mis¬ 
sile  use,  in  relatively  short  order.  These  results  provide  insight  into  the  combat 
problem  and  support  development  of  the  SLAACM  algorithms. 

NASA  has  three  tools  that  can  be  applied  to  fighter  combat.  Two  are  Markov  compu¬ 
tational  analysis  programs.  The  third  is  a  sophisticated  utility  program  that  generates 
the  inputs  (models)  for  the  analysis  programs.  The  tools  were  developed  by  NASA  to 
estimate  failure  probabilities  in  highly  reliable,  reconfigurable  avionics  and  space 
electronics. 

The  two  analysis  programs  are  Scaled  Taylor  Exponential  Matrix  (STEM)1  and 
Semi-Markov  Unreliability  Range  Estimator  (SURE).2  The  former  is  a  pure 
Markov  analysis  tool  in  which  all  failure  rates  are  constant  with  state  probabilities 


1  NASA,  The  PAWS  and  STEM  Reliability’  Analysis  Programs,  Technical  Memorandum 
100572,  R.  Butler  and  P.  Stevenson,  March  1988. 

2  NASA,  SURE  Reliability’  Analysis,  Program  and  Mathematics,  Technical  Paper  2764, 
R.  Butler  and  A.  White,  1998. 
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having  generalized  exponential  fonn.  The  latter  is  a  semi-Markov  model  that  al¬ 
lows  use  of  non-exponential  reconfiguration  probabilities.  STEM  and  SURE  use 
identical  input  files,  but  some  commands  are  processed  only  by  SURE.  For  rea¬ 
sons  discussed  below,  STEM  is  the  tool  used  for  fighter  combat  calculations.  The 
utility  program  ASSIST  (Abstract  Semi-Markov  Specification  Interface  to  the 
SURE  Tool)3  generates  the  STEM  and  SURE  input  files,  or  “models.”  ASSIST 
allows  the  straightforward  generation  of  extremely  complex  Markov  models.4 

Reliability  analyses  are  typically  conducted  for  short  model  times  compared  to  the 
failure  rates  of  the  components  studied,  such  as  a  10-hour  flight  and  a  failure  rate 
of  0.0001  failure  per  hour.  For  air  combat  models,  we  are  interested  in  running  the 
engagement  to  completion  so  we  will  run  model  time  units  on  the  order  of  the  re¬ 
ciprocal  of  the  lowest  kill  rate,  e.g.,  1,000  units  for  a  kill  rate  of  0.001.  STEM  has 
no  problem  running  such  models,  but  SURE’s  mathematical  algorithms  that  han¬ 
dle  non-Markov  recovery  rates  are  not  designed  to  run  to  completion.  Conse¬ 
quently,  the  remainder  of  the  discussion  will  focus  on  ASSIST  and  STEM. 

All  three  programs  were  developed  for  UNIX  platforms  and  have  been  converted 
to  Windows.  STEM  is  also  available  on  LINUX.  The  programs  are  available  at  no 
charge  from  NASA  Langley  Research  Center.  Documentation  includes  an 
ASSIST  users  guide,  a  report  on  modeling  techniques,  and  reports  on  SURE  and 
STEM  mathematics  and  performance.  Figure  B-l  shows  the  basic  relationships  of 
the  tools. 


Figure  B-1.  Basic  Relationships  of  the  Tools 


HOL:  Higher  Order  Language 

ASSIST:  Abstract  Semi-Markov  Specification  Interface  to  the  SURE  Tool  -  Input  model  generator 

Model:  Designation  for  the  input  files  for  the  Markov  analysis  programs,  containing  variables,  states,  transitions,  and  rates 
STEM:  Scaled  Taylor  Exponential  Matrix  -  Markov  analysis  program 

SURE:  Semi-Markov  Unreliability  Range  Estimator  -  Markov  and  Semi-Markov  analysis  program 


3  NASA  Langley  Research  Center.  ASSIST  User  Manual,  S.  Johnson  and  D.  Boerschlein, 
September  1993. 

4  NASA,  Techniques  for  Modeling  the  Reliability  of  Fault-Tolerant  Systems  with  the  Markov 
State-Space  Approach,  Reference  Publication  1348,  R.  Butler  and  S.  Johnson,  September  1995. 
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Fighter  Combat  Analysis 


In  this  section,  we  use  case  examples  of  increasing  complexity  to  describe  the  ca¬ 
pabilities  and  limitations  of  the  tools  for  analysis  of  fighter  combat.  All  the  analy¬ 
ses  use  the  common  structure  of  an  M  vs.  N,  Blue  vs.  Red  engagement  with 
constant  kill  rates. 

Example  1 — Simple  2  vs.  2  Engagement 

The  2  vs.  2  engagement  demonstrates  the  basic  structure  of  the  models  without 
generating  extensive  output.  The  Blue  fighter  kill  rate  is  0.1  per  unit  time,  and  the 
Red  fighter  kill  rate  is  0.01  per  unit  time.  Listing  1  is  the  ASSIST  code  that  gener¬ 
ates  the  2  vs.  2  model. 

Listing  1.  ASSIST  Code  for  Generating  the  2  vs.  2  Model 

(*  ASSIST  model  for  Fighter  Combat  *) 

LIST  =  3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

(*  equipment  list  *) 

nblue  =2;  (*  blue  fighters  *) 

nred  =2;  (*  red  fighters  *) 

(*  kill  probabilities  *) 
k_rate_blue  =  0.1; 
k_rate_red  =  0.01; 

SPACE  =  (blue:  0. nblue,  red:  O.nred); 

START  =  (nblue,  nred) ; 

DEATH IF  (blue  =  0)  AND  (red>0); 

DEATH IF  (red  =  0)  AND  (blue>0); 

(*  Blue  kills  *) 

IF  (blue  >  0)and  (red>0)  TRANTO  red  =  red-1  BY  blue*k_rate_blue; 

(*  Red  kills  *) 

IF  (red  >  0  )  and  (blue>0)  TRANTO  blue  =  blue-1  BY  red*k_rate_red; 


ASSIST  uses  a  higher  order  definition  language  that  allows  algebraic  manipulation 
of  variables  and  compact  description  of  Markov  state  transfers.5  Several  features 
are  noteworthy  in  the  listing  above.  The  listing  starts  with  the  editing  commands 
LIST  and  ONEDEATH  that  control  SURE/STEM  output.  Next,  the  number  of  air¬ 
craft  and  their  kill  rates  are  input  parameters  defined  as  constant  types.  The  SPACE 
statement  defines  the  range  of  the  Markov  state  space,  i.e.,  in  this  case,  the  first 
state  can  vary  from  0  to  nblue  where  nblue  equals  2.  The  START  statement  identi¬ 
fies  the  initial  state  populations.  The  DEATHIF  statements  define  the  absorbing 
states  in  the  model;  in  this  case,  the  two  DEATHIF  conditions  represent  Blue  and 
Red  wins.  If  the  ONDEATH  OFF  command  is  enabled,  STEM  and  SURE  list  all 


5  ASSIST  inputs  are  fully  defined  in  the  citations  of  Notes  1-3. 
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the  individual  absorbing  states;  otherwise,  they  combine  the  results  for  absorbing 
states  into  the  appropriate  Death  states.6 

We  now  come  to  the  powerful  TRANTO  statements,  which  define  the  conditions, 
nature,  and  rate  of  Markov  state  transfers.  In  this  example,  the  transfers  occur 
whenever  both  Blue  and  Red  fighters  are  available;  kills  are  made  one  at  a  time, 
and  the  kill  rates  are  proportional  to  the  number  of  killers.  The  two  TRANTO 
statements  above  are  adequate  to  generate  the  models  for  any  combination  of  Red 
and  Blue  aircraft. 

The  state  diagram  for  the  2  vs.  2  engagement,  Figure  B-2,  is  relatively  simple.  It 
shows  that  there  are  8  total  states,  of  which  4  are  absorbing  (Death)  states  and  8 
are  transitions.7 


Figure  B-2.  2  vs.  2  Engagement  Diagram 


Designates  an  intermediate  state  with  X  Blue  fighters  and  Y  Red  fighters 


Designates  an  absorbing  (Death)  state  with  X  Blue  fighters  no  Red  fighters  or  visa  versa 


XB  and 


Designate  Blue  and  Red  kill  rates 


6  SURE  mathematics  requires  the  existence  of  absorbing  states,  while  STEM  can  handle  mod¬ 
els  that  include  recovery  from  all  states. 

7  There  is  no  (0,  0)  absorbing  state  because  we  chose  not  to  include  the  case  in  which  the  last 
two  opponents  simultaneously  shoot  each  other  down.  Such  a  state  could  be  modeled  with  an  ad¬ 
ditional  TRANTO  statement,  but  determination  of  the  controlling  kill  rate  would  take  some  care¬ 
ful  thought.  The  pure  statistical  rate  for  independent  events  suggests  that  it  would  be  equal  to 

K  RATE  BLUE  *  K  RATE  RED,  but  simultaneous  kills  can  happen  only  in  certain  conditions 
such  as  ramming,  head-to-head  gun  attacks,  and  head-to-head  missile  attacks.  These  are  only  a 
subset  of  the  configurations  that  make  up  the  basic  Blue  and  Red  kill  rates.  We  assessed  the  likely 
rate  to  be  sufficiently  low  to  justify  ignoring  the  (0,  0)  state  until  better  data  become  available. 
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Running  ASSIST  with  the  listing  above  generates  the  STEM/SURE  input  file  or 
“model”  shown  in  Listing  2.  Note  the  8  states,  the  4  absorbing  states  (0,  N  or  M, 
0),  and  the  8  transitions.  The  file  extension  is  “.mod.” 


Listing  2.  2  vs.  2  Engagement  Model 


LIST  =  3; 

TIME  =  100; 

NBLUE  =  2; 

NRED  =  2; 

MIN_B  =  0; 

MIN_R  =  0; 
K_RATE_BLUE  =  0.1; 
K  RATE  RED  =  0.01; 


1  (* 

2,2 

*) 

2  (* 

t — 1 

CXI 

*) 

=  2*K 

RATE_ 

_BLUE ; 

1  (* 

2,2 

*) 

3  ( * 

1,2 

*) 

=  2*K 

RATE 

_RED; 

2  (* 

2,1 

*) 

4  (* 

2,0 

*) 

=  2*K 

RATE 

BLUE ; 

2  (* 

2,1 

*) 

5  (* 

t — 1 

t — 1 

*) 

* 
t — 1 

II 

RATE 

_RED; 

3  ( * 

1,2 

*) 

5  (* 

t — 1 

t — 1 

*) 

=  1*K 

RATE 

BLUE ; 

3(* 

1,2 

*) 

6  (* 

0,2 

*) 

=  2*K 

RATE 

_RED; 

5  (* 

1,1 

*) 

7  (* 

o 

5 - 1 

*) 

=  1*K 

RATE 

BLUE ; 

5  (* 

1,1 

*) 

8  (* 

5 - 1 

o 

*) 

II 

I—4 

* 

RATE_ 

_RED; 

(*  NUMBER  OF  STATES  IN  MODEL  =  8  *) 

(*  NUMBER  OF  TRANSITIONS  IN  MODEL  =  8  *) 


The  model  file  contains  the  input  definitions  and  transition  descriptions  needed  by 
SURE  and  STEM.  Transitions  are  defined  by  source  state  number,  destination 
state  number,  and  rate  of  transfer  between  the  source  and  destination.  The  state 
descriptions,  e.g.,  (*2,2*),  are  included  as  comments;  “(*”  and  “*)”  are  comment 
delimiters.  Note  that  the  absorbing  states,  such  as  state  4,  (*2,0*),  appear  only  in 
the  destination  column.  The  DEATHIF  definitions  and  ONEDEATH  command 
are  not  printed.  Listing  3  is  the  STEM  output  file  for  the  model  above,  run  for  10 
time  units. 


Listing  3.  2  vs.  2  Engagement  STEM  Output 


Model  =  C:\Markovl\fighters\Report2x2.mod 
-  RUN  #1 

D-STATE  PROBABILITY  ACCURACY 


4  8.65800856430E-0001 

6  1. 5151310349 5E-0002 

7  1.08204982318E-0001 

8  1.08204982318E-0002 


TOTAL  9.99977647330E-0001  11  DIGITS 
STATE  PROBABILITY 


1  2.78946809287E-0010 

2  9. 5861846700 9E-0009 

3  1. 2287866813 0E-0006 

4  8. 6580085643 0E-0001 

5  2 . 11140182554E-0005 

6  1. 5151310349 5E-0002 

7  1.08204982318E-0001 

8  1.08204982318E-0002 
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Figure  B-3  shows  the  STEM  results  for  increasing  units  of  combat  time.  As  dis¬ 
cussed  in  Chapter  2,  we  use  the  reciprocal  of  the  lowest  rate  as  a  rule  of  thumb  for 
setting  the  model  run  time.  In  this  case,  the  lowest  rate  is  KRATERED  =  0.01 
and  its  reciprocal  is  100.  By  Time  =  70,  the  total  probability  of  being  in  an  ab¬ 
sorbing  state  is  0.9995,  and  at  Time  =  100,  it  is  0.99998,  indicating  that  our  rule 
of  thumb  is  good.  The  state  probabilities  at  Time  =100  for  the  individual  absorb¬ 
ing  states  are  (2B,  OR)  =  0.8658,  (IB,  OR)  =  0.1082,  (0B,  2R)  =  0.0152,  and  (0B, 
1R)  =  0.0108. 

Figure  B-3.  2  vs.  2  Engagement  Results  versus  Engagement  Time 


Elapsed  Time  Units 


Example  2 — 2  vs.  2  Engagement  with  Blue  Missiles 

This  example  adds  the  tracking  of  Blue  missiles  to  the  simple  2  vs.  2  engagement. 
We  assume  that  Blue  aircraft  each  carry  6  missiles,  each  having  a  known  single 
shot  probability  of  kill,  P^.  It  is  important  to  understand  in  this  discussion  that 
Blue  and  Red  kills  are  still  dependent  only  on  the  kill  rate  ratio,  and  that  missile 
Pk  is  used  only  to  calculate  missile  consumption.  For  simplicity  and  clarity,  we 
assume,  for  this  example,  that  Red  aircraft  have  unlimited  missiles. 


Our  basic  approach  is  to  expand  the  state  space  to  (bO,  bl,  b2,  b3,  b4,  b5,  b6,  r)  to 
include  separate  states  for  Blue  aircraft  having  0  to  6  missiles.  We  still  need  only 
one  state  vector  for  Red  aircraft.  At  the  start  of  an  engagement,  we  have  2  Blue 
aircraft  in  the  b6  6-missile  state  and  2  Red  aircraft  in  the  r  Red  state. 


We  still  transition  from  Red  states  in  only  two  ways:  Red  can  kill  a  Blue,  or  vice 
versa.  Now,  however,  we  can  transition  from  a  Blue  state  in  three  ways:  Blue  can 
kill  a  Red,  Blue  can  fire  and  miss,  or  Red  can  kill  a  Blue.  To  formulate  the  trans¬ 
fer  statements,  we  need  to  determine  the  missile  miss  rate.  We  know  the  Blue  kill 
rate,  X,  and  we  know  the  missile  single  shot  probability  of  kill,  P^.  We  want  to 
estimate  the  missile  usage,  including  misses.  We  assume,  for  now,  that  the  suc¬ 
cessful  missile  usage  is  one  missile  per  kill,  i.e.,  we  do  not  fire  salvos.  We  also 
assume  the  kill  rate  includes  both  the  missile  that  hit  and  any  that  missed. 
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Now,  we  define  /u  to  be  the  missile  firing  rate,  and  let  n  be  the  number  of  Blue 
aircraft  engaged.  The  successful  missile  rate  is  then:  n*Pk  *  |u,  and  the  unsuccess¬ 
ful  missile  rate  (miss  rate)  is  n*(l-Pk)*p. 

Based  on  our  assumptions,  the  successful  missile  rate  must  equal  the  kill  rate,  i.e., 

n*Pk*p  =  n*k  [Eq.  B-l] 

and,  therefore,  the  missile  firing  rate  is 

n  *  l i  =  n  *  — ,  and  p  =  —  .  [Eq.  B-2] 

Pk  Pk 

Substituting  for  p,  the  miss  rate,  n*(l-Pk)*|u,  becomes 


n*(l-Pk)*—  or  n*(1  Pk)*k.  [Eq.  B-31 

Pk  Pk 

Figure  B-4  is  the  state  transition  diagram  for  the  single  missile  per  kill  case. 

Figure  B-B-4.  State  Diagram  for  a  Single  Missile  per  Kill  Engagement 


Missile  Firing  Rate  =  p 

Missile  Kill  Rate  =  PkP 

Missile  Miss  Rate  =  (l-Pk)p 

Engagement  Kill  Rates  =  Xb,Xr 
Aircraft  in  Engagement  =  n  Blue,  m  Red 


We  can  generalize  the  derivation  above  in  a  straightforward  manner  for  cases  in 
which  multiple  missiles  are  fired  per  engagement  by  noting  that  the  failure  prob¬ 
ability  for  x  failures  is  (1  -  P£f  and  substituting  [1  -  (l-Pk)x]  for  Pk: 

♦  Define  //  to  be  the  missile  firing  rate. 

♦  Let  n  be  the  number  of  Blue  aircraft  engaged. 
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♦  Let  x  be  the  number  of  missiles  fired  per  kill. 

The  successful  missile  rate  is  n*[  1-  (1  -  Pk)A]  *  p,  and  the  unsuccessful  missile 
rate  (miss  rate)  is  n*(l  -  Pk)x*p. 

Based  on  our  assumptions,  the  successful  missile  rate  must  equal  the  kill  rate,  i.e., 

n*[l  -(1  -Pk)x]*p  =  n*L  [Eq.  B-4] 

and,  therefore,  the  missile  firing  rate  is 


n*p  =  n* 


X 


i ,  and  p  —  r 


X 


[i-(i-Pk)" 

The  miss  rate,  n*(l  -  Pk)x*p,  now  becomes 

a-pkr 


n*(l-Pk)? 


or  n ' 


1-(1-Pk)x  l-(l-Pk) 


*X. 


For  x  =  2  missiles  per  kill,  we  have  a  miss  rate  of 


1  —  2  P,  +K 


n  *- — — ^ 

2Pk-Pk 


[Eq.  B-5] 


[Eq.  B-6] 


[Eq.  B-7] 


In  addition  to  missile  counting,  we  also  add  the  capability  for  either  Red  or  Blue 
to  exit  the  combat  when  losses  reach  a  preset  level  by  assigning  minimum  values 
(min  b,  minr)  to  the  DEATFIIF  conditions.  If  minb  =  0  and  minr  =  0,  the 
model  corresponds  to  fighting  to  annihilation;  higher  values  will  generate  models 
reflecting  disengagement  breakpoints.  Listing  4  is  the  ASSIST  code  for  the  2  vs.  2 
case  with  missile  tracking. 

Listing  4.  ASSIST  Code  for  a  2  vs.  2  Engagement  with  Blue  Missile  Tracking 

(*  ASSIST  code  for  2  Blue  vs.  2  Red  with  6  missiles  per  Blue  aircraft  *) 

LIST  =3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

Time  =  100; 

(*  equipment  list  *) 

nblue  =2;  (*  blue  fighters  *) 

nred  =2;  (*  red  fighters  *) 

(*  kill  rates  *) 

k_rate_blue  =0.1;  (*  blue  kill  rate  *) 

k_rate_red  =  0.01;  (*  red  kill  rate  *) 

pk  =  0.85;  (*  missile  kill  probability  *) 

miss_rate  =  ( 1-pk) /pk*k_rate_blue;  (*  missile  miss  rate  *) 
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(*  minimum  number  of  combatants  for  state  pruning  *) 
min_b  =0;  (*  minimum  number  of  blue  aircraft  *) 
min_r  =0;  (*  minimum  number  of  red  aircraft  *) 

SPACE  = 

(nO : 0 . nblue, nl : 0 . nblue, n2 : 0 . nblue , n3 : 0 . nblue, n4 : 0 . nblue, n5 : 0 . nblue , n6 : 0 . nblue, 
red: 0 . nred) ; 

START  =  (0, 0, 0, 0, 0, 0, nblue, nred)  ; 

DEATH IF  (n0+nl+n2+n3+n4+n5+n6) <=min_b; ; 

DEATH IF  red  <=  min_r; 

(*  transition  cases:  Blue  kill  then  Red  kill*) 

IF  (n6>0)  AND  (red>0)  THEN 

TRANTO  n6=n6-l,  n5=n5+l,  red=red-l  BY  n6*k_rate_blue; 

TRANTO  n6=n6-l,  n5=n5+l  BY  n6*miss_rate; 

TRANTO  n6=n6-l  BY  (n6/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n5>0)  AND  (red>0)  THEN 

TRANTO  n5=n5-l,  n4=n4+l,  red=red-l  BY  n5*k_rate_blue; 

TRANTO  n5=n5-l,  n4=n4+l  BY  n5*miss_rate ; 

TRANTO  n5=n5-l  BY  (n5/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n4>0)  AND  (red>0)  THEN 

TRANTO  n4=n4-l ,  n3=n3+l,  red=red-l  BY  n4*k_rate_blue; 

TRANTO  n4=n4-l ,  n3=n3+l  BY  n4*miss_rate; 

TRANTO  n4=n4-l  BY  (n4/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n3>0)  AND  (red>0)  THEN 

TRANTO  n3=n3-l ,  n2=n2+l,  red=red-l  BY  n3*k_rate_blue; 

TRANTO  n3=n3-l,  n2=n2+l  BY  n3*miss_rate; 

TRANTO  n3=n3-l  BY  (n3/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n2>0)  AND  (red>0)  THEN 

TRANTO  n2=n2-l ,  nl=nl+l,  red=red-l  BY  n2*k_rateJolue; 

TRANTO  n2=n2-l ,  nl=nl+l  BY  n2*miss_rate; 

TRANTO  n2=n2-l  BY  (n2/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (nl>0)  AND  (red>0)  THEN 

TRANTO  nl=nl-l ,  n0=n0+l,  red=red-l  BY  nl*k_rate_blue; 

TRANTO  nl=nl-l ,  n0=n0+l  BY  nl*miss_rate; 

TRANTO  nl=nl-l  BY  (nl/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

(*  This  line  is  commented  out  based  on  Blues  ability  to  disengage  after  firing 
all  his  missiles  *) 

( *  I F  (red>0)  AND  (n0>0)  TRANTO  n0=n0-l  BY 
(n0/ (n0+nl+n2+n3+n4+n5+n6 ) ) *red*k_rate_red; * ) 


Adding  the  missile  Pk,  missile  miss  rate,  Blue  and  Red  aircraft  minimums,  the 
expanded  state  space,  and  the  starting  conditions  is  straightforward.  The  expan¬ 
sion  of  the  transfer  statements  deserves  some  discussion. 

The  basic  transition  logic  is  as  follows: 

♦  Separate  sets  of  transfers  are  established  for  each  Blue  missile  state. 

♦  State  transfers  occur  only  when  Blue  aircraft  are  in  the  state  and  Red  air¬ 
craft  are  available. 
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♦  Blue  kills  at  a  rate  proportional  to  the  total  number  of  Blue  aircraft.  When 
Blue  makes  a  kill,  one  aircraft  is  removed  from  the  Blue  state,  one  aircraft 
is  added  to  the  Blue  state  having  one  fewer  missiles,  and  one  aircraft  is 
removed  from  the  Red  state. 

♦  Blue  misses  at  a  rate  proportional  to  the  number  of  aircraft  in  the  Blue 
state.  When  Blue  misses,  one  aircraft  is  removed  from  the  Blue  state,  one 
aircraft  is  added  to  the  Blue  state  having  one  fewer  missiles,  but  no  aircraft 
are  removed  from  the  Red  state. 

♦  Red  kills  Blue  at  a  rate  proportional  to  the  number  of  Red  aircraft  and  the 
fraction  of  Blue  aircraft  in  the  Blue  state.  When  Red  makes  a  kill,  one  air¬ 
craft  is  removed  from  the  Blue  state. 

♦  Only  Red  kills  are  possible  when  Blue  has  no  missiles  left.  This  TRANTO 
statement  is  commented  out  in  the  example  listing  using  “(*”  and  “*)”  op¬ 
erators  based  on  the  assumption  that  Blue  can  disengage  at  will.  Another 
alternative  would  be  to  give  Blues  with  no  missiles  a  “guns  only”  kill  rate. 


Running  the  ASSIST  program  above  generates  a  model  having  103  states,  includ¬ 
ing  34  absorbing  (Death)  states,  and  301  transitions.  Listing  5  shows  the  first  and 
last  few  lines  of  the  model  file.  The  state  notation  in  the  model  shows  the  state 
number  and  the  contents  of  each  element  in  the  state.  This  model  has  eight  ele¬ 
ments  in  each  state.  The  first  seven  elements  are  Blue  aircraft  with  0  to  6  missiles, 
and  the  last  element  is  Red  aircraft.  Thus,  State  1  is  1(0, 0,0, 0,0, 0,2, 2),  indicating  2 
Blues  with  6  missiles  and  2  Reds.  The  first  transition  is  to  2(0, 0,0, 0,0, 1,1,1),  indi¬ 
cating  that  one  Blue  has  fired  a  missile  and  killed  one  Red.  A  Blue  absorbing  state 
corresponds  to  0s  in  all  of  the  first  seven  elements,  and  a  Red  absorbing  state  has 
0  in  the  eighth  element. 

Listing  5.  Model  File  Segments  for  2  vs.  2  Engagement 
with  Blue  Missile  Tracking 


LIST  =  3; 

NBLUE  =  2; 

NRED  =  2; 

K_RATE_BLUE  =  0.1; 

K_RATE_RED  =  0.01; 

PK  =  0.85; 

MISS_RATE  =  (1-PK) / PK*K_RATE_BLUE ; 

MIN_B  =  0; 

MIN  R  =  0; 


1  (* 

CM 

CM 

O 

O 

O 

O 

O 

O 

*) 

2  (* 

o 

o 

o 

o 

o 

*) 

=  2*K  RATE  BLUE; 

1  (* 

0,0, 0,0, 0,0, 2, 2 

*) 

3  ( * 

0,0, 0,0, 0,1, 1,2 

*) 

=  2*MISS  RATE; 

1  (* 
*2*: 

0,0, 0,0, 0,0, 2, 2 
K  RATE  RED; 

*) 

4  (* 

CM 

O 

O 

O 

O 

O 

O 

*) 

=  (2/ (0  +  0  +  0  +  0  +  0  +  0+2) 

2  (* 

0,0, 0,0, 0,1, 1,1 

*) 

5  (* 

0,0, 0,0, 0,2, 0,0 

*) 

=  1*K  RATE  BLUE; 

2  (* 

t — 1 

O 

o 

o 

o 

o 

*) 

6  (* 

o 

CM 

O 

O 

O 

O 

O 

*) 

=  1*MISS  RATE; 

2  (* 
*  *■ 

0,0, 0,0, 0,1, 1,1 
K  RATE  RED; 

*) 

7  (* 

O 

O 

o 

o 

o 

o 

*) 

=  (1/ (0  +  0  +  0  +  0  +  0  +  1  +  1) 

2  (* 

0,0, 0,0, 0,1, 1,1 

*) 

8  (* 

0,0, 0,0, 1,0, 1,0 

*) 

=  1*K  RATE  BLUE; 

2  (* 

o 

o 

o 

o 

o 

*) 

9  (* 

o 

o 

o 

o 

o 

*) 

=  1*MISS  RATE; 

2  (* 

0,0, 0,0, 0,1, 1,1 

*) 

10  r 

*  0,0, 0,0, 0,0, 1,1 

*■  ’ 

)  =  (1/ (0+0+0+0+0+1+1 
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*1*K_RATE_RED; 

3 ( *  0,0, 0,0, 0,1, 1,2  *) 
3 ( *  0,0, 0,0, 0,1, 1,2  *) 
3 ( *  0,0, 0,0, 0,1, 1,2  *) 
*2*K_RATE_RED; 

3 ( *  0,0, 0,0, 0,1, 1,2  *) 
3 ( *  0,0, 0,0, 0,1, 1,2  *) 
3 ( *  0,0, 0,0, 0,1, 1,2  *) 


6 ( *  0,0, 0,0, 0,2,0, 1  *) 
11 (*  0,0, 0,0, 0,2, 0,2  *) 
12 (*  0,0, 0,0, 0,1, 0,2  *) 

9 ( *  0,0, 0,0, 1,0, 1,1  *) 
13 (*  0,0, 0,0, 1,0, 1,2  *) 
4 (*  0,0, 0,0, 0,0, 1,2  *) 


=  1*K_RATE_BLUE; 

=  1*MISS_RATE; 

=  (1/ (0+0+0+0+0+1+1) ) 

=  1*K_RATE_BLUE; 

=  1*MISS_RATE; 

=  (1/ (0+0+0+0+0+1+1) ) 


97 (*  1,0, 1,0, 0,0, 0,2  *) , 
97 (*  1,0, 1,0, 0,0, 0,2  *) , 
97 (*  1,0, 1,0, 0,0, 0,2  *) , 
*2*K_RATE_RED; 

99 (*  1,  1,0, 0,0, 0,0, 1  *)  , 
99 (*  1, 1,0, 0,0, 0,0, 1  *)  , 
99 (*  1, 1,0, 0,0, 0,0, 1  *) , 
*1*K_RATE_RED; 

100 (*  1, 1,0, 0,0, 0,0, 2  *) 
100 ( *  1,1, 0,0, 0,0, 0,2  *) 
100 ( *  1,1, 0,0, 0,0, 0,2  *) 
*2*K  RATE  RED; 


99  ( *  1,1, 0,0, 0,0, 0,1  *) 
100 (*  1,1, 0,0, 0,0, 0,2  *) 
76 (*  1,0, 0,0, 0,0, 0,2  *) 

101  (*  2, 0,0, 0,0, 0,0,0  *) 
102  (*  2, 0, 0, 0, 0, 0, 0,  1  *) 
72 (*  1,0, 0,0, 0,0, 0,1  *) 

102  (*  2, 0, 0, 0, 0, 0, 0, 1  * 
103 (*  2, 0,0, 0,0, 0,0, 2  * 
76 (*  1,0, 0,0, 0,0, 0,2  *) 


=  1*K_RATE_BLUE; 

=  1*MISS_RATE; 

=  (1/ (1+0+1+0+0+0+0) ) 

=  1*K_RATE_BLUE; 

=  1*MISS_RATE; 

=  (1/ (1+1+0+0+0+0+0) ) 

)  =  1*K_RATE_BLUE; 

)  =  1*MISS_RATE; 

=  (1/ (1+1+0+0+0+0+0) ) 


(*  NUMBER  OF  STATES  IN  MODEL  =  103  *) 

(*  NUMBER  OF  TRANSITIONS  IN  MODEL  =  285  *) 


At  100  time  units,  the  engagement  is  0.99994  complete;  the  sums  of  the  absorbing 
missile  tracking  states  agree  within  four  significant  figures  with  the  sums  of  the 
absorbing  states  for  the  simple  2  vs.  2  example.  Figure  B-5  combines  the  prob¬ 
abilities  for  the  absorbing  states  to  show  the  marginal  probabilities  of  Red  victory 
and  the  marginal  probabilities  of  Blue  victory  with  specific  numbers  of  aircraft 
and  missiles  after  a  2  vs.  2  engagement  of  100  time  units. 

Figure  B-5.  2  vs.  2  Engagement  Aircraft  and  Missile  Configuration  Probabilities 

Marginal  Probabilities  of  Red  Win  and  of  Blue  Wins  for  N  Blue  with  X  Missiles  Remaining 
(Time=100,  Kb=0.1,  Kr=0.01,  Missile  Pk=0.85) 
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Example  3 — 4  vs.  4+4+4  Sequential  Engagement 

The  standard  SLAACM  campaign  scenario  has  4  Blue  defenders  engaging  4  Red 
long-range  escorts  (LEs),  4  close  escorts  (CEs),  and  4  bombers  sequentially.  This 
sequential  engagement  can  be  modeled  using  ASSIST  by  imposing  on  the 
TRANTO  statements  the  conditions  that  Blue  cannot  engage  the  close  escorts  or 
bombers  until  all  the  LEs  are  dead  and  cannot  engage  the  bombers  until  all  the 
CEs  are  also  dead.  We  include  individual  kill  rates  for  the  three  pairs  of  combat¬ 
ants.  Listing  6  is  the  ASSIST  code  for  the  sequential  engagement.  The  resulting 
model  has  38  states  and  48  transitions  when  the  Blue  breakpoint  is  2  losses 
(minb  =  2  in  the  listing),  and  has  64  states  and  96  transitions  when  Blue  fights  to 
annihilation  (min  b  =  0).  Figure  B-6  shows  the  marginal  probabilities  for  the  out¬ 
comes  for  the  2  Blue  breakpoint,  and  Figure  B-7  shows  the  corresponding  ex¬ 
pected  number  of  surviving  aircraft  types. 

Listing  6.  4  Blue  Sequentially  Engaging  4  Red  LE,  4  Red  CE, 
and  4  Red  Bombers 


(*  ASSIST  model  for  Fighter  Combat  *) 

(*  This  is  multistage  combat  where  blue  aircraft  sequentially  engage*) 
(*  three  classes  of  red  aircraft  *) 

(*  This  version  does  not  include  missiles  *) 

LIST  =3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

(*"Time  =  5  to  100  By  5;"*) 

(*  the  time  is  set  here  -  no  matter  what  the  STEM  GUI  shows  *) 

Time  =  1000; 

(*  equipment  list  *) 

nblue  =4;  (*  blue  fighters  *) 

nredl  =4;  (*  red  fighters  -  lead  escorts  *) 

nred2  =4;  (*  red  fighters  -  close  escorts  *) 

nred3  =4;  (*  red  bombers  or  fighter  bombers  *) 

min_b  =2; 
min_rl  =  0; 
min_r2  =  0; 
min  r3  =  0; 


(*  kill  probabilities  *) 

(*  Note:  the  1st  12  characters  of  a  variable  must  be  unique  *) 
blue_redl_k_rate  =  0.1;  (*  0.1  *) 

blue_red2_k_rate  =  0.1;  (*  0.2  *) 
blue_red3_k_rate  =  0.1;  (*  0.4  *) 

k_rate_redl  =  0.02;  (*  0.05  *) 

k_rate_red2  =  0.01;  (*  0.01  *) 

k_rate_red3  =  0.005;  (*  0.005  *) 

SPACE  =  (blue:  0. nblue,  redl :  0. nredl,  red2 :  0.nred2,  red3 :  0.nred3); 

START  =  (nblue,  nredl,  nred2,  nred3) ; 

DEATH IF  (blue  =  min_b)  AND  ( (redl>0)  OR  (red2>0)  OR  (red3>0) ) ; 

DEATH IF  (redl  =  min_rl)  AND  (red2  =  min_r2)  AND  (red3  =  min_r3)  AND  (blue>0) ; 

(*  Phase  1  combat  *) 

(*  Blue  kills  *) 
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IF  (blue  >  min_b)  and  (redl  >  min_rl)  TRANTO  redl  =  redl-1  BY 
blue*blue_redl_k_rate; 

(*  Red  kills  *) 

IF  (redl  >  min_rl  )  and  (blue>min_b)  TRANTO  blue  =  blue-1  BY  redl*k_rate_redl ; 

(*  Phase  2  combat  *) 

(*  Blue  kills  *) 

IF  (blue  >  min_b)  and  (redl=min_rl )  and  (red2  >  min_r2)  TRANTO  red2  =  red2-l 
BY  blue*blue_red2_k_rate; 

(*  Red  kills  *) 

IF  (red2  >  min_r2  )  and  (redl=min_rl )  and  (blue>min_b)  TRANTO  blue  =  blue-1  BY 
red2  *  k_rate_red2 ; 

(*  Phase  3  combat  *) 

(*  Blue  kills  *) 

IF  (blue  >  min_b)  and  (red3  >  min_r3)  TRANTO  red3  =  red3-l  BY 
blue*blue_red3_k_rate; 

(*  Red  kills  *) 

IF  (red3  >  min_r3  )  and  (redl=min_rl )  and  (red2=min_r2 )  and  (blue>min_b) 

TRANTO  blue  =  blue-1  BY  red3*k_rate_red3 ; 

Figure  B-6.  Sequential  Engagement  Marginal  Probabilities 

Aircraft  Count  per  Type  Probabilities 
4B  x  4R1 ,  then  x  4R2  then  x  4R3 
(B:R  Kill  Rates  -  5:1,  10:1,  20:1,  Blue  Min  =  2) 
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Figure  B-7.  Sequential  Engagement  Surviving  Aircraft 


Post  Engagement  Aircraft 
4B  x  4R1,  then  x  4R2  then  x  4R3 
(B:R  Kill  Rates  -  5:1,  10:1,  20:1,  Blue  Min  =  2) 


3.50 
3.00 

2.50 
2.00 

1.50 
1.00 
0.50 
0.00 


3.17- 


1.01 

_ 0. 77 _ _ 


Blue  Red  1  Red  2  Red  3 


B-13 


Example  4 — 4  vs.  4+4+4  Sequential  Engagement  with  Missiles 

The  last  example  is  the  standard  sequential  engagement  with  missiles  on  the  Blue 
aircraft.  Listing  7  is  the  ASSIST  code  for  this  case.  This  listing  generates  a  model 
having  3,695  states  and  21,729  transitions  for  the  2  Blue  breakpoint  case  and 
3,823  states  and  23,454  transitions  for  the  Blue  fight  to  annihilation  case.  Figure 
B-8  shows  the  marginal  probabilities  of  missile  counts  combined  with  information 
on  how  many  aircraft  are  carrying  the  missiles  for  the  case  in  which  Blue  breaks 
after  two  losses.  Figure  B-9  shows  the  remaining  aircraft  count  by  type.  Table  B-l 
contains  the  summary  results  for  both  the  2  Blue  breakpoint  and  Blue  fight  to  an¬ 
nihilation  cases. 

Listing  7.  4  Blue  vs.  4  Red  LE,  4  Red  CE,  and  4  Bombers  Sequentially 

with  Blue  Missiles 


(*  ASSIST  model  for  *) 

(*  This  version  includes  sequential  combat  with  3  types  of  red  fighters  *) 
(*  and  includes  Blue  missiles  *) 

(*  Fighter  Combat  using  Dave  Lee's  construct  for  blue  missile  counting*) 

(*  7/20/05  This  version  includes  missiles  for  blues  *) 

(*  7/20/05  THIS  VERSION  INCLUDES  6  BLUE  MISSILES 

(*  THIS  VERSION  CAN  LIMIT  STATES  BASED  ON  A  MINIMUM  NUMBER  OF  COMBATANTS  *) 
(*  The  kill  rate  is  for  one-on-one  and  is  multiplied  for  N  killers  *) 

(*  Assume  one  missile  per  kill  even  with  N  killers  *) 

(*  Therefore,  the  kill  rate  includes  one  hit  and  all  misses  *) 

(*  Find  miss  rate  based  on  ratio  of  Pmiss:Pk  *) 

LIST  =3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

(*"Time=10  TO+  100  BY  10;"*) 

Time  =  1000; 

(*  equipment  list  *) 
nblue  =4;  (*  blue  fighters  *) 

nredl  =4;  (*  red  fighters  *) 

nred2  =  4 ; 
nred3  =  4 ; 


(*  missile  equipage  *) 
b_msls  =6;  (*  blue  missiles  *) 

(*  the  first  12  characters  in  a  variable  name  are  significant  and  must  be 
unique  *) 

(*  kill  rates  *) 

redl_k_rate_blue  =  0.1;  (*  B:R1  =  5:1  *) 

red2_k_rate_blue  =  0.1;  (*  B:R2  =  10:1  *) 

red3_k_rate_blue  =  0.1;  (*  B:R3  =  20:1  *) 

k_rate_redl  =  0.02; 
k_rate__red2  =  0.01; 
k  rate  red3  =  0.005; 


(*  blue  missiles  miss  rates  *) 
bpk  =  0.85;  (*  missile  kill  probability  *) 

redl_bmiss_rate  =  ( 1-bpk) /bpk*redl_k_rate_blue;  (*  missile  miss  rate  against 
redl  *) 

red2_bmiss_rate  =  ( 1-bpk) /bpk*red2_k_rate_blue;  (*  missile  miss  rate  against 
red2  *) 
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red3_bmiss_rate  =  ( 1-bpk) /bpk*red3_k_rate_blue;  (*  missile  miss  rate  against 
red3  *) 

(*  minimum  number  of  combatants  for  state  pruning  *) 
min_b  =2;  (*  minimum  number  of  blue  aircraft  *) 

min_rl  =0;  (*  minimum  number  of  redl  aircraft  *) 

min_r2  =  0; 
min_r3  =  0; 

SPACE  =  (b : array [0 .b_msls] of  O.nblue,  redl:  O.nredl,  red2 :  0.nred2,  red3 : 

0 . nred3) ; 

(*  states  0  to  max  missiles  *) 

(*  THIS  STATEMENT  MUST  BE  CHANGED  WHEN  THE  MISSILE  COUNT  IS  CHANGED!  *) 

START  =  (0, 0, 0, 0, 0, 0, nblue, nredl, nred2, nred3)  ; 

DEATH IF  (sum (b) <=  min_b)  OR  (b [ 0 ] =4 ) ; 

DEATH IF  (redl  =  min_rl)  AND  (red2  =  min_r2)  AND  (red3  =  min_r3) ; 

(*  state  transitions  and  kills  *) 

(*  advanced  escorts  *) 

FOR  I  IN  [l.b_msls] 

IF  (b[I]>0)  AND  (sum (b) >min_b)  AND  ( redl>min_rl )  THEN 

TRANTO  b [ I ] =b [ I ] - 1 ,  b [1-1] =b [1-1] +1,  redl  =redl-l  BY  b [ I ] *redl_k_rate_blue; 
TRANTO  b [ I ] =b [ I ] -1 ,  b [ I— 1 ] =b [1-1] +1  BY  b [ I ] *redl_bmiss_rate; 

TRANTO  b [ I ] =b [ I ] — 1  BY  (b [ I ] /sum (b) ) *redl *k_rate_redl ; 

ENDIF; 

ENDFOR; 

(*  no  missile  case  against  advanced  escorts  *) 

IF  (b [ 0 ] >0 )  AND  (sum (b) >min_b)  AND  (redl  >  min_rl)  THEN 

(*  TRANTO  redl=redl-l  BY  b [ 0 ] *redl_k_rate_blue;  *)  (*  NO  BLUE  KILLS  WITHOUT 

MISSILES  *) 

(*  TRANTO  b [ 0 ] =b [ 0 ] — 1  BY  (b [ 0 ] /sum (b) ) *redl *k_rate_redl ;  *)  (*  BLUE  ESCAPES  IF 

NO  MISSILES  *) 

ENDIF; 

(*  close  escorts  *) 

IF  redl  =  min_rl  THEN 
FOR  I  IN  [l.b_msls] 

IF  (b [I] >0)  AND  (sum (b) >min_b)  AND  (red2  >  min_r2)  THEN 

TRANTO  b [ I ] =b [ I ] -1 ,  b[I-l]=b[I-l]+l,  red2  =red2-l  BY  b [ I ] *red2_k_rate_blue; 
TRANTO  b [ I ] =b [ I ] -1 ,  b [ I— 1 ] =b [1-1] +1  BY  b [ I ] *red2_bmiss_rate; 

TRANTO  b [ I ] =b [ I ] — 1  BY  (b [ I ] /sum (b) ) *red2 *k_rate_red2 ; 

ENDIF; 

ENDFOR; 

(*  no  missile  case  against  close  escorts  *) 

IF  (b [0] >0)  AND  (sum (b) >min_b)  AND  (red2  >  min_r2)  THEN 

(*  TRANTO  red2=red2-l  BY  b [ 0 ] *red2_k_rate_blue;  *)  (*  NO  BLUE  KILLS  WITHOUT 

MISSILES  *) 

(*  TRANTO  b [ 0 ] =b [ 0 ] -1  BY  (b [ 0 ] /sum (b) ) *red2 *k_rate_red2 ;  *)  (*  BLUE  ESCAPES  IF 

NO  MISSILES  *) 

ENDIF; 

ENDIF; 

(*  bombers  *) 

IF  (red2  =  min_r2)  AND  (redl  =  min_rl)  THEN 
FOR  I  IN  [l.b_msls] 

IF  (b [I] >0)  AND  (sum (b) >min_b)  AND  (red3  >  min_r3)  THEN 

TRANTO  b [ I ] =b [ I ] -1 ,  b [ I  — 1 ] =b [1-1] +  1,  red3  =red3-l  BY  b [ I ] * red3_k_rate_blue ; 
TRANTO  b [ I ] =b [ I ] -1 ,  b [ I — 1 ] =b [ I — 1 ] + 1  BY  b [ I ] *red3_bmiss_rate; 

TRANTO  b [ I ] =b [ I ] — 1  BY  (b [ I ] / sum (b) ) *red3*k_rate_red3 ; 

ENDIF; 

ENDFOR; 

(*  no  missiles  case  against  bombers  *) 

IF  (b [0] >0)  AND  (sum (b) >min_b)  AND  (red3  >  min_r3)THEN 

(*  TRANTO  red3=red3-l  BY  b [ 0 ] *red3_k_rate_blue;  *)  (*  NO  BLUE  KILLS  WITHOUT 

MISSILES  *) 

(*  TRANTO  b [0] =b [0] -1  BY  (b [ 0 ] /sum (b) ) *red3*k_rate_red3 ;  *)  (*  BLUES  ESCAPES 

IF  NO  MISSILES  *) 

ENDIF; 

ENDIF; 
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Marginal  Probability 


(*  FOR  NO  MISSILE  CASES:  *) 

(*  use  only  red  rates  if  blue  has  no  missiles  left  and  cannot  escape  *) 
(* —  comment  out  all  statements  if  blue  can  escape  *) 

(* —  include  blue  kill  statements  to  continue  combat  without  missiles  *) 
(*  possibly  add  rates  above  and  TRANTO  statments  here  for  Blue  guns  *) 


Figure  B-8.  Sequential  Engagement  Missile  Marginal  Probabilities 
with  Carrier  Aircraft  Information 


4v12  Sequential  Engagement  Aircraft  and  Missiles  Remaining 
(BvLE-5:1,  BvCE-10:1,  BvBm-20:1,  Min  Blue  =  2,  T=1000  units) 
(Blue  with  0  missiles  escape) 


24  23  22  21  20  19  18  17  16  15  14  13  12  11  10  9  8 

Missiles  Remaining  (Count) 


□  2  Blue 
■  3  Blue 

□  4  Blue 


Figure  B-9.  Sequential  Engagement  with  Missiles,  Remaining  Aircraft 


4v12  Sequential  Engagement  Aircraft  Remaining 
(BvLE-5:1,  BvCE-10:1,  BvBm-20:1,  Min  Blue  =  2,  T=1000  units) 
(Blue  with  0  missiles  escape) 
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Analyzing  Large  State-Space  Engagement  Models  with  ASSIST  and  STEM 


Table  B-1.  Summary  Engagement  Results 


Scenario 

Blue 

remaining 

Red 

remaining 

Missiles 

used 

2  Blue  breakpoint 

3.2 

2.1 

16 

Blue  fight  to  annihilation 

3.0 

0.8 

18 

Summary 


In  this  appendix,  we  have  shown  how  the  NASA-developed  Markov  tools 
ASSIST  and  STEM  can  be  used  for  analysis  of  increasingly  complex  air  combat 
engagements.  We  have  shown  through  examples  how  ASSIST  is  particularly  use¬ 
ful  for  supporting  parametric  analysis  involving  large  state-space  problems.  The 
NASA  tools  have  proved  useful  for  standalone  analyses,  for  prototyping  of 
SLAACM  engagements,  and  for  independent  confirmation  of  SLAAM  calcula¬ 
tions. 
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